Microsolvation effects on the electron capturing ability of thymine: thymine-water clusters

The effects of solvation on the stability of thymine and its negative ion have been investigated by explicitly considering the structures of complexes of thymine with up to five water molecules and the respective anions at the B3LYP/DZP++ level of theory. The vertical detachment energy of thymine was predicted to increase gradually with the hydration number, consistent with experimental observations from a photodetachment-photoelectron spectroscopy study J. Schiedt et al., [Chem. Phys. 239, 511 (1998)]. The adiabatic electron affinity of thymine was also found to increase with the hydration number, which implies that while the conventional valence anion of thymine is only marginally bound in the gas phase, it may form a stable anion in aqueous solution.     

Photoelectron spectrum of valence anions of uracil and first-principles calculations of excess electron binding energies

The photoelectron spectrum (PES) of the uracil anion is reported and discussed from the perspective of quantum chemical calculations of the vertical detachment energies (VDEs) of the anions of various tautomers of uracil. The PES peak maximum is found at an electron binding energy of 2.4 eV, and the width of the main feature suggests that the parent anions are in a valence rather than a dipole-bound state. The canonical tautomer as well as four tautomers that result from proton transfer from an NH group to a C atom were investigated computationally. At the Hartree-Fock and second-order Moller-Plesset perturbation theory levels, the adiabatic electron affinity (AEA) and the VDE have been converged to the limit of a complete basis set to within +/-1 meV. Post-MP2 electron-correlation effects have been determined at the coupled-cluster level of theory including single, double, and noniterative triple excitations. The quantum chemical calculations suggest that the most stable valence anion of uracil is the anion of a tautomer that results from a proton transfer from N1H to C5. It is characterized by an AEA of 135 meV and a VDE of 1.38 eV. The peak maximum is as much as 1 eV larger, however, and the photoelectron intensity is only very weak at 1.38 eV. The PES does not lend support either to the valence anion of the canonical tautomer, which is the second most stable anion, and whose VDE is computed at about 0.60 eV. Agreement between the peak maximum and the computed VDE is only found for the third most stable tautomer, which shows an AEA of approximately -0.1 eV and a VDE of 2.58 eV. This tautomer results from a proton transfer from N3H to C5. The results illustrate that the characteristics of biomolecular anions are highly dependent on their tautomeric form. If indeed the third most stable anion is observed in the experiment, then it remains an open question why and how this species is formed under the given conditions.     

Effects of microsolvation on uracil and its radical anion: uracil(H2O)n (n = 1-5)

Microsolvation effects on the stabilities of uracil and its anion have been investigated by explicitly considering the structures of complexes of uracil with up to five water molecules at the B3LYPDZP++ level of theory. For all five systems, the global minimum of the neutral cluster has a different equilibrium geometry from that of the radical anion. Both the vertical detachment energy (VDE) and adiabatic electron affinity (AEA) of uracil are predicted to increase gradually with the number of hydrating molecules, qualitatively consistent with experimental results from a photodetachment-photoelectron spectroscopy study [J. Schiedt et al., Chem. Phys. 239, 511 (1998)]. The trend in the AEAs implies that while the conventional valence radical anion of uracil is only marginally bound in the gas phase, it will form a stable anion in aqueous solution. The gas-phase AEA of uracil (0.24 eV) was higher than that of thymine by 0.04 eV and this gap was not significantly affected by microsolvation. The largest AEA is that predicted for uracil(H2O)5, namely, 0.96 eV. The VDEs range from 0.76 to 1.78 eV.     

Barrier-free proton transfer in the valence anion of 2'-deoxyadenosine-5'-monophosphate. II. A computational study

The propensity of four representative conformations of 2(')-deoxyadenosine-5(')-monophosphate (5(')-dAMPH) to bind an excess electron has been studied at the B3LYP6-31++G(d,p) level. While isolated canonical adenine does not support stable valence anions in the gas phase, all considered neutral conformations of 5(')-dAMPH form adiabatically stable anions. The type of an anionic 5(')-dAMPH state, i.e., the valence, dipole bound, or mixed (valence/dipole bound), depends on the internal hydrogen bond(s) pattern exhibited by a particular tautomer. The most stable anion results from an electron attachment to the neutral syn-south conformer. The formation of this anion is associated with a barrier-free proton transfer triggered by electron attachment and the internal rotation around the C4(')-C5(') bond. The adiabatic electron affinity of the a_south-syn anion is 1.19 eV, while its vertical detachment energy is 1.89 eV. Our results are compared with the photoelectron spectrum (PES) of 5(')-dAMPH(-) measured recently by Stokes et al., [J. Chem. Phys. 128, 044314 (2008)]. The computational VDE obtained for the most stable anionic structure matches well with the experimental electron binding energy region of maximum intensity. A further understanding of DNA damage might require experimental and computational studies on the systems in which purine nucleotides are engaged in hydrogen bonding.     

Temporary anion states and dissociative electron attachment in diphenyl disulfide

The temporary anion states of gas-phase diphenyl disulfide are characterized by means of electron transmission (ET) and dissociative electron attachment (DEA) spectroscopies. The measured energies of vertical electron attachment are compared to the virtual orbital energies of the neutral state molecule supplied by MP2 and B3LYP calculations with the 6-31G basis set. The calculated energies, scaled with empirical equations, reproduce satisfactorily the attachment energies measured in the ET spectrum. The first anion state of diphenyl disulfide is stable, thus escaping detection in ETS. The vertical and adiabatic electron affinities, evaluated with B3LYP/6-31+G calculations as the energy difference between the neutral and anion states, are predicted to be 0.37 and 1.38 eV, respectively. The anion current displayed in the DEA spectrum has a sharp and intense peak at zero energy, essentially due to the C6H5S- negative fragment. In agreement, according to the calculations, the localization properties of the first anion state are strongly S-S antibonding, and the energetic requirement for its dissociation along the S-S bond is fulfilled even at zero energy.     

Extrapolation to the limit of a complete basis set for electronic structure calculations on the N2 molecule

Results obtained from nonrelativistic electronic structure calculations using finite Gaussian basis sets are extrapolated to the limit of a complete basis set, employing the results of explicitly correlated coupled-cluster calculations including singles and doubles substitutions (CCSD). For N₂, the basis-set limits for the electronic binding energy, equilibrium bond length and harmonic vibrational wave number are established for the CCSD model including a perturbative correction for triples substitutions and for the internally contracted multireference configuration interaction method. The resulting numbers are in good agreement with experimental values. Received: 2 December 1997 / Accepted: 3 February 1998 / Published online: 17 June 1998     

Valence anion of thymine in the DNA pi-stack

Most of theoretical data on the stability of radical anions supported by nucleic acid bases have been obtained for anions of isolated nucleobases, their nucleosides, or nucleotides. This approach ignores the hallmark forces of DNA, namely, hydrogen bonding and pi-stacking interactions. Since these interactions might be crucial for the electron affinities of nucleobases bound in DNA, we report for the first time on the stability of the thymine valence anion in trimers of complementary bases possessing the regular B-DNA geometry but differing in base sequence. In order to estimate the energetics of electron attachment to a trimer, we developed a thermodynamic cycle employing all possible two-body interaction energies in the neutral and anionic duplex as well as the adiabatic electron affinity of isolated thymine. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. The two-body interaction energies were corrected for the basis set superposition error, and in benchmark systems, they were extrapolated to the basis set limit and supplemented with correction for higher order correlation terms calculated at the CCSD(T) level. We have demonstrated that the sequence of nucleic bases has a profound effect on the stability of the thymine valence anion: the anionic 5'-CTC-3' (6.0 kcal/mol) sequence is the most stable configuration, and the 5'-GTG-3' (-8.0 kcal/mol) trimer anion is the most unstable species. On the basis of obtained results, one can propose DNA sequences that are different in their vulnerability to damage by low energy electron.     

Theoretical study of the spectroscopically relevant parameters for the detection of HNP(q) and HPN(q) (q = 0, +1, -1) in the gas phase

High level ab initio electronic structure calculations at different levels of theory have been performed on HNP and HPN neutrals, anions, and cations. This includes standard coupled cluster CCSD(T) level with augmented correlation-consistent basis sets, internally contacted multi-reference configuration interaction, and the newly developed CCSD(T)-F12 methods in connection with the explicitly correlated basis sets. Core-valence correction and scalar relativistic effects were examined. We present optimized equilibrium geometries, harmonic vibrational frequencies, rotational constants, adiabatic ionization energies, electron affinities, vertical detachment energies, and relative energies. In addition, the three-dimensional potential energy surfaces of HNP(-1,0,+1) and of HPN(-1,0,+1) were generated at the (R)CCSD(T)-F12b∕cc-pVTZ-F12 level. The anharmonic terms and fundamentals were derived using second order perturbation theory. For HNP, our best estimate for the adiabatic ionization energy is 7.31 eV, for the adiabatic electron affinity is 0.47 eV. The higher energy isomer, HPN, is 23.23 kcal∕mol above HNP. HPN possesses a rather large adiabatic electron affinity of 1.62 eV. The intramolecular isomerization pathways were computed. Our calculations show that HNP(-) to HPN(-) reaction is subject to electron detachment.     

Valence and dipole-bound anions of the most stable tautomers of guanine

Anionic states of guanine, which is the only nucleic acid base of which the anions have not yet been studied in either photoelectron spectroscopic (PES) or Rydberg electron transfer (RET) experiments, have been characterized for the four most stable tautomers of neutral guanine using a broad spectrum of electronic structure methods from the density functional theory, with the B3LYP exchange-correlation functional, to the coupled-cluster method, with single, double, and perturbative triple excitations. Both valence and dipole-bound anionic states were addressed. We identified some of the difficulties facing future PES or RET experiments on the anion of guanine. Even if guanine is successfully transferred to the gas phase without thermal decomposition, it is critical to have the canonical amino-oxo (G) and both amino-hydroxy (GH and GHN7H) tautomers in the beam, not only the most stable, a noncanonical, amino-oxo tautomer (GN7H), as the latter does not support an adiabatically bound anionic state. We also suggested a scheme for enrichment of gas-phase guanine with the canonical tautomer, which is not the most stable in the gas phase, but which is of main interest due to its biological relevance. The tautomers G, GN7H, and GHN7H support vertically bound valence anionic states with the CCSD(T) value of vertical detachment energy of +0.58, +0.21, and +0.39 eV, respectively. These anionic states are, however, adiabatically unbound and thus metastable. The vertical electronic stability of these valence anionic states is accompanied by serious "buckling" of the molecular skeleton. The G and GHN7H tautomers support dipole-bound states with the CCSD(T) values of adiabatic electron affinity of 65 and 36 meV, respectively. A contribution from higher-than-second-order correlation terms represents, respectively, 48 and 68% of the total vertical electron detachment energy determined at the CCSD(T) level.     

Coupled-cluster methods with perturbative inclusion of explicitly correlated terms: a preliminary investigation

We propose to account for the large basis-set error of a conventional coupled-cluster energy and wave function by a simple perturbative correction. The perturbation expansion is constructed by L?wdin partitioning of the similarity-transformed Hamiltonian in a space that includes explicitly correlated basis functions. To test this idea, we investigate the second-order explicitly correlated correction to the coupled-cluster singles and doubles (CCSD) energy, denoted here as the CCSD(2)(R12) method. The proposed perturbation expansion presents a systematic and easy-to-interpret picture of the "interference" between the basis-set and correlation hierarchies in the many-body electronic-structure theory. The leading-order term in the energy correction is the amplitude-independent R12 correction from the standard second-order M?ller-Plesset R12 method. The cluster amplitudes appear in the higher-order terms and their effect is to decrease the basis-set correction, in accordance with the usual experience. In addition to the use of the standard R12 technology which simplifies all matrix elements to at most two-electron integrals, we propose several optional approximations to select only the most important terms in the energy correction. For a limited test set, the valence CCSD energies computed with the approximate method, termed , are on average precise to (1.9, 1.4, 0.5 and 0.1%) when computed with Dunning's aug-cc-pVXZ basis sets [X = (D, T, Q, 5)] accompanied by a single Slater-type correlation factor. This precision is a roughly an order of magnitude improvement over the standard CCSD method, whose respective average basis-set errors are (28.2, 10.6, 4.4 and 2.1%). Performance of the method is almost identical to that of the more complex iterative counterpart, CCSD(R12). The proposed approach to explicitly correlated coupled-cluster methods is technically appealing since no modification of the coupled-cluster equations is necessary and the standard M?ller-Plesset R12 machinery can be reused.     

Explicit and implicit solvation of radical ions: the cycloreversion of CPD dimers

The electron transfer catalyzed cycloreversion of cyclobutane pyrimidine dimers is the key step in repair of light-induced DNA lesions catalyzed by the enzyme CPD photolyase. The formation of the CPD radical anion was found to be strongly solvent dependent due to a specific hydrogen bond that stabilizes the valence bound state over the dipole bound state of the additional electron. The effect of solvation on the vertical and adiabatic electron affinity of uracil and uracil dimers as well as on the mechanism of the cycloreversion of the uracil dimer radical anion is explored for three model systems that include explicit solvent molecules at the B3LYP/6-311++G^∗∗/B3LYP/6-31+G^∗ level of theory. The second solvation shell is described using the implicit C-PCM solvation model. These calculations indicate an effectively barrierless mechanism. These results are in agreement with the available experimental data for the reaction energies and isotope effects. It is also shown that a single hydrogen bond donor is a sufficient minimal model for the first solvation shell by adequately describing the stabilization of the valence bound state of the radical anion through hydrogen bonding. The relationship of these model systems with the enzymatic reaction catalyzed by DNA photolyase is also discussed.     

Stability and spectroscopic properties of singly and doubly charged anions

Using density functional theory and hybrid B3LYP exchange-correlation energy functional we have studied the structure, stability, and spectroscopic properties of singly and doubly charged anions composed of simple metal atoms (Na, Mg, Al) decorated with halogens such as Cl and pseudohalogens such as CN. Since pseudohalogens mimic the chemistry of halogen atoms, our objective is to see if pseudohalogens can also form superhalogens much as halogens do and if the critical size for a doubly charged anion depends upon the ligand. The electron affinities of MCl(n) (M = Na, Mg, Al) exceed the value of Cl for n ≥ (k + 1), where k is the normal valence of the metal atom. However, for M(CN)(n) complexes this is only true when n = k + 1. In addition, while the electron affinities and vertical detachment energies of MCl(n) complexes are close to each other, they are markedly different when Cl is replaced by pseudohalogen, CN. The origin of these anomalous results is found to be due to the large binding energy of cyanogen, (NCCN) molecule. Because of the tendency of CN molecules to dimerize, the ground state geometries of the neutral and anionic M(CN)(n) complexes are very different when their number exceed the normal valence of the metal atom. While our calculations support the conclusion of Skurski and co-workers that pseudohalogens can form the building blocks of superhalogens, we show that there is a limitation on the number of CN moieties where this is true. Equally important, we find large differences between the ground state geometries of the neutral and anionic M(CN)(n) complexes for n ≥ (k + 2) which could play an important role in interpreting future experimental data on M(CN)(n) complexes. This is because the electron affinity defined as the energy difference between the ground states of the anion and neutral can be very different from the adiabatic detachment energy defined as the energy difference between the ground state of the anion and its structurally similar neutral isomer.     

A G3 study of the structure of carbon-nitrogen nanoclusters

Possible structures of the carbon-nitrogen clusters of the form C(m)N(n) (m = 1-4, n = 1-4, m + n = 2-5) were predicted for the neutral, anion, and cation species in the singlet, doublet, and triplet states, whenever appropriate. The calculations were performed at the G3, MP2(fc)/6-311+G*, and B3LYP/6-311+G* levels of theory. Several molecular properties related to the experimental data--such as the electronic energy, equilibrium geometry, binding energy, HOMO-LUMO gap (HLG), and spin contamination --were calculated. In addition the vertical electron attachment, the adiabatic electron affinity, and vertical ionization energy, of the neutral clusters were calculated. Most of the predicted lowest energy structures were linear, whereas bent structures became more stable with the increase of the cluster size and increase of the number of the N atoms. In most of the predicted lowest energy structures, the N atom prefers the terminal position with acetylenic bond. The calculated BE of the predicted clusters increases with the increase of the cluster size for the neutral and cation clusters but decreases with the increase of the cluster size for the anion clusters. The predicted clusters are characterized by high HLG of about 11 eV on the average, with that of the anion clusters is smaller than that for the neutral and cation clusters. It is concluded then that the anion clusters are less stable than the corresponding neutral and cation clusters. Finally, the N(2) loss reaction is treated.     

Electron attachment to dye-sensitized solar cell components: cyanoacetic acid

The energies of electron attachment associated with temporary occupation of the lower-lying virtual orbitals of cyanoacetic acid (CAA), proposed as a possible component of dye-sensitized solar cells, and its derivative methyl cyanoacetate (MCA) are measured in the gas phase with electron transmission spectroscopy (ETS). The corresponding orbital energies of the neutral molecule, supplied by B3LYP/6-31G(d) calculations and scaled using an empirically calibrated linear equation, are compared with the experimental vertical attachment energies (VAEs). The vertical and adiabatic electron affinities are also evaluated at the B3LYP/6-31+G(d) level as the anion/neutral total energy difference. Dissociative electron attachment spectroscopy (DEAS) is used to measure the total anion current as a function of the incident electron energy in the 0-4 eV energy range, and the negative fragments generated through the dissociative decay channels of the molecular anion are detected with a mass filter. In both compounds only two intense fragment anion currents are observed, that due to loss of a hydrogen atom from the molecular anion ([M - H](-)) and that due to formation of CN(-). In CAA the former signal displays a very sharp feature at 0.68 eV, assigned to a vibrational Feshbach resonance arising from coupling between a dipole bound anion state and a temporary σ* anion state.     

Highly accurate equilibrium structure of the C2h symmetric N1‐to‐O2 hydrogen‐bonded uracil‐dimer

The highly accurate ab initio equilibrium geometry of the hydrogen‐bonded uracil dimer is derived using a composite geometry extrapolation scheme based on all‐electron, complete basis set extrapolated Møller–Plesset perturbation theory using the jun‐pwCV[T,Q]Z basis sets combined with a valence CCSD(T)/cc‐pVTZ high‐level correction. Geometrical changes on dimerization are discussed and the performance of the several density functional approximations (among others SCAN, ωB97M‐V, DSD‐PBEP86‐D3(BJ), and DSD‐PBEP86‐NL) is evaluated. Orbital‐optimized MP2.5 is discussed as a reduced‐cost alternative to the CCSD(T) gradient in the composite scheme. A new reference interaction energy is calculated with explicitly correlated F12‐CCSD(T).     

A combined ab initio and Franck-Condon factor simulation study on the photodetachment spectrum of ScO2(-)

Restricted-spin coupled-cluster single-double plus perturbative triple excitation {RCCSD(T)} potential energy functions (PEFs) of the X(2)B2 state of ScO2 and the 1A1 state of ScO2(-) were computed, employing the augmented correlation-consistent polarized-weighted core-valence quadruple-zeta (aug-cc-pwCVQZ) basis set for Sc and augmented correlation-consistent polarized valence quadruple-zeta (aug-cc-pVQZ) basis set for O, and with the outer core Sc 3s(2)3p(6) electrons being explicitly correlated. Franck-Condon factors, which include allowance for Duschinsky rotation and anharmonicity, were calculated using the computed RCCSD(T) PEFs, and were used to simulate the first photodetachment band of ScO2(-). The simulated spectrum matches well with the corresponding experimental 355 nm photodetachment spectrum of Wu and Wang, J Phys Chem A 1998, 102, 9129, confirming the assignment of the photodetachment spectrum and the reliability of the RCCSD(T) PEFs used. Further calculations on low-lying electronic states of ScO2 gave adiabatic relative electronic energies (T(e)'s) of, and vertical excitation energies (T(v)'s) to, the 2A1, 2B1, and 2A2 states of ScO2 (from the X(2)B2 state of ScO2), as well as electron affinities (EAs) and vertical detachment energies (VDEs) to these neutral states from the 1A1 state of ScO2(-).     

Ab initio molecular-orbital study of structures and energetics of Si3H3 neutral and anion

The geometric structures and isomeric stabilities of various stationary points in Si(3)H(3) neutral and its anion are investigated at the coupled-cluster singles, doubles (triples) [CCSD(T)] level of theory. For geometrical surveys, the basis sets used are of the Dunning's correlation consistent basis sets of triple-zeta quality for the neutral. To the anions, the Dunning's correlation consistent basis sets of double-zeta quality with diffuse functions are applied. For the three lower-lying anion isomers, the Dunning's correlation consistent basis sets of triple-zeta quality with diffuse functions (aug-cc-pVTZ) are also used. The final energies for the optimized stationary points are calculated at the CCSD(T) level of theory with the aug-cc-pVTZ basis sets. The basis sets of 6-311++G(3df,2pd) were also used for the lower-lying anion isomers. The Gaussian-2 method was performed only for the lower-lying anion isomers to clarify the relative stabilities. The global minimum neutral 1 (C(1):(2)A) has an unsymmetrical hydrogen-bridged bond; the conformer 2 in C(s) symmetry is a saddle point connecting the two equivalent isomers 1. Two lower-lying isomers (3 and 4) are also predicted within the energy range of 20 kJmol. In the anion, however, the conformer 4 (C(s):(1)A(')) with five formal valence electrons is a global minimum. Two more isomers (2 and 3) lie within 20 kJmol as in the neutral; the conformer 1 converts to the isomer 2. The quartets for the neutrals and diradical triplets for the anions were further studied; lower-lying quartets and triplets, competing with the corresponding doublet and singlet, respectively, were not found in the present systems. The vertical and adiabatic electron affinities of the global minimum neutral 1, producing the second lowest-lying anion isomer 2, amount to 2.18 and 2.35 eV, respectively, at the CCSD(T)/aug-cc-pVTZ level of theory. The electron addition to the third lowest-lying neutral isomer 4 produces the largest vertical electron affinities of 2.48 eV. The D(3h) structure, being the global minimum in the corresponding Si(3)H(3) (+) cation (trisilacyclopropenyl cation), converts to the isomer 8 (C(s)) or 11 (C(2)) due to the Jahn-Teller effect in the Si(3)H(3) neutral.     

Accurate ab initio potential energy surface and vibration‐rotation energy levels of hydrogen peroxide

The accurate ground‐state potential energy surface of hydrogen peroxide, H₂O₂, has been determined from ab initio calculations using the coupled‐cluster approach in conjunction with the correlation‐consistent basis sets up to septuple‐zeta quality. Results obtained with the conventional and explicitly correlated coupled‐cluster methods were compared. The core–electron correlation, scalar relativistic, and higher‐order valence–electron correlation effects were taken into account. The adiabatic effects were also discussed. The vibration–rotation energy levels of the H₂O₂, D₂O₂, and HOOD isotopologues were predicted, and the experimental vibrational fundamental wavenumbers were reproduced to 1 cm^?1 (“spectroscopic”) accuracy. © 2012 Wiley Periodicals, Inc.     

Singlet-triplet energy splitting and excited states of phenylnitrene

The vertical and adiabatic singlet-triplet energy splittings (Delta E ST) of phenylnitrene were computed by a variety of multireference configuration interaction and perturbation theory methods employing basis sets of up to quadruple-xi quality and extrapolation to the complete basis set limit. The vertical and adiabatic energy gaps are 18.9 and 15.9 kcal mol (-1), respectively, the latter in reasonable agreement with the revised experimental value of 15.1 +/- 0.2 kcal mol (-1). The energy difference between both states at the geometry of the a (1)A 2 singlet state was also considered and amounts to 13.8 kcal mol (-1). In obtaining accurate state energy splittings, basis set completeness turns out to be a more important issue than the level of dynamical electron correlation treatment. Density functional theory that is frequently employed to investigate phenylnitrenes and their rearrangements yields varying results and, depending on the functional, gives adiabatic energy differences between 9 and 16 kcal mol (-1). The b (1)A 1 state has a similar geometry as the ground state of 1 and is 31 kcal mol (-1) higher in energy. According to best estimates, the next higher singlet states, c (1)A 1 and d (1)B 1, are 57 and 72 kcal mol (-1) above the ground state. In the triplet manifold, vertical excitation energies to the A (3)B 1 and B (3)A 2 states are 71 and 77 kcal mol (-1), respectively.     

An ab initio MO study on structures and energetics of CH2Si2, CH2Si2, and CH2Si2

The geometric structures and isomeric stabilities of various stationary points in CH₂Si₂ neutral, cation and anion are investigated at the coupled-cluster singles, doubles (triples) (CCSD(T)) level of theory. For the geometrical survey, the basis sets used are of the cc-pVTZ for the neutral and cation. The final energies are calculated by the use of the CCSD(T) level of theory with the aug-cc-pVTZ basis set at their optimized geometries. To the competitive two-anion isomers, the aug-cc-pVTZ basis sets are applied. The global minimum (N-1) of the CH₂Si₂ neutral has a quite different framework from those of the C₃H₂ (cyclopropenylidene) and Si₃H₂ (trisilacyclopropenylidene) neutrals. No competitive low-lying isomers are found in the CH₂Si₂ neutral. The attractive conformer (C-1) is predicted for the most stable cation, where its framework is quite different from that of the neutral N-1. Both H atoms are connected to the same C atom, but each C–H bond length is different from each other. Two competitive anion isomers with positive (real) electron affinities are predicted. The framework of the most stable anion A-1 is quite similar to that of the cation C-1, whereas both H atoms are equally connected to the same C atom. The framework of the anion isomer A-2 is the same as that in the neutral N-1. The vertical and adiabatic ionization potentials from the most stable neutral N-1 are 9.02 and 8.71 eV, respectively. The adiabatic electron affinity of the lowest lying isomer N-1 is only 0.43 eV and the vertical electron detachment energy form the global minimum anion (A-1) is 2.02 eV. The multi-centered Si–H–Si bonds are found in the neutral, cation, and anion.