Breakthrough characteristics of volatile organic compounds in the -10 to + 170 degrees C temperature range on Tenax TA determined by microtrap technology

In this work the breakthrough volumes (BTVs) of volatile organic compounds (VOCs) on Tenax TA were determined in the -10 to 170 degrees C temperature range by using microtrap (MT) technology. The MT technology allowed experimental investigation of the temperature dependence of BTVs. Along with the BTV data, we also discuss the thermodynamics of the temperature dependence of the BTV through a two-parameter equation In (BTV) = A1/T + C1 where T is temperature (K), A1 = -deltaH/R where deltaH is enthalpy of sorption and C1 is constant. This equation fitted well the experimental results with R2 values between 0.9737 (acetone) and 0.9995 (dimethyl disulfide), with n between 6 and 11. However, for n-pentane, n-hexane and 1-hexene it proved that a three parameter equation In (BTV) = A2/ T + BTB + C2 fitted better to the experimental results, with A2 = -deltaHT0/R, B = deltaCp/R, TB = 1n(T/T0) + (T0 - T)/T, C2 a constant, deltaHT0 the adsorption enthalpy at reference temperature T0 and deltaCp the difference in the molar heat capacity of compound under investigation between the sorbed and the free gas phase state. The statistical analysis showed for example for n-pentane now R2 = 0.9969 instead of R2 = 0.9746, and Fisher statistics F = 487 instead of F = 153, with a significance level P = 0.018 for the third parameter. The results show that microtrap technology well serves as a technology to get information on temperature dependence of BTVs in an extended range. Simultaneously, it turns out that MT technology, extending the operational temperature range, is well served by a careful investigation of the temperature dependence models of BTVs.     

A dual-level ab initio and hybrid density functional theory dynamics study on the unimolecular decomposition reaction C2H5O-->CH2O + CH3

We present a direct ab initio and hybrid density functional theory dynamics study of the thermal rate constants of the unimolecular decomposition reaction of C2H5O-->CH2O + CH3 at a high-pressure limit. MPW1K/6-31+G(d,p), MP2/6-31+G(d,p), and MP2(full)/6-31G(d) methods were employed to optimize the geometries of all stationary points and to calculate the minimum energy path (MEP). The energies of all the stationary points were refined at a series of multicoefficient and multilevel methods. Among all methods, the QCISD(T)/aug-cc-pVTZ energies are in good agreement with the available experimental data. The rate constants were evaluated based on the energetics from the QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p) level of theory using both microcanonical variational transition state theory (microVT) and RRKM theory with the Eckart tunneling correction in the temperature range of 300-2500 K. The calculated rate constants at the QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p) level of theory are in good consistent with experimental data. The fitted three-parameter Arrhenius expression from the microVT/Eckart rate constants in the temperature range 200-2500 K is k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1). The falloff curves of pressure-dependent rate constants are performed using master-equation method within the temperature range of 391-471 K. The calculated results are in good agreement with the available experimental data.     

Experimental and theoretical studies of rate coefficients for the reaction O(3P)+C2H5OH at high temperatures

Rate coefficients of the reaction O(3P)+C2H5OH in the temperature range 782-1410 K were determined using a diaphragmless shock tube. O atoms were generated by photolysis of SO2 at 193 nm with an ArF excimer laser; their concentrations were monitored via atomic resonance absorption. Our data in the range 886-1410 K are new. Combined with previous measurements at low temperature, rate coefficients determined for the temperature range 297-1410 K are represented by the following equation: k(T)=(2.89+/-0.09)x10(-16)T1.62 exp[-(1210+/-90)/T] cm3 molecule(-1) s(-1); listed errors represent one standard deviation in fitting. Theoretical calculations at the CCSD(T)/6-311+G(3df, 2p)//B3LYP/6-311+G(3df) level predict potential energies of various reaction paths. Rate coefficients are predicted with the canonical variational transition state (CVT) theory with the small curvature tunneling correction (SCT) method. Reaction paths associated with trans and gauche conformations are both identified. Predicted total rate coefficients, 1.60 x 10(-22)T3.50 exp(16/T) cm3 molecule(-1) s(-1) for the range 300-3000 K, agree satisfactorily with experimental observations. The branching ratios of three accessible reaction channels forming CH3CHOH+OH (1a), CH2CH2OH+OH (1b), and CH3CH2O+OH (1c) are predicted to vary distinctively with temperature. Below 500 K, reaction 1a is the predominant path; the branching ratios of reactions 1b,c become approximately 40% and approximately 11%, respectively, at 2000 K.     

Effect of boron on the surface tension of molten silicon and its temperature coefficient

The influence of boron concentration (C(B)/mass%) on the surface tension of molten silicon has been investigated with the sessile drop method under oxygen partial pressure P(O(2))=1.62x10(-25)-2.63x10(-22) MPa, and the results can be summarized as follows. The surface tension increases with C(B) in the range below 2.09 mass%, and the maximum increase rate of the surface tension is about 30 mN m(-1)(mass% C(B))(-1). The temperature coefficient of the surface tension, ( partial differential sigma/ partial differential T)C(B), was found to increase with the boron concentration in molten silicon. At the interface between molten silicon and the BN substrate, a discontinuous Si(3)N(4) layer was reckoned to form and the layer might prevent BN from dissolving into the molten silicon. Since dissolved boron from the BN substrate into the molten silicon is below 0.054 mass% and the associated increase in surface tension is below 1.5 mN m(-1), the contamination from the BN substrate on the surface tension can be ignored. The relation between the surface tension and C(B) indicates negative adsorption of boron and can be well described by combining the Gibbs adsorption isotherm with the Langmuir isotherm.     

Experimental and theoretical study of the carbon-13 and deuterium kinetic isotope effects in the Cl and OH reactions of CH3F

A laser flash photolysis-resonance fluorescence technique has been employed to determine absolute rate coefficients for the CH3F + Cl reaction in N2 bath gas in the temperature range of 200-700 K and pressure range of 33-133 hPa. The data were fitted to a modified Arrhenius expression k(T) = 1.14 x 10(-12) x (T/298)2.26 exp{-313/T}. The OH and Cl reaction rates of (13)CH3F and CD3F have been measured by long-path FTIR spectroscopy relative to CH3F at 298 +/- 2 K and 1013 +/- 10 hPa in purified air. The FTIR spectra were fitted using a nonlinear least-squares spectral fitting method including line data from the HITRAN database and measured infrared spectra as references. The relative reaction rates defined by alpha = k(light)/k(heavy) were determined to be k(OH+CH3F)/k(OH+CD3F) = 4.067 +/- 0.018, k(OH+CH3F)/k(OH+(13)CH3F) = 1.067 +/- 0.006, k(Cl+CH3F)/k(Cl+CD3F) = 5.11 +/- 0.07, and k(Cl+CH3F)/k(Cl+(13)CH3F) = 1.016 +/- 0.006. The carbon-13 and deuterium kinetic isotope effects in the OH and Cl reactions of CH3F have been further investigated by quantum chemistry methods and variational transition state theory.     

Shock tube study on the thermal decomposition of ethanol

The thermal decomposition of C(2)H(5)OH highly diluted in Ar (1 and 3 ppm) has been studied by monitoring H atoms using the atomic resonance absorption spectrometry (ARAS) technique behind reflected shock waves over the temperature range 1450-1760 K at fixed pressure: 1, 1.45, and 2 atm. The rate constant and the product branching fractions have been determined by analyzing temporal profiles of H atoms; the effect of the secondary reactions on the results has been examined by using a detailed reaction mechanism composed of 103 elementary reactions. The apparent rate constant of ethanol decomposition can be expressed as k(1)/s(-1) = (5.28 ± 0.14) × 10(10) exp[-(23,530 ± 980)/T] (T = 1450-1670 K, P = 1-2 atm) without a detectable pressure dependence within the tested pressure range of this study. Branching fractions for producing CH(3) + CH(2)OH (1a) and H(2)O + C(2)H(4) (1b) have been examined by a quantitative measurement of H atoms produced in the successive decompositions of the products CH(2)OH (1a): the pressure dependence of the branching fraction for channel 1a is obtained by a linear least-squares analysis of the experimental data and can be expressed as φ(1a) = (0.71 ± 0.07) - (826 ± 116)/T, (0.92 ± 0.04) - (1108 ± 70)/T, and (1.02 ± 0.10) - (1229 ± 168)/T for T = 1450-1760 K, at P = 0.99, 1.45, and 2.0 atm, respectively. The rate constant obtained in this study is found to be consistent with previous theoretical and experimental results; however, the pressure dependence of the branching fraction obtained in this study is smaller than those of previous theoretical works. Modification of the parameters for the decomposition rate in the falloff region is suggested to be important to improve the practical modeling of the pyrolysis and combustion of ethanol.     

Experimental and theoretical study of the reaction of the ethynyl radical with nitrous oxide, C2H + N2O

We investigated the rate constants and reaction mechanism of the gas phase reaction between the ethynyl radical and nitrous oxide (C(2)H + N(2)O) using both experimental methods and electronic structure calculations. A pulsed-laser photolysis/chemiluminescence technique was used to determine the absolute rate coefficient over the temperature range 570 K to 836 K. In this experimental temperature range, the measured temperature dependence of the overall rate constants can be expressed as: k(T) (C(2)H + N(2)O) = 2.93 × 10(-11) exp((-4000 ± 1100) K/T) cm(3) s(-1) (95% statistical confidence). Portions of the C(2)H + N(2)O potential energy surface (PES), containing low-energy pathways, were constructed using the composite G3B3 method. A multi-step reaction route leading to the products HCCO + N(2) is clearly preferred. The high selectivity between product channels favouring N(2) formation occurs very early. The pathway corresponds to the addition of the terminal C atom of C(2)H to the terminal N atom of N(2)O. Refined calculations using the coupled-cluster theory whose electronic energies were extrapolated to the complete basis set limit CCSD(T)/CBS led to an energy barrier of 6.0 kcal mol(-1) for the entrance channel. The overall rate constant was also determined by application of transition-state theory and Rice-Ramsperger-Kassel-Marcus (RRKM) statistical analyses to the PES. The computed rate constants have similar temperature dependence to the experimental values, though were somewhat lower.     

Kinetics of the O + ICN reaction

The kinetics of the O + ICN reaction was studied using a relative rate method, with O + C(2)H(2) as the competing reaction. Carbon monoxide products formed in the competing reaction and subsequent secondary chemistry were detected as a function of reagent ICN pressure to obtain total rate constants for the O + ICN reaction. Analysis of the experimental data yields rate constants of k(1) = (3.7 ± 1.0 to 26.2 ± 4.0) × 10(-14) cm(3) molecule(-1) s(-1) over the total pressure range 1.5-9.5 Torr. Product channel NCO + I, the only bimolecular exothermic channel of the reaction, was investigated by detection of N(2)O in the presence of NO and found to be insignificant. An ab initio calculation of the potential energy surface (PES) of the reaction at the CCSD(T)/CEP-31G//DFT-B3LYP/CEP-31G level of theory was also performed. The pathways leading to bimolecular product channels are kinetically unfavorable. Formation and subsequent stabilization of an ICNO adduct species appears to dominate the reaction, in agreement with the experimentally observed pressure dependent rate constants.     

Ab initio study on mechanisms and kinetics for reaction of NCS with NO

The mechanisms and kinetics of the reaction of a thiocyanato radical (NCS) with NO were investigated by a high-level ab initio molecular orbital method in conjunction with variational RRKM calculations. The species involved were optimized at the B3LYP/6-311++G(3df,2p) level, and their single-point energies were refined by the CCSD(T)/aug-cc-PVQZ//B3LYP/6-311+G(3df,2p) method. Our calculated results indicate favorable pathways for the formation of several isomers of an NCSNO complex. Formation of OCS + N 2 also is possible, although this pathway involves a substantial energy barrier. The predicted total rate constants, k total, at a 2 torr He pressure can be represented by the following equations: k total = 9.74 x 10 (26) T (-13.88) exp(-6.53 (kcal mol (-1))/ RT) at T = 298-950 K and 1.17 x 10 (-22) T (2.52) exp(-6.86 (kcal mol (-1))/ RT) at T = 960-3000 K, in units of cm (3) molecule (-1) s (-1), and the predicted values are in good agreement with the experimental results in the temperature range of 298-468 K. The calculated results clearly indicate that the branching ratio for R M1 in the temperature range of 298-950 K has the largest value ( R M1 accounts for 0.53-0.39). However, in the higher temperature range (960-3000 K), the formation of OCS + N 2 ( P5) with branching ratio R P5 (0.40-0.79) becomes dominant. The rate constants for key individual product channels are provided for different temperature and pressure conditions.     

A molecular dynamics study of the thermodynamic properties of calcium apatites. 1. Hexagonal phases

Structural and thermodynamic properties of crystal hexagonal calcium apatites, Ca10(PO4)6(X)2 (X = OH, F, Cl, Br), were investigated using an all-atom Born-Huggins-Mayer potential by a molecular dynamics technique. The accuracy of the model at room temperature and atmospheric pressure was checked against crystal structural data, with maximum deviations of ca. 4% for the haloapatites and 8% for hydroxyapatite. The standard molar lattice enthalpy, delta(lat)H298(o), of the apatites was calculated and compared with previously published experimental results, the agreement being better than 2%. The molar heat capacity at constant pressure, C(p,m), in the range 298-1298 K, was estimated from the plot of the molar enthalpy of the crystal as a function of temperature, H(m) = (H(m,298) - 298C(p,m)) + C(p,m)T, yielding C(p,m) = 694 +/- 68 J x mol(-1) x K(-1), C(p,m) = 646 +/- 26 J x mol(-1) x K(-1), C(p,m) = 530 +/- 34 J x mol(-1) x K(-1), and C(p,m) = 811 +/- 42 J x mol(-1) x K(-1) for hydroxy-, fluor-, chlor-, and bromapatite, respectively. High-pressure simulation runs, in the range 0.5-75 kbar, were performed in order to estimate the isothermal compressibility coefficient, kappaT, of those compounds. The deformation of the compressed solids is always elastically anisotropic, with BrAp exhibiting a markedly different behavior from those displayed by HOAp and ClAp. High-pressure p-V data were fitted to the Parsafar-Mason equation of state with an accuracy better than 1%.     

Theoretical and experimental investigation of the thermochemistry of CrO2(OH)2(g)

In this paper, we report the results of equilibrium pressure measurements designed to identify the volatile species in the Cr-O-H system and to resolve some of the discrepancies in existing experimental data. In addition, ab initio calculations were performed to lend confidence to a theoretical approach for predicting the thermochemistry of chromium-containing compounds. Equilibrium pressure data for CrO2(OH)2 were measured by the transpiration technique for the reaction 0.5Cr2O3(s) + 0.75O2(g) + H2O(g) = CrO2(OH)2(g) over a temperature range of 573 to 1173 K at 1 bar total pressure. Inductively coupled plasma atomic emission spectroscopy (ICP-AES) was used to analyze the condensate in order to quantify the concentration of Cr-containing volatile species. The resulting experimentally measured thermodynamic functions are compared to those computed using B3LYP density functional theory and the coupled-cluster singles and doubles method with a perturbative correction for connected triple substitutions [CCSD(T)].     

Evaluation of thermal analysis equipment for the determination of vapor pressure and heat of vaporization

A rapid, relatively simple method for determining vapor pressure and heat of vaporization on small amounts of organic compounds is described. A DuPont 900 differential thermal analyzer (DTA), a Perkin—Elmer Model DSC-1B differential scanning calorimeter (DSC), and a Thomas—Hoover (TH) melting point apparatus were evaluated in this work. Vapor pressure data for a wide variety of organic liquids were obtained by measuring the boiling points of the liquids at pressures ranging from 20 to 735 torr. A computer was used to rapidly plot the experimental data. The average deviations of boiling points from the literature values were 2.3°C for the DTA 1.2°C for the DSC, and 1.5°C for the TH. The vapor pressure data were used to solve the Haggenmacher equation for heat of vaporization (ΔH_v). The deviations of the experimental values for ΔH_v. from the literature values were 5.5%, 8.3%. and 3.3% for the DTA, DSC, and TH methods, respectively.     

Quantum chemical and theoretical kinetics study of the O(3P) + C2H2 reaction: a multistate process

The potential energy surfaces of the two lowest-lying triplet electronic surfaces 3A' and 3A' for the O(3P) + C2H2 reaction were theoretically reinvestigated, using various quantum chemical methods including CCSD(T), QCISD, CBS-QCI/APNO, CBS-QB3, G2M(CC,MP2), DFT-B3LYP and CASSCF. An efficient reaction pathway on the electronically excited 3A' surface resulting in H(2S) + HCCO(A2A') was newly identified and is predicted to play an important role at higher temperatures. The primary product distribution for the multistate multiwell reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. Allowing for nonstatistical behavior of the internal rotation mode of the initial 3A' adducts, our computed primary-product distributions agree well with the available experimental results, i.e., ca. 80% H(2S) + HCCO(X2A' + A2A') and 20% CH2(X3B1) + CO(X1sigma+) independent of temperature and pressure over the wide 300-2000 K and 0-10 atm ranges. The thermal rate coefficient k(O + C2H2) at 200-2000 K was computed using multistate transition state theory: k(T) = 6.14 x 10(-15)T (1.28) exp(-1244 K/T) cm3 molecule(-1) s(-1); this expression, obtained after reducing the CBS-QCI/APNO ab initio entrance barriers by 0.5 kcal/mol, quasi-perfectly matches the experimental k(T) data over the entire 200-2000 K range, spanning 3 orders of magnitude.     

A theoretical study of the effect of a tetraalkylammonium counterion on the hydrogen bond strength in Z-hydrogen maleate.

High-level ab initio calculations (B3LYP/6-31+G and QCISD(T)/6-311+G**) were carried out to resolve the disagreement between recent experimental and computational estimates of the relative strength of the intramolecular hydrogen bond in Z-hydrogen maleate anion with respect to the normal hydrogen bond in maleic acid. The computational estimates for the strength of the intramolecular hydrogen bond in the gas-phase maleate anion are in a range of 14-28 kcal/mol depending on the choice of the reference structure. Computational data suggest that the electrostatic influence of a counterion such as a tetraalkylammonium cation can considerably weaken the hydrogen bonding interaction (by 1.5-2 times) in the complexed hydrogen maleate anion relative to that in the naked anion. The estimated internal H-bonding energies for a series of Z-maleate/R4N+ salts (R = CH3, C2H5, CH3CH2CH2CH2) range from 8 to 13 kcal/mol. The calculated energy differences between the E- and Z-hydrogen maleates complexed to Me4N+, Et4N+, and Bu4N+ cation are 4.9 (B3LYP/6-31+G(d,p)) and 5.7 and 5.8 kcal/mol (B3LYP/6-31G(d)). It is also demonstrated that the sodium cation exerts a similar electrostatic influence on the hydrogen bond strength in bifluoride anion (FHF-). The present study shows that while low-barrier short hydrogen bonds can exist in the gas phase (the barrier for the hydrogen transfer in maleate anion is only 0.2 kcal/mol at the QCISD(T)/6-311+G//QCISD/6-31+G level), whether they can also be strong in condensed media or not depends on how their interactions with their immediate environment affect their strength.     

Ionization of five aqueous fluorobenzoic acids: Conductance and thermodynamics

The three monofluorobenzoic acids together with 2,4-difluoro and 2,6-difluorobenzoic acids in aqueous solution are the subject of precision conductance measurements. The experimental data are analyzed to give ionization constants and limiting conductances at temperatures from 0 to 100°C. Walden products for the acid anions are derived from the limiting conductances while the ionization consatants are fitted by statistical methods to the function pK _a (m)=A+B/T+ C logT+DT. Only the 2,6- acid requires the fourth term of the function to fit the data to a precision of better than 0.03%. Mathematical analysis of the pK function gives the standard changes in enthalpy, entropy, and heat capacity. All the acids studied are more acidic than the parent, benzoic acid, as well as more acidic than the isoelectronic methylbenzoic acids. In general the increased acidity is tied to decreases in enthalpy while entropy changes on ionization differn little from those found for the parent acid.     

Parametric effects of the potential energy function on the geometrical features of ternary lennard-jones clusters

The impact of parameters in potential function for describing atomic or molecular clusters is complex due to the complicated potential energy surface. Ternary Lennard-Jones (TLJ) A(l)B(m)C(n) clusters with two-body potential are investigated to study the effect of parameters. In the potential, the size parameter (σ(AA)) of A atoms is fixed, and corresponding parameters of B and C atoms (relative to A atoms), i.e., σ(BB)/σ(AA) and σ(CC)/σ(AA) > 1.00, are used to control the atomic interaction among A, B, and C atoms in TLJ clusters. The minimum energy configurations of A(l)B(m)C(n) clusters with different species are optimized by adaptive immune optimization algorithm. Ternary cluster structures, bonds, and energies of the putative minima are studied. The results show that two different structures based on double-icosahedra are found in 30-atom TLJ clusters. Furthermore, with increasing potential size parameters of B and C atoms, A atoms tend to be more compact for the increasing numbers of A-A bonds, but the short-range attractive part in TLJ clusters becomes weaker. To lower the potential energy, B and C atoms grow around the A atoms in pursuit of a compact configuration. The results are also approved in A(l)B(m)C(n) (l + m + n = 9-55) clusters and A(l)B(m)C(n) (l = 13, m + n = 42) clusters.     

用玻璃电极测定直到15 m的浓盐酸体系的活度系数

本文在25 ℃下, 对0.01844—15.02 mol kg~(-1)浓度范围内的盐酸体系, 用自制Ag|AgGl电极和氢电极, 测定了Harned电池的电动势。在此基础上, 以玻璃电极代替氢电极, 测量上述盐酸体系的电动势, 计算了HCl的平均活度系数和H_2O的活度。实验结果表明, 玻璃电极在浓盐酸体系的“标准电势”是稳定的, 且遵从Nernst响应, 重现性良好, 响应速度快; 在高浓度和稀溶液中交替测量时, 电势无明显偏差。这说明电极有较好的可复性。求得的HCl的活度系数与文献值偏差小于1.2％。玻璃电极的“标准电势”的标准偏差为0.016％。     

Effects of antimony on the surface tension of molten silicon

The effect of antimony concentration (C(Sb)/mass%) on the surface tension of molten silicon has been determined with the sessile drop method in the temperature range from 1693 to 1773 K and in the range of the oxygen partial pressure, Po(2), in an Ar atmosphere from 10(-23) to 10(-21) MPa. The results show that the surface tension of molten silicon decreases with increasing Sb concentration in the range of C(Sb)<0.9 mass%, which indicates positive adsorption of Sb in molten silicon and can be fairly described with the Szyszkowski's equation. The maximum decrease rate of surface tension is about 65 mN m(-1) (mass% C(Sb))(-1), and the temperature coefficient of surface tension, (partial differential sigma/ partial differential T)C(Sb), increases with increasing C(Sb). The evaporation of the systems was only observed between the melting points of antimony (904 K) and silicon (1683 K), and the surface tension presents no dependence on measuring time above the melting point of silicon.     

Metal electrodes for continuous amperometric measurement of free hydrofluoric acid in acidic solutions containing complexing ions

Titanium-based electrodes have been demonstrated to be useful for measuring free HF concentration in the range 0.01-0.07M at temperatures between 19 and 61 degrees C in acidic, complexed HF solutions, by measurement of the current density resulting when the electrode potential is approximately +0.6 V with respect to an inert electrode. The current density, i, is linearly related to the HF concentration when the solution is adequately stirred to remove diffusion effects at the electrode, and is independent of HNO(3) concentration up to 1.75M. The dependence on the absolute temperature is ln i = - A/T + B + C ln T, and the activation energy at an absolute temperature T(r) is R(A + T(r)C), where A,B and C are constants and R is the gas constant. For a Ti-2% Pd electrode, the activation energy for anodic dissolution is found to be 43.7 +/- 0.2 kJ/mole at 312 K. The equations relating the current density to HF concentration at different temperatures can be combined to give a single equation for [HF] as a function of both i and T.     

NH0-1+2-3离解能等的高级ab initio计算与评价

A recent experimental determination[1] of the dissociation energies (D0) for H2N-H, H2N+-H and H2N-H+, the ionization energies for NH3 and NH2 resulted in large deviations when compared with those of the earlier values and the QCISD(T)/6-311+G(3df,2p) ab initio calculations. We have performed some higher level ab initio calculations on these data by utilizing the Gaussian 92/DFT and Gaussian 94 pakages of programs and have assessed the available experimental values. Our calculations were carried out at the QCISD (TQ)/aug-cc-pVDZ, G2(QCI), QCISD(T)/6-311 ++G(3df,3pd) and QCISD(T)/aug-cc-pVTZ levels of theory. Geometries were optimized at both of the MP2(full)/6-31G(d) and the MP2(full)/6-31(d,p) levels, and were compared with those of the experiments if available. The MP2(full)/6-31G(d,p) tight-optimized geometries for the neutrals are closer to those of the experiments than those of the MP2 (full)/6-31G(d), and are in excellent agreement with the experimental results as shown in Table 1. In this case, we assumed that the optimized geometries for the cations would be better if p polarization functions are added to the hydrogen atoms. We firstly noted that the symmetry of the NH3+ cation was D3h, other than Cs. as reported in ref.[1]. All of the zero-point energies and the final geometries are calculated at the MP2(full)/6-31G(d,p) level of theory. We have also repeated the QCISD(T )/6-311 + G(3df,2p) calculations of ref. [1], because we could not identify their level of goemetry optimization. It is found that the total energy, -55.244 19 Hartrees, for NH2+(1A1 ) in ref.[1] might be in error. Our result is -55.336 29 Hartrees at the same level of theory. At our highest level [QCISD(T)/aug-cc-pVTZ] of calculations as shown in Table 3, the D0 (temperature at zero Kelvin) values of H2N-H, H2N+-H(3B1for NH2+ ) and H2N- H+ are 4.51, 5.49 and 8.00 eV, respectively. These data reported in re f.[1] were 4.97, 5.59 and 8.41 eV, respectively. Our result on D0(H2N-H) supports the work of ref.[2,3,5,6]. The ionization energies (IE) for NH3 and NH2 (3B1 for NH2+) at our highest level are 10.11 and 11.09 eV while in ref.[1] were 10.16 and 10.78 eV, respectively. For the latter, our result supports the experiment of ref.[3]. Our predicted D0 for HN2+-H and IE for NH2 (1A1 for each NH2+) are 6.80 and 12.39 eV, respectively. These values differ greatly from the predicted values (9.29 and 14.88 eV) of ref.[1] where the total energy of NH2+(1A1) might be in error. The D0 value for HN-H has not been found in ref.[1]. Our result supports the work of ref.[3]. We have also derived all of these values at the temperature of 298K and under the pressure of 101kPa at several levels of thoery as shown in Table 3. On examining the experiment of ref.[1] in detail, it is easy to find that all of the larger deviations might be from a too high value of the appearance potential of proton AP(H+). Indeed, ref.[1] has mentioned that the determintion of AP(H+), due to kinetic shift, would lead to a hihger value for the dissociation energy as has been pointed out by Berkowitz and Ruscic. In this work, we concluded that, besides some mistakes in the theoretical calculations of ref.[1], the dissociation energies for H2N-H and H2N-H+,the IE for NH2 (3B1 for NH2+) might also be unreliable and need to be re-examined.