Search PubMed⌕ Search

Biomedical subjects

Henry F Schaefer

Publications and source records attributed to Henry F Schaefer.

At least 73 records · Page 4Linked to original sources
PubMed2004

Hypervalency avoided: simple substituted BrF3 and BrF5 molecules. Structures, thermochemistry, and electron affinities of the bromine hydrogen fluorides HBrF2 and HBrF4.

Five different pure density functional theory (DFT) and hybrid Hartree-Fock/DFT methods have been used to search for the molecular structures, thermochemistry, and electron affinities of the bromine hydrogen fluorides HBrF(n)/HBrF(n)(-) (n = 2, 4). The basis sets used in this work are of double-zeta plus polarization quality in conjunction with s- and p-type diffuse functions, labeled as DZP++. Structures with Br-F and Br-H normal bonds, that is, HBrF(2)/HBrF(2)(-) with C(2v) or C(s) symmetry and HBrF(4)/HBrF(4)(-) with C(4v) or C(s) symmetry, are genuine minima. However, unlike the original BrF(3) and BrF(5) molecules, the global minima for HBrF(n)/HBrF(n)(-) (n = 2, 4) species are predicted to be complexes, some of which contain hydrogen bonds. The demise of the hypervalent structures is due to the availability of favorable dissociation products involving HF, which has a much larger dissociation energy than F(2). Similar reasoning suggests that PF(4)H, SF(3)H, SF(5)H, ClF(2)H, ClF(4)H, AsF(4)H, SeF(3)H, and SeF(5)H will all be hydrogen bond structures incorporating diatomic HF. The most reasonable theoretical values of the adiabatic electron affinities (EA(ad)) are 3.69 (HBrF(2)) and 4.38 eV (HBrF(4)) with the BHLYP method. These electron affinities are comparable to those of the analogous molecules: Br(2)F(n), ClBrF(n), and BrF(n)(+1) systems. The first F-atom dissociation energies for the neutral global minima are 60 (HBrF(2)) and 49 kcal/mol (HBrF(4)) with the B3LYP method. The first H-atom dissociation energies for the same systems are 109 (HBrF(2)) and 116 kcal/mol (HBrF(4)). The large Br-H bond energies are not sufficient to render the hypervalent structures energetically tenable. The dissociation energies for the complexes to their fragments are relatively small.

Journal Article↗
PubMed2004

The ability of silylenes to bind excess electrons: electron affinities of SiX(2), and SiXY species (X,Y=H,CH(3),SiH(3),F,Cl,Br).

Recently, Ishida and co-workers have isolated silylene radical anions via the one-electron reduction of isolable cyclic dialkylsilylenes, discovering these corresponding radical anions to be relatively stable at low temperatures. Herein we report theoretical predictions of the adiabatic electron affinities (AEA), vertical electron affinities, and vertical detachment energies of a series of methyl, silyl, and halosubstituted silylene compounds. This research utilizes the carefully calibrated DZP++ basis with the combination of the popular nonhybrid and hybrid DFT functionals, BLYP, B3LYP, and BHHLYP. The level of theory employed and the ensemble of species under study confirm the ability of silylenes to bind excess electrons with Si(SiH(3))(2) being the most effective, having a predicted AEA of 1.95 eV. While it is known that methyl substituents have a diminishing effect on the computed electron affinities (EAs), it is shown that fluorine shows an analogous negative effect. Similarly, previous suggestions that Si(CH(3))(2) will not bind an electron appear incorrect, with EA[Si(CH(3))(2)] predicted here to be 0.46 eV.

Journal Article↗
PubMed2004

Dihydrogen and methane elimination from adducts formed by the interaction of carbenium and silylium cations with nucleophiles.

Stationary points for reactions R'R' 'HX(+) + YH --> [R'R' 'X-Y](+) + H(2) (I) and R'(CH(3))HX(+) + YH -->[R'HX-Y](+) + CH(4) (II) (R', R' ' = CH(3), H; X = C, Si; Y = CH(3)O, (CH(3))(2)N, and C(6)H(5)) are located and optimized by the B3LYP/aug-cc-pVDZ method. A similar mechanism was found to be operative for both types of reactions with X = C and X = Si. Formation of the intermediate (adduct) results in the transfer of electron density from the electron-rich bases to the X atoms and in the growth of a positive charge on a hydrogen atom attached to Y. This mobile proton may shift from Y to X, and the relative energies of transition states for elimination reactions (Delta) depend on the ability of the X atom to retain this proton. Therefore, Deltagrows on going from Si to C and with increasing numbers of methyl substituents. For X = C, the Deltavalue for both reactions correlates well with the population of the valence orbitals of X in a wide range from -44 kcal/mol (methyl cation/benzene) to 31 kcal/mol (isopropyl cation/methanol). For X = Si this range is more narrow (from -19 to -5.0 kcal/mol), but all Delta values are negative with the exclusion of silylium ion/benzene systems, adducts of which are pi- rather than sigma-complexes. The energy minima for product complexes for H(2) elimination are very shallow, and several are dissociative. However, complexes with methane which exhibit bonding between X and the methane hydrogen are substantially stronger, especially for systems with X = Si. The latter association energy may reach 8 kcal/mol (Si...H distance is 2 A).

Journal Article↗
PubMed2004

Thermochemistry of disputed soot formation intermediates C4H3 and C4H5.

Accurate isomeric energy differences and standard enthalpies of formation for disputed intermediates in soot formation, C(4)H(3) and C(4)H(5), have been determined through systematic extrapolations of ab initio energies. Electron correlation has been included through second-order Z-averaged perturbation theory (ZAPT2), and spin-restricted, open-shell coupled-cluster methods through triple excitations [ROCCSD, ROCCSD(T), and ROCCSDT] utilizing the correlation-consistent hierarchy of basis sets, cc-pVXZ (X = D, T, Q, 5, and 6), followed by extrapolations to the complete basis set limit via the focal point method of Allen and co-workers. Reference geometries were fully optimized at the ROCCSD(T) level with a TZ(2d1f,2p1d) basis set. Our analysis finds that the resonance-stabilized i-C(4)H(3) and i-C(4)H(5) isomers lie 11.8 and 10.7 kcal mol(-1) below E-n-C(4)H(3) and E-n-C(4)H(5), respectively, several kcal mol(-1) (more, less) than reported in recent (diffusion Monte Carlo, B3LYP density-functional) studies. Moreover, in these systems Gaussian-3 (G3) theory suffers from large spin contamination in electronic wave functions, poor reference geometries, and anomalous vibrational frequencies, but fortuitous cancellation of these sizable errors leads to isomerization energies apparently accurate to 1 kcal mol(-1). Using focal-point extrapolations for isodesmic reactions, we determine the enthalpies of formation (delta(f)H(0) (composite function)) for i-C(4)H(3), Z-n-C(4)H(3), E-n-C(4)H(3), i-C(4)H(5), Z-n-C(4)H(5), and E-n-C(4)H(5) to be 119.0, 130.8, 130.8, 78.4, 89.7, and 89.1 kcal mol(-1), respectively. These definitive values remove any remaining uncertainty surrounding the thermochemistry of these isomers in combustion models, allowing for better assessment of whether even-carbon pathways contribute to soot formation.

Journal Article↗
PubMed2004

Spontaneous gas-phase generation of needle-shaped clusters which violate the isolated square rule: a facile road to GaN nanorods?

Structures and thermodynamic properties of the imidogallanes [HGaNH](n) with high oligomerization degrees n = 7-16 and related amido-imido compounds have been investigated theoretically at the B3LYP level of theory with all-electron pVDZ and effective core potential LANL2DZ(d,p) basis sets. Needle-shaped oligomers which violate the "isolated square rule" were found to be more stable than cage isomers. The needle-shape oligomer with n = 16 is predicted to be exceptionally stable at low temperatures. Octamer and hexamer clusters dominate the gas phase at higher temperatures. The highest oligomerization degree of the spontaneous cluster formation has been estimated. It is concluded that generation of the gas-phase [HGaNH](n) clusters with oligomerization degree n >or= 60 is viable. This makes these species possible intermediates involved in the gas-phase generation of GaN nanoparticles. A case study of methyl-substituted analogues suggests that formation of the gaseous [MeGaNH](n) oligomers is even more favorable compared to that of [HGaNH](n). We predict that spontaneous growth of GaN oligomers is favorable thermodynamically. Laser-assisted generation of GaN nanorods at low-temperature conditions appears to be feasible.

Journal Article↗
PubMed2004

The low-lying electronic excited states of NiCO.

Highly correlated coupled cluster methods with single and double excitations (CSSD) and CCSD with perturbative triple excitations were used to predict molecular structures and harmonic vibrational frequencies for the electronic ground state X 1Sigma+, and for the 3Delta, 3Sigma+, 3Phi, 1 3Pi, 2 3Pi, 1Sigma+, 1Delta, and 1Pi excited states of NiCO. The X 1Sigma+ ground state's geometry is for the first time compared with the recently determined experimental structure. The adiabatic excitation energies, vertical excitation energies, and dissociation energies of these excited states are predicted. The importance of pi and sigma bonding for the Ni-C bond is discussed based on the structures of excited states.

Journal Article↗
PubMed2004

Toward subchemical accuracy in computational thermochemistry: focal point analysis of the heat of formation of NCO and [H,N,C,O] isomers.

In continuing pursuit of thermochemical accuracy to the level of 0.1 kcal mol(-1), the heats of formation of NCO, HNCO, HOCN, HCNO, and HONC have been rigorously determined using state-of-the-art ab initio electronic structure theory, including conventional coupled cluster methods [coupled cluster singles and doubles (CCSD), CCSD with perturbative triples (CCSD(T)), and full coupled cluster through triple excitations (CCSDT)] with large basis sets, conjoined in cases with explicitly correlated MP2-R12/A computations. Limits of valence and all-electron correlation energies were extrapolated via focal point analysis using correlation consistent basis sets of the form cc-pVXZ (X=2-6) and cc-pCVXZ (X=2-5), respectively. In order to reach subchemical accuracy targets, core correlation, spin-orbit coupling, special relativity, the diagonal Born-Oppenheimer correction, and anharmonicity in zero-point vibrational energies were accounted for. Various coupled cluster schemes for partially including connected quadruple excitations were also explored, although none of these approaches gave reliable improvements over CCSDT theory. Based on numerous, independent thermochemical paths, each designed to balance residual ab initio errors, our final proposals are DeltaH(f,0) ( composite function )(NCO)=+30.5, DeltaH(f,0) ( composite function )(HNCO)=-27.6, DeltaH(f,0) ( composite function )(HOCN)=-3.1, DeltaH(f,0) ( composite function )(HCNO)=+40.9, and DeltaH(f,0) ( composite function )(HONC)=+56.3 kcal mol(-1). The internal consistency and convergence behavior of the data suggests accuracies of +/-0.2 kcal mol(-1) in these predictions, except perhaps in the HCNO case. However, the possibility of somewhat larger systematic errors cannot be excluded, and the need for CCSDTQ [full coupled cluster through quadruple excitations] computations to eliminate remaining uncertainties is apparent.

Journal Article↗
PubMed2004

Cyclic perfluorocarbon radicals and anions having high global warming potentials (GWPs): structures, electron affinities, and vibrational frequencies.

Adiabatic electron affinities, optimized molecular geometries, and IR-active vibrational frequencies have been predicted for small cyclic hydrocarbon radicals C(n)H(2)(n)(-)(1) (n = 3-6) and their perfluoro counterparts C(n)F(2)(n)(-)(1) (n = 3-6). Total energies and optimized geometries of the radicals and corresponding anions have been obtained using carefully calibrated (Chem. Rev. 2002, 102, 231) density functional methods, namely, the B3LYP, BLYP, and BP86 functionals in conjunction with the DZP++ basis set. The predicted electron affinities show that only the cyclopropyl radical tends to bind electrons among the hydrocarbon radicals studied. The trend for the perfluorocarbon (PFC) radicals is quite different. The electron affinities increase with expanding ring size until n = 5 and then slightly decrease at n = 6. Predicted electron affinities of the hydrocarbon radicals using the B3LYP hybrid functional are 0.24 eV (C(3)H(5)/C(3)H(5)(-)), -0.19 eV (C(4)H(7)/C(4)H(7)(-)), -0.15 eV (C(5)H(9)/C(5)H(9)(-)), and -0.11 eV (C(6)H(11)/C(6)H(11)(-)). Analogous electron affinities of the perflurocarbon radicals are 2.81 eV (C(3)F(5)/C(3)F(5)(-)), 3.18 eV (C(4)F(7)/C(4)F(7)(-)), 3.34 eV (C(5)F(9)/C(5)F(9)(-)), and 3.21 eV (C(6)F(11)/C(6)F(11)(-)).

Journal Article↗
PubMed2004

The diazocarbene (CNN) molecule: characterization of the X 3Sigma- and A 3Pi electronic states.

The ground (X (3)Sigma(-)) and first excited triplet (A (3)Pi) electronic states of diazocarbene (CNN) have been investigated systematically starting from the self-consistent-field theory and proceeding to the coupled cluster with single, double, and full triple excitations (CCSDT) method with a wide range of basis sets. While the linear X (3)Sigma(-) ground state of CNN has a real degenerate bending vibrational frequency, the A (3)Pi state of CNN is subject to the Renner-Teller effect and presents two distinct real vibrational frequencies along the bending coordinate. The bending vibrational frequencies of the A (3)Pi state were evaluated via the equation-of-motion coupled cluster (EOM-CC) techniques. The significant sensitivity to level of theory in predicting the ground-state geometry, harmonic vibrational frequencies, and associated infrared intensities has been attributed to the fact that the reference wave function is strongly perturbed by the excitations of 1pi-->3pi followed by a spin flip. At the highest level of theory with the largest basis set, correlation-consistent polarized valence quadruple zeta (cc-pVQZ) CCSDT, the classical X-A splitting (T(e) value) was predicted to be 68.5 kcal/mol (2.97 eV, 24 000 cm(-1)) and the quantum mechanical splitting (T(0) value) to be 69.7 kcal/mol (3.02 eV, 24 400 cm(-1)), which are in excellent agreement with the experimental T(0) values, 67.5-68.2 kcal/mol (2.93-2.96 eV, 23 600-23 900 cm(-1)). With the EOM-CCSD method the Renner parameter (epsilon) and averaged bending vibrational frequency (omega(2)) for the A (3)Pi state were evaluated to be epsilon=-0.118 and omega(2)=615 cm(-1), respectively. They are in fair agreement with the experimental values of epsilon=-0.07 and nu(2)=525 cm(-1).

Journal Article↗
PubMed2004

Structural and thermodynamic properties of group 13 imidometallanes and their heavier analogues.

Systematic theoretical studies of the [XMYH](n) inorganic rings and clusters (M = Al, Ga, In; Y = N, P, As; X = H, F, Cl, Br, I; n = 1-6) have been carried out using hybrid Hartree-Fock density functional theory. A consistent set of the structural and thermodynamic properties has been obtained. The stability of the MY bond decreases in the order Al > Ga >or= In; N >> P > As. Terminal groups X have a minor influence on the subsequent elimination enthalpies of the clusters. In the case of X = H, hydrogen elimination makes formation of the [HMYH](6) oligomers from MH(3) and YH(3) thermodynamically favorable; while in the case of halide substituents, formation of [XMYH](6) is thermodynamically unfavorable, except for the system with the strongest MY bond (AlN). Substitution of the acidic hydrogen by X is favorable energetically for all [HMYH](6) clusters, but is complicated by the processes of cluster destruction to form the [X(2)MYH(2)](2) dimers. The high stability of the [HMNH](6) clusters makes them attractive single-source precursors for the production of 13-15 composites.

Journal Article↗
PubMed2004

Low-lying electronic states of FeNC and FeCN: a theoretical journey into isomerization and quartet/sextet competition.

With several levels of multireference and restricted open-shell single-reference electronic structure theory, optimum structures, relative energetics, and spectroscopic properties of the low-lying (6)Delta, (6)Pi, (4)Delta, (4)Pi, and (4)Sigma(-) states of linear FeNC and FeCN have been investigated using five contracted Gaussian basis sets ranging from Fe[10s8p3d], C/N[4s2p1d] to Fe[6s8p6d3f2g1h], C/N[6s5p4d3f2g]. Based on multireference configuration interaction (MRCISD+Q) results with a correlation-consistent polarized valence quadruple-zeta (cc-pVQZ) basis set, appended with core correlation and relativistic corrections, we propose the relative energies: T(e)(FeNC), (6)Delta(0)<(6)Pi (2300 cm(-1))<(4)Delta (2700 cm(-1))<(4)Pi (4200 cm(-1))<(4)Sigma(-); and T(e)(FeCN), (6)Delta(0)<(6)Pi (1800 cm(-1))<(4)Delta (2500 cm(-1))<(4)Pi (2900 cm(-1))<(4)Sigma(-). The (4)Delta and (4)Pi states have massive multireference character, arising mostly from 11sigma-->12sigma promotions, whereas the sextet states are dominated by single electronic configurations. The single-reference CCSDT-3 (coupled cluster singles and doubles with iterative partial triples) method appears to significantly overshoot the stabilization of the quartet states provided by both static and dynamical correlation. The (4,6)Delta and (4,6)Pi states of both isomers are rather ionic, and all have dipole moments near 5 D. On the ground (6)Delta surface, FeNC is predicted to lie 0.6 kcal mol(-1) below FeCN, and the classical barrier for isocyanide/cyanide isomerization is about 6.5 kcal mol(-1). Our data support the recent spectroscopic characterization by Lei and Dagdigian [J. Chem. Phys. 114, 2137 (2000)] of linear (6)Delta FeNC as the first experimentally observed transition-metal monoisocyanide. Their assignments for the ground term symbol, isotopomeric rotational constants, and the Fe-N omega(3) stretching frequency are confirmed; however, we find rather different structural parameters for (6)Delta FeNC:r(e)(Fe-N)=1.940 A and r(N-C)=1.182 A at the cc-pVQZ MRCISD+Q level. Our results also reveal that the observed band of FeNC originating at 27 236 cm(-1) should have an analog in FeCN near 23 800 cm(-1) of almost equal intensity. Therefore, both thermodynamic stability and absorption intensity factors favor the eventual observation of FeCN via a (6)Pi<--(6)Delta transition in the near-UV.

Journal Article↗
PubMed2004

Binding energies of small lithium clusters (Li(n)) and hydrogenated lithium clusters (Li(n)H).

Large coupled cluster computations utilizing the Dunning weighted correlation-consistent polarized core-valence (cc-pwCVXZ) hierarchy of basis sets have been conducted, resulting in a panoply of internally consistent geometries and atomization energies for small Li(n) and Li(n)H (n=1-4) clusters. In contrast to previous ab initio results, we predict a monotonic increase in atomization energies per atom with increasing cluster size for lithium clusters, in accordance with the historical Knudsen-effusion measurements of Wu. For hydrogenated lithium clusters, our results support previous theoretical work concerning the relatively low atomization energy per atom for Li(2)H compared to LiH and Li(3)H. The CCSD(T)/cc-pwCVQZ atomization energies for LiH, Li(2)H, Li(3)H, and the most stable isomer of Li(4)H, including zero-point energy corrections, are 55.7, 79.6, 113.0, and 130.6 kcal/mol, respectively. The latter results are not consistent with the most recent experiments of Wu.

Journal Article↗
PubMed2004

The microwave and infrared spectroscopy of benzaldehyde: conflict between theory and experimental deductions.

Recently, it has been proposed that ab initio calculations cannot accurately treat molecules comprised of a benzene ring with a pi-conjugated substituent, for example, benzaldehyde. Theoretical predictions of the benzaldehyde barrier to internal rotation are typically a factor of 2 too high in comparison to the experimental values of 4.67 (infared) and 4.90 (microwave) kcal mol(-1). However, both experiments use Pitzer's 1946 model to compute the reduced moment of inertia and employ the experimentally observed torsional frequency to deduce benzaldehyde's rotational barrier. When Pitzer's model is applied to a system with a nonconjugated functional group, such as phenol, the model and theoretical values are in close agreement. Therefore, we conclude the model may not account for conjugation between the substituent and the pi-system of benzene. The experimental values of the benzaldehyde rotational barrier are therefore misleading. The true rotational barrier lies closer to the theoretically extrapolated limit of 7.7 kcal mol(-1), based on coupled cluster theory.

Journal Article↗
PubMed2004

The arsenic clusters Asn (n = 1-5) and their anions: structures, thermochemistry, and electron affinities.

The molecular structures, electron affinities, and dissociation energies of the As(n)/As(-) (n) (n = 1-5) species have been examined using six density functional theory (DFT) methods. The basis set used in this work is of double-zeta plus polarization quality with additional diffuse s- and p-type functions, denoted DZP++. These methods have been carefully calibrated (Chem Rev 2002, 102, 231) for the prediction of electron affinities. The geometries are fully optimized with each DFT method independently. Three different types of the neutral-anion energy separations reported in this work are the adiabatic electron affinity (EA(ad)), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). The first dissociation energies D(e)(As(n-1)-As) for the neutral As(n) species, as well as those D(e)(As(-) (n-1)-As) and D(e) (As(n-1)-As(-)) for the anionic As(-) (n) species, have also been reported. The most reliable adiabatic electron affinities, obtained at the DZP++ BLYP level of theory, are 0.90 (As), 0.74 (As(2)), 1.30 (As(3)), 0.49 (As(4)), and 3.03 eV (As(5)), respectively. These EA(ad) values for As, As(2), and As(4) are in good agreement with experiment (average absolute error 0.09 eV), but that for As(3) is a bit smaller than the experimental value (1.45 +/- 0.03 eV). The first dissociation energies for the neutral arsenic clusters predicted by the B3LYP method are 3.93 eV (As(2)), 2.04 eV (As(3)), 3.88 eV (As(4)), and 1.49 eV (As(5)). Compared with the available experimental dissociation energies for the neutral clusters, the theoretical predictions are excellent. Two dissociation limits are possible for the arsenic cluster anions. The atomic arsenic results are 3.91 eV (As(-) (2) --> As(-) + As), 2.46 eV (As(-) (3) --> As(-) (2) + As), 3.14 eV (As(-) (4) --> As(-) (3) + As), and 4.01 eV (As(-) (5) --> As(-) (4) + As). For dissociation to neutral arsenic clusters, the predicted dissociation energies are 2.43 eV (As(-) (3) --> As(2) + As(-)), 3.53 eV (As(-) (4) --> As(3) + As(-)), and 3.67 eV (As(-) (5) --> As(4) + As(-)). For the vibrational frequencies of the As(n) series, the BP86 and B3LYP methods produce good results compared with the limited experiments, so the other predictions with these methods should be reliable.

Journal Article↗
PubMed2003

Nonplanarity at tri-coordinated aluminum and gallium: cyclic structures for X3Hn(m) (X = B, Al, Ga).

Structures and energies of X3H3(2-), X3H4-, X3H5, and X3H6+ (X = B, Al and Ga) were investigated theoretically at B3LYP/6-311G(d) level. The global minimum structures of B are not found to be global minima for Al and Ga. The hydrides of the heavier elements Al and Ga have shown a total of seven, six and eight minima for X3H3(2-), X3H(4-), and X3H5, respectively. However, X3H(6+) has three and four minima for Al and Ga, respectively. The nonplanar arrangements of hydrogens with respect to X3 ring is found to be very common for Al and Ga species. Similarly, species with lone pairs on heavy atoms dominate the potential energy surfaces of Al and Ga three-ring systems. The first example of a structure with tri-coordinate pyramidal arrangement at Al and Ga is found in X3H(4-) (2g), contrary to the conventional wisdom of C3H3+, B3H3, etc. The influence of pi-delocalization in stabilizing the structures decreases from X3H3(2-) to X3H6+ for heavier elements Al and Ga. In general, minimum energy structures of X3H4-, X3H5, and X3H6+ may be arrived at by protonating the minimum energy structures sequentially starting from X3H3(2-). The resonance stabilization energy (RSE) for the global minimum structures (or nearest structures to global minimum which contains pi-delocalization) is computed using isodesmic equations.

Journal Article↗
PubMed2003

Mono- and dibridged isomers of Si2H3 and Si2H4: the true ground state global minima. Theory and experiment in concert.

Highly correlated ab initio coupled-cluster theories (e.g., CCSD(T), CCSDT) were applied on the ground electronic states of Si(2)H(3) and Si(2)H(4), with substantive basis sets. A total of 10 isomers, which include mono- and dibridged structures, were investigated. Scalar relativistic corrections and zero-point vibrational energy corrections were included to predict reliable energetics. For Si(2)H(3), we predict an unanticipated monobridged H(2)Si-H-Si-like structure (C(s), (2)A'') to be the lowest energy isomer, in constrast to previous studies which concluded that either H(3)Si-Si (C(s), (2)A'') or near-planar H(2)Si-SiH (C(1), (2)A) is the global minimum. Our results confirm that the disilene isomer, H(2)Si-SiH(2), is the lowest energy isomer for Si(2)H(4) and that it has a trans-bent structure (C(2)(h), (1)A(g)). In addition to the much studied silylsilylene, H(3)Si-SiH, we also find that a new monobridged isomer H(2)Si-H-SiH (C(1), (1)A, designated 2c) is a minimum on the potential energy surface and that it has comparable stability; both isomers are predicted to lie about 7 kcal/mol above disilene. By means of Fourier transform microwave spectroscopy of a supersonic molecular beam, the rotational spectrum of this novel Si(2)H(4) isomer has recently been measured in the laboratory, as has that of the planar H(2)Si-SiH radical. Harmonic vibrational frequencies as well as infrared intensities of all 10 isomers were determined at the cc-pVTZ CCSD(T) level.

Journal Article↗
PubMed2003

On the viability of small endohedral hydrocarbon cage complexes: X@C4H4, X@C8H8, X@C8H14, X@C10H16, X@C12H12, AND X@C16H16.

Small hydrocarbon complexes (X@cage) incorporating cage-centered endohedral atoms and ions (X = H(+), H, He, Ne, Ar, Li(0,+), Be(0,+,2+), Na(0,+), Mg(0,+,2+)) have been studied at the B3LYP/6-31G(d) hybrid HF/DFT level of theory. No tetrahedrane (C(4)H(4), T(d)()) endohedral complexes are minima, not even with the very small hydrogen atom or beryllium dication. Cubane (C(8)H(8), O(h)()) and bicyclo[2.2.2]octane (C(8)H(14), D(3)(h)()) minima are limited to encapsulating species smaller than Ne and Na(+). Despite its intermediate size, adamantane (C(10)H(16), T(d)()) can enclose a wide variety of endohedral atoms and ions including H, He, Ne, Li(0,+), Be(0,+,2+), Na(0,+), and Mg(2+). In contrast, the truncated tetrahedrane (C(12)H(12), T(d)()) encapsulates fewer species, while the D(4)(d)() symmetric C(16)H(16) hydrocarbon cage (see Table of Contents graphic) encapsulates all but the larger Be, Mg, and Mg(+) species. The host cages have more compact geometries when metal atoms, rather than cations, are inside. This is due to electron donation from the endohedral metals into C-C bonding and C-H antibonding cage molecular orbitals. The relative stabilities of endohedral minima are evaluated by comparing their energies (E(endo)) to the sum of their isolated components (E(inc) = E(endo) - E(cage) - E(x)) and to their exohedral isomer energies (E(isom) = E(endo) - E(exo)). Although exohedral binding is preferred to endohedral encapsulation without exception (i.e., E(isom) is always exothermic), Be(2+)@C(10)H(16) (T(d)(); -235.5 kcal/mol), Li(+)@C(12)H(12) (T(d)(); 50.2 kcal/mol), Be(2+)@C(12)H(12) (T(d)(); -181.2 kcal/mol), Mg(2+)@C(12)H(12) (T(d)(); -45.0 kcal/mol), Li(+)@C(16)H(16) (D(4)(d)(); 13.3 kcal/mol), Be(+)@C(16)H(16) (C(4)(v)(); 31.8 kcal/mol), Be(2+)@C(16)H(16) (D(4)(d)(); -239.2 kcal/mol), and Mg(2+)@C(16)H(16) (D(4)(d)(); -37.7 kcal/mol) are relatively stable as compared to experimentally known He@C(20)H(20) (I(h)()), which has an E(inc) = 37.9 kcal/mol and E(isom) = -35.4 kcal/mol. Overall, endohedral cage complexes with low parent cage strain energies, large cage internal cavity volumes, and a small, highly charged guest species are the most viable synthetic targets.

Journal Article↗