Structure and dynamics of Au+ ion in aqueous solution: ab initio QM/MM MD simulations

Structure and dynamics of hydrated Au(+) have been investigated by means of molecular dynamics simulations based on ab initio quantum mechanical molecular mechanical forces at Hartree-Fock level for the treatment of the first hydration shell. The outer region of the system was described using a newly constructed classical three-body corrected potential. The structure was evaluated in terms of radial and angular distribution functions and coordination number distributions. Water exchange processes between coordination shells and bulk indicate a very labile structure of the first hydration shell whose average coordination number of 4.7 is a mixture of 3-, 4-, 5-, 6-, and 7-coordinated species. Fast water exchange reactions between first and second hydration shell occur, and the second hydration shell is exceptionally large. Therefore, the mean residence time of water molecules in the first hydration shell (5.6 ps/7.5 ps for t*= 0.5 ps/2.0 ps) is shorter than that in the second shell (9.4 ps/21.2 ps for t*= 0.5 ps/2.0 ps), leading to a quite specific picture of a "structure-breaking" effect.     

MP2, DFT and ab initio calculations on thioxanthone

Density functional theory (DFT), HF and MP2 calculations have been carried out to investigate thioxanthone molecule using the standard 6-31+G(d,p) basis set. The results of MP2 calculations show a butterfly structure for thioxanthone. The calculated results show that the predicted geometry can well reproduce the structural parameters. The predicted vibrational frequencies were assigned and compared with experimental IR spectra. A good harmony between theory and experiment is found. The theoretical electronic absorption spectra have been calculated using CIS method. ¹³C and ¹H NMR of the title compound have been calculated by means of B3LYP density functional method with 6-31+G(d,p) basis set. The comparison of the experimental and the theoretical results indicate that density functional B3LYP method is able to provide satisfactory results for predicting NMR properties.     

Treatment of dilute clusters of methanol and water by ab initio quantum mechanical calculations

Large molecular clusters can be considered as intermediate states between gas and condensed phases, and information about them can help us understand condensed phases. In this paper, ab initio quantum mechanical methods have been used to examine clusters formed of methanol and water molecules. The main goal was to obtain information about the intermolecular interactions and the structure of methanol/water clusters at the molecular level. The large clusters (CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10)) containing one molecule of one component (methanol or water) and many (12, 10) molecules of the other component were considered. M?ller-Plesset perturbation theory (MP2) was used in the calculations. Several representative cluster geometries were optimized, and nearest-neighbor interaction energies were calculated for the geometries obtained in the first step. The results of the calculations were compared to the available experimental information regarding the liquid methanol/water mixtures and to the molecular dynamics and Monte Carlo simulations, and good agreement was found. For the CH(4)O...(H(2)O)(12) cluster, it was shown that the molecules of water can be subdivided into two classes: (i) H bonded to the central methanol molecule and (ii) not H bonded to the central methanol molecule. As expected, these two classes exhibited striking energy differences. Although they are located almost the same distance from the carbon atom of the central methanol molecule, they possess very different intermolecular interaction energies with the central molecule. The H bonding constitutes a dominant factor in the hydration of methanol in dilute aqueous solutions. For the H(2)O...(CH(4)O)(10) cluster, it was shown that the central molecule of water has almost three H bonds with the methanol molecules; this result differs from those in the literature that concluded that the average number of H bonds between a central water molecule and methanol molecules in dilute solutions of water in methanol is about two, with the water molecules being incorporated into the chains of methanol. In contrast, the present predictions revealed that the central water molecule is not incorporated into a chain of methanol molecules, but it can be the center of several (2-3) chains of methanol molecules. The molecules of methanol, which are not H bonded to the central water molecule, have characteristics similar to those of the methane molecules around a central water molecule in the H(2)O...(CH(4))(10) cluster. The ab initio quantum mechanical methods employed in this paper have provided detailed information about the H bonds in the clusters investigated. In particular, they provided full information about two types of H bonds between water and methanol molecules (in which the water or the methanol molecule is the proton donor), including information about their energies and lengths. The average numbers of the two types of H bonds in the CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10) clusters have been calculated. Such information could hardly be obtained with the simulation methods.     

The mechanism of formamide hydrolysis in water from ab initio calculations and simulations

The neutral hydrolysis of formamide in water is a suitable reference to quantify the efficiency of proteolytic enzymes. However, experimental data for this reaction has only very recently been obtained and the kinetic constant determined experimentally is significantly higher than that predicted by previous theoretical estimations. In this work, we have investigated in detail the possible mechanisms of this reaction. Several solvent models have been considered that represent a considerable improvement on those used in previous studies. Density functional and ab initio calculations have been carried out on a system which explicitly includes the first solvation shell of the formamide molecule. Its interaction with the bulk has been treated with the aid of a dielectric continuum model. Molecular dynamics simulations at the combined density functional/molecular mechanics level have been carried out in parallel to better understand the structure of the reaction intermediates in aqueous solution. Overall, the most favored mechanism predicted by our study involves two reaction steps. In the first step, the carbonyl group of the formamide molecule is hydrated to form a diol intermediate. The corresponding transition structure involves two water molecules. From this intermediate, a water-assisted proton transfer occurs from one of the hydroxy groups to the amino group. This reaction step may lead either to the formation of a new reaction intermediate with a marked zwitterionic character or to dissociation of the system into ammonia and formic acid. The zwitterionic intermediate dissociates quite easily but its lifetime is not negligible and it could play a role in the hydrolysis of substituted amides or peptides. The predicted pseudo-first-order kinetic constant for the rate-limiting step (the first step) of the hydrolysis reaction at 25 degrees C (3.9x10(-10) s(-1)) is in excellent agreement with experimental data (1.1x10(-10) s(-1)).     

Interaction of neutral and zwitterionic glycine with Zn2+ in gas phase: ab initio and SIBFA molecular mechanics calculations

The interaction of Zn²⁺ with glycine (Gly) in the gas phase is studied by a combination of ab initio and molecular mechanics techniques. The structures and energetics of the various isomers of the Gly–Zn²⁺ complex are first established via high‐level ab initio calculations. Two low‐energy isomers are characterized: one in which the metal ion interacts with the carboxylate end of zwitterionic glycine, and another in which it chelates the amino nitrogen and the carbonyl oygen of neutral glycine. These calculations lead to the first accurate value of the gas‐phase affinity of glycine for Zn²⁺. Ab initio calculations were also used to evaluate the performance of various implementations of the SIBFA force field. To assess the extent of transferability of the distributed multipoles and polarizabilities used in the SIBFA computations, two approaches are followed. In the first, approach (a), these quantities are extracted from the ab initio Hartree–Fock wave functions of glycine or its zwitterion in its entirety, and for each individual Zn²⁺‐binding conformation. In the second, approach (b), they are assembled from the appropriate constitutive fragments, namely methylamine and formic acid for neutral glycine, and protonated methylamine and formate for the zwitterion; they undergo the appropriate vector or matrix rotation to be assembled in the conformation studied. The values of the Zn²⁺–glycine interaction energies are compared to those resulting from ab initio SCF and MP2 computations using both the all‐electron 6‐311+G(2d,2p) basis set and an effective core potential together with the valence CEP 4‐31G(2d) basis set. Approach (a) values closely reproduce the ab initio ones, both in terms of the total interaction energies and of the individual components. Approach (b) can provide a similar match to ab initio interaction energies as does approach (a), provided that the two constitutive Gly building blocks are considered as separate entities having mutual interactions that are computed simultaneously with those occurring with Zn²⁺. Thus, the supermolecule is treated as a three‐body rather than a two‐body system. These results indicate that the current implementation of the SIBFA force field should be adequate to undertake accurate studies on zinc metallopeptides. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 963–973, 2000     

Dissociation of carbonic acid: gas phase energetics and mechanism from ab initio metadynamics simulations

A comprehensive metadynamics study of the energetics, stability, conformational changes, and mechanism of dissociation of gas phase carbonic acid, H2CO3, yields significant new insight into these reactions. The equilibrium geometries, vibrational frequencies, and conformer energies calculated using the density functional theory are in good agreement with the previous theoretical predictions. At 315 K, the cis-cis conformer has a very short life time and transforms easily to the cis-trans conformer through a change in the O=C-O-H dihedral angle. The energy difference between the trans-trans and cis-trans conformers is very small (approximately 1 kcal/mol), but the trans-trans conformer is resistant to dissociation to carbon dioxide and water. The cis-trans conformer has a relatively short path for one of its hydroxyl groups to accept the proton from the other end of the molecule, resulting in a lower activation barrier for dissociation. Comparison of the free and potential energies of dissociation shows that the entropic contribution to the dissociation energy is less than 10%. The potential energy barrier for dissociation of H2CO3 to CO2 and H2O from the metadynamics calculations is 5-6 kcal/mol lower than in previous 0 K studies, possibly due to a combination of a finite temperature and more efficient sampling of the energy landscape in the metadynamics calculations. Gas phase carbonic acid dissociation is triggered by the dehydroxylation of one of the hydroxyl groups, which reorients as it approaches the proton on the other end of the molecule, thus facilitating a favorable H-O-H angle for the formation of a product H2O molecule. The major atomic reorganization of the other part of the molecule is a gradual straightening of the O=C=O bond. The metadynamics results provide a basis for future simulation of the more challenging carbonic acid-water system.     

Rotationally invariant ab initio evaluation of Coulomb and exchange parameters for DFT+U calculations

Conventional density functional theory (DFT) fails for strongly correlated electron systems due to large intra-atomic self-interaction errors. The DFT+U method provides a means of overcoming these errors through the use of a parametrized potential that employs an exact treatment of quantum mechanical exchange interactions. The parameters that enter into this potential correspond to the spherically averaged intra-atomic Coulomb (U) and exchange (J) interactions. Recently, we developed an ab initio approach for evaluating these parameters on the basis of unrestricted Hartree-Fock (UHF) theory, which has the advantage of being free of self-interaction errors and does not require experimental input [Mosey and Carter, Phys. Rev. B 76, 155123 (2007)]. In this work, we build on that method to develop a more robust and convenient ab initio approach for evaluating U and J. The new technique employs a relationship between U and J and the Coulomb and exchange integrals evaluated using the entire set of UHF molecular orbitals (MOs) for the system. Employing the entire set of UHF MOs renders the method rotationally invariant and eliminates the difficulty in selecting unambiguously the MOs that correspond to localized states. These aspects overcome two significant deficiencies of our earlier method. The new technique is used to evaluate U and J for Cr(2)O(3), FeO, and Fe(2)O(3). The resulting values of U-J are close to empirical estimates of this quantity for each of these materials and are also similar to results of constrained DFT calculations. DFT+U calculations using the ab initio parameters yield results that are in good agreement with experiment. As such, this method offers a means of performing accurate and fully predictive DFT+U calculations of strongly correlated electron materials.     

Using MD snapshots in ab initio and DFT calculations: OH vibrations in the first hydration shell around Li+(aq)

The average OH stretching vibrational frequency for the water molecules in the first hydration shell around a Li(+) ion in a dilute aqueous solution was calculated by a hybrid molecular dynamics + quantum-mechanical ("MD + QM") approach. Using geometry configurations from a series of snapshots from an MD simulation, the anharmonic, uncoupled OH stretching frequencies were calculated for 100 first-shell OH oscillators at the B3LYP and HF/6-31G(d,p) levels of theory, explicitly including the first shell and the relevant second shell water molecules into charge-embedded supermolecular QM calculations. Infrared intensity-weighting of the density-of-states (DOS) distributions by means of the squared dipole moment derivatives (which vary by a factor of 20 over the OH stretching frequency band at the B3LYP level), changes the downshift from approximately -205 to -275 cm(-1) at the B3LYP level. Explicit inclusion of relevant third-shell water molecules in the supermolecular cluster leads to a further downshift by approximately -30 cm(-1). Our final estimated average downshift is approximately -305 cm(-1). The experimental value lies somewhere in the range between -290 and -420 cm(-1). Also, the absolute nu(OH) frequency is well reproduced in our calculations. "In-liquid" instantaneous correlation curves between nu(OH) and various typical H-bond strength parameters such as R(O...O), R(H...O), the intramolecular OH bond length, and the IR intensity are presented. Some of these correlations are robust and persist also for the rather distorted instantaneous geometries in the liquid; others are less so.     

Modeling SN2 reactions in methanol solution by ab initio calculation of nucleophile solvent-substrate clusters

[Structure: see text]. Ab initio calculations were used to study the S(N)2 reactions of the CH3OCH2I molecule with a methoxide ion (CH3O-) and a methanol molecule by systematically building up the reaction system with explicit incorporation of the methanol solvent molecules. For the reaction of CH3OCH2I with a methoxide ion, the explicit incorporation of the methanol molecules to better solvate the methoxide ion led to an increase in the barrier to reaction. For the reaction of CH3OCH2I with a methanol molecule, the explicit incorporation of the methanol molecules led to a decrease in the barrier to reaction because of an inclination of this reaction to proceed with the nucleophilic displacements accompanied by proton transfer through the H-bonding chain. The H-bonding chain served as both acid and base catalysts for the displacement reaction. A ca. 10(15)-fold acceleration of the methanol tetramer incorporated S(N)2 reaction was predicted relative to the corresponding methanol monomer reaction. The properties of the reactions examined are discussed briefly.     

MgH^+,Mg2H^+量子化学从头计算的研究

用量子化学自洽场分子轨道从头计算方法,采用STO-3G基组计算了MgH~+,Mg_2H~+离子的位能曲线与位能面;给出电子波函数与电子集居数;得到MgH~+的平衡构型键长γ_(Mg-H)1.60A.Mg_2H~+的平衡构型有C_(∞v)与D_(∞h)两类,前者为[Mg—Mg—H]~+,键长γ_(Mg-Mg)2.41,γ_(Mg-H)1.63A;后者为[Mg—H—Mg]~+,γ_(Mg-H)1.73A.前者的总能量比后者低50kcal/mol.还讨论了它们的稳定性.MgH~+的键长计算结果和实测值较吻合,Mg_2H~+的平衡构型计算结果支持了Porter从热力学研究提出的假设.     

Communication: Thermodynamics of water modeled using ab initio simulations

We regularize the potential distribution framework to calculate the excess free energy of liquid water simulated with the BLYP-D density functional. Assuming classical statistical mechanical simulations at 350 K model the liquid at 298 K, the calculated free energy is found in fair agreement with experiments, but the excess internal energy and hence also the excess entropy are not. The utility of thermodynamic characterization in understanding the role of high temperatures to mimic nuclear quantum effects and in evaluating ab initio simulations is noted.     

Microscopic solvation of NaBO2 in water: anion photoelectron spectroscopy and ab initio calculations

We investigated the microscopic solvation of NaBO(2) in water by conducting photoelectron spectroscopy and ab initio studies on NaBO(2)(-)(H(2)O)(n) (n = 0-4) clusters. The vertical detachment energy (VDE) of NaBO(2)(-) is estimated to be 1.00 ± 0.08 eV. The photoelectron spectra of NaBO(2)(-)(H(2)O)(1) and NaBO(2)(-)(H(2)O)(2) are similar to that of bare NaBO(2)(-), except that their VDEs shift to higher electron binding energies (EBE). For the spectra of NaBO(2)(-)(H(2)O)(3) and NaBO(2)(-)(H(2)O)(4), a low EBE feature appears dramatically in addition to the features observed in the spectra of NaBO(2)(-)(H(2)O)(0-2). Our study shows that the water molecules mainly interact with the BO(2)(-) unit in NaBO(2)(-)(H(2)O)(1) and NaBO(2)(-)(H(2)O)(2) clusters to form Na-BO(2)(-)(H(2)O)(n) type structures, while in NaBO(2)(-)(H(2)O)(3) and NaBO(2)(-)(H(2)O)(4) clusters, the water molecules can interact strongly with the Na atom, therefore, the Na-BO(2)(-)(H(2)O)(n) and Na(H(2)O)(n)···BO(2)(-) types of structures coexist. That can be seen as an initial step of the transition from a contact ion pair (CIP) structure to a solvent-separated ion pair (SSIP) structure for the dissolution of NaBO(2).     

Theoretical study for the HNSi–HSiN isomerization: ab initio and DFT calculations

Ab initio and density functional theory (DFT) calculations using the GAUSSIAN 94 program have been performed to investigate the molecular structures of HNSi and HSiN in the ground state as well as the transition state for the HNSi–HSiN isomerization reaction at the 6-311G(d,p), 6-311+G(2d,p) and 6-311+G(2df,p) basis sets. The results show that DFT calculations at higher levels of theory reproduce experimental vibrational frequencies of both HNSi and HSiN better than ab initio methods including electron correlation effects. Those calculated geometries are accurate enough to predict the rotational constant of HNSi. The barrier height for the isomerization reaction is found to be about 10 kcal/mol.     

The singlet-triplet separation in CF2: state-of-the-art ab initio calculations and Franck-Condon simulations including anharmonicity.

Geometrical parameters, vibrational frequencies and relative electronic energies of the X2B1 state of CF2- and the X1A1 and ?3B1 states of CF2 have been calculated. Core-electron effects on the computed minimum-energy geometries and relative electronic energies have been investigated, and relativistic contributions to the computed relative electronic energies calculated. Potential energy functions of the X2B1 state of CF2- and the X1A1 and ?3B1 states of CF2 have been determined, and anharmonic vibrational wavefunctions of these states calculated variationally. Franck-Condon factors including anharmonicity and Duschinsky rotation have been computed and used to simulate the ?-X emission spectrum of CF2 determined by S. Koda [Chem. Phys. Lett. 1978, 55, 353] and the 364 nm laser photodetachment spectrum of CF2- obtained by R. L. Schwartz et al. [J. Phys. Chem. A 1999, 103, 8213]. Comparison between theory and experiment shows that the theoretical approach benchmarked in the present study is able to give highly reliable positions for the CF2(X1A1) + e <-- CF2-(X2B1) and CF2(?3B1) + e <-- CF2-(X2B1) bands in the photoelectron spectrum of CF2- and a reliable singlet-triplet gap for CF2. It is therefore concluded that the same theoretical approach should give reliable simulated CCl2(X1A1) + e <-- CCl2-(X2B1) and CCl2(?3B1) + e <-- CCl2-(X2B1) bands in the photodetachment spectrum of CCl2- and a reliable singlet-triplet gap for CCl2.     

Organolithium complex as a gas sensing material for oxides from ab initio calculations and molecular dynamics simulations

Gas sensing study of C₂H₄Li complex toward oxides viz. CO, CO₂, NO, NO₂, SO, and SO₂ gas molecules has been carried out using ab initio method. Different possible configurations of gas molecule adsorption on C₂H₄Li complex are considered. The structural parameters of most stable configuration of gas molecule adsorbed complexes are thoroughly analysed. Electronic properties are studied using total density of states (DOS) plot. Charge transferred between the gas molecule and the substrate is studied using NBO charge analysis. Gas sensing of all the six gas molecules is possible at ambient conditions. Atom centred density matrix propagation (ADMP) molecular dynamics simulations confirmed that all the gas molecules remain adsorbed on C₂H₄Li complex at room temperature during the simulation. This study suggests that the C₂H₄Li complex acts as a novel gas sensing material for CO, CO₂, NO, NO₂, SO, and SO₂ gas molecules at ambient conditions, below room temperature as well as at high pressure.     

DFT and ab initio study on the reaction mechanism of CH2SH + O2

The mechanism for the CH₂SH + O₂ reaction was investigated by DFT and ab initio chemistry methods. The geometries of all possible stationary points were optimized at the B3LYP/6-311+G(d,p) level, and the single point energy was calculated at the CCSD(T)/cc-pVXZ(X = D and T), G3MP2 and BMC-CCSD levels. The results indicate that the oxidation of CH₂SH by O₂ to form HSCH₂OO is a barrierless process. The most favorable channel is the rearrangement of the initial adduct HSCH₂OO (IM1) to form another intermediate H₂C(S)OOH (IM3) via a five-center transition state, and then the C–O bond fission in IM3 leads to a complex CH₂S. . .HO₂ (MC1), which finally gives out to the major product CH₂S + HO₂. Due to high barriers, other products including cis- and trans-HC(O)SH + HO could be negligible. The direct abstraction channel was also determined to yield CH₂S + HO₂, with the barrier height of 22.3, 18.1 and 15.0 kcal/mol at G3MP2, CCSD(T)/cc-pVTZ and BMC-CCSD levels, respectively, it is not competitive with the addition channel, in which all stationary points are lower than reactant energetically. The other channels to produce cis- and trans-CHSH + HO₂ are also of no importance.     

The interaction of water and dibromine in the gas phase: an investigation of the complex H(2)O...Br(2) by rotational spectroscopy and ab initio calculations

The ground-state rotational spectra of eight isotopomers of a complex formed by water and dibromine in the gas phase were observed by pulsed-jet, Fourier transform microwave spectroscopy. The spectroscopic constants B(0), C(0), delta(J), delta(JK), chi(aa)(Br(x)) (x=i for inner, o for outer), [chi(bb)(Br(x))-chi(cc)(Br(x))] and M(bb)(Br(x)) were determined for H(2)O...(79)Br(79)Br, H(2)O...(81)Br(79)Br, H(2)O...(79)Br(81)Br, H(2)O...(81)Br(81)Br, D(2)O...(79)Br(81)Br and D(2)O...(81)Br(81)Br. For the isotopomers HDO...(79)Br(81)Br and HDO...(81)Br(81)Br, only (B(0) + C(0))/2, delta(J), the chi(aa)(Br(x)) and M(bb)(Br(x)) were determinable. The spectroscopic constants were interpreted on the basis of several models of the complex to give information about its geometry, binding strength and the extent of electronic rearrangement on complex formation. The molecule H(2)O...Br(2) has C(s) symmetry with a pyramidal configuration at O. The zero-point effective quantities r(O...Br(i))=2.8506(1) A and phi(0)=46.8(1), where phi is the angle between the C(2) axis of H(2)O and the O...Br-Br internuclear axis, were obtained under the assumption of monomer geometries unchanged by complexation. Ab initio calculations, carried out at the aug-cc-pVDZ/MP2 level of theory, gave the equilibrium values r(e)(O...Br(i))=2.7908 A and phi(e)=45.7 degrees and confirmed the collinearity of the O...Br-Br nuclei. The potential energy function V(phi), also determined ab initio, showed that the wavenumber required for inversion of the configuration at O in the zero-point state is only 9 cm(-1). By interpreting the Br nuclear quadrupole coupling constants, the fractions delta(O-->Br(i))=0.004(5) and delta (Br(i)-->Br(o))=0.050(2) of an electron were determined to be transferred from O to Br(i) and Br(i) to Br(o), respectively, when the complex is formed. The complex is relatively weak, as indicated by the small value k(sigma)=9.8(2) N m(-1) of the intermolecular stretching force constant obtained from delta(J). A comparison of the properties, similarly determined, of H(2)O...F(2), H(2)O...Cl(2), H(2)O...Br(2), H(2)O...BrCl, H(2)O...ClF and H(2)O...ICl is presented.     

Probing mechanistic photochemistry of glyoxal in the gas phase by ab initio calculations of potential-energy surfaces and adiabatic and nonadiabatic rates

In the present work, the wavelength-dependent mechanistic photochemistry of glyoxal in the gas phase has been explored by ab initio calculations of potential-energy surfaces, surface crossing points, and adiabatic and nonadiabatic rates. The CHOCHO molecules in S1 by photoexcitation at 393-440 nm mainly decay to the ground state via internal conversion, which is followed by molecular eliminations to form CO, H2CO,H2, and HCOH. Upon photodissociation of CHOCHO at 350-390 nm, intersystem crossing to T1 followed by the C-C bond cleavage is the dominant process in this wavelength range, which is responsible for the formation of the CHO radicals. The C-C and C-H bond cleavages along the S1 pathway are energetically accessible upon photodissociation of CHOCHO at 290-310 nm, which can compete with the S1-->T1 intersystem crossing process. The present study predicts that the C-H bond cleavage on the S1 surface is probably a new photolysis pathway at high excitation energy, which has not been observed experimentally. In addition, the trans-cis isomerization is predicted to occur more easily in the ground state than in the excited states.     

Structure of the photodissociation products of CCl4, CBr4, and CI4 in solution studied by DFT and ab initio calculations

Various molecular species that can be populated during the photoreaction of carbon tetrahalides CX(4) (X = Cl, Br, I) in the gas phase and in solution have been studied by ab initio and density functional theory (DFT) calculations. Geometries, energies, and vibrational frequencies of CX(4), CX(3), CX(2), C(2)X(6), C(2)X(5), C(2)X(4), X(2), and the isomer X(2)CX-X were calculated and transition states connecting these species were characterized. Spin-orbit DFT (SODFT) computations were also performed to include the relativistic effects, which cannot be neglected for Br and I atoms. The calculated potential energy surfaces satisfactorily describe the reactions of the photoexcited CX(4) molecules. In the gas phase, the initial C-X bond rupture in CX(4) is followed by secondary C-X breakage in the CX(3) radical, leading to CX(2) and 2X, and the formation of C(2)X(6) or C(2)X(4) through bimolecular recombination of the CX(3) or CX(2) radicals is favored thermodynamically. In solution, by contrast, the X(2)CX-X isomer is formed via X-X binding, and two CX(3) radicals recombine nongeminately to form C(2)X(6), which then dissociates into C(2)X(4) and X(2) through C(2)X(5). The Raman intensities and the vibrational frequencies, as well as the absorption spectra and oscillator strengths of the Br(2)CBr-Br isomer in the gas phase and in various solvents were computed and the calculated absorption and Raman spectra of the Br(2)CBr-Br isomer in various solutions are in good agreement with the experimental data. The natural population analysis indicates that the Br(2)CBr-Br isomer corresponds to the recently reported solvent-stabilized solvated ion pair (CBr(3)(+)//Br(-))(solv) in the highly polar alcohol solvent. The singlet-triplet energy separations of the CX(2) radicals in the gas phase and in solution were evaluated with high level computational methods, and the optimized geometric parameters are in good agreement with the experimental results. The geometric and energetic differences between the singlet and triplet states were explained by the electronic properties of the CX(2) radicals. C(2)X(4), C(2)X(5), and C(2)X(6) (X = Br, I) in the gas phase and in solution were optimized at different computational levels, and the optimized geometric parameters of C(2)I(4) are in very good agreement with the experimental data.