Search PubMed⌕ Search

Biomedical subjects

Chia-Chung Sun

Publications and source records attributed to Chia-Chung Sun.

At least 37 records · Page 2Linked to original sources

Theoretical study on the [Si, C, N, O] potential energy surface.

The structures, energetics, spectroscopies, and stabilities of the doublet [Si, C, N, O] radical are explored at the density functional theory and ab initio levels. Sixteen isomers are located, connected by 29 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest lying isomer is a linear SiNCO 1 (0.0 kcal/mol) mainly featuring a cumulene | . Si = N = C = O. The second and third low-lying isomers are bent OSiCN 2 (8.8) and bent OSiNC 3 (11.1), respectively. All the three low-lying isomers 1, 2, 3, and another high-lying species 5 (75.4) with a linear SiCNO structure are shown to have considerable kinetic stability and may be experimentally observable. The predicted results of isomers 1 and 2 are consistent with the previous mass spectrometry experiments. Moreover, the fourth low-lying species SiOCN 4 (23.9) with bent structure is expected to be observable in low-temperature environments. The bonding nature of the five isomers 1, 2, 3, 4, and 5 is analyzed. The calculated results are compared with those of the analogous molecules C(2)NO and Si(2)NO. Implications in interstellar space and N,O-doped SiC vaporization processes are also discussed.

Journal Article↗

Special stability of cationic MPb12+ clusters and superalkali character of neutral MPb12 clusters (M = B, Al, Ga, In, and Tl).

The electronic structures and stabilities of cationic MPb12+ clusters (M = B, Al, Ga, In, and Tl) with 50 valence electrons are investigated within density functional theory. It is shown that, at the B3LYP/cc-pVDZ(-PP) and BPW91/cc-pVDZ(-PP) levels of theory, the structures of MPb12+ with icosahedra (I(h)) symmetry are energetically favorable, and their high stabilities may arise from the closed-shell nature of the pi subsystems which are subject to the 2(N(pi + 1)2 rule with N(pi = 1). In addition, the possessing of large nucleus-independent chemical shifts of the five kinds of clusters reflects the common aromatic character of these clusters. From the comparison of our studies on the binding energies and the highest occupied molecular orbital and the lowest unoccupied molecular orbital energy gaps, the cluster AlPb12+ has higher stability than the others and this is consistent with the recent mass-spectrometric discovery of Al-doped Pb(n)+ clusters, in which AlPb12+ is highly abundant. The same methods are used to search for the structures of the neutral MPb12 clusters. The calculations reveal that the most stable geometries of the BPb12 and GaPb12 clusters have I(h) symmetry, the AlPb12 and InPb12 clusters have T(h) symmetry, and the TlPb12 cluster has C5v symmetry. Furthermore, the vertical ionization potentials of the neutral MPb12 clusters are smaller than that of some alkali atoms, indicating that the neutral MPb12 clusters possess superalkali character.

Journal Article↗

Ylidic radical-molecule reactions.

The known reactions between a neutral radical and a neutral closed-shell sigma-type molecule are generally associated with the atom/group-abstraction (e.g., CN, OH reactions) or the carbenoid insertion processes (e.g., CH reactions). In this article, we describe a new type of neutral radical-molecule reaction that feature the nucleophilic addition and elimination mechanism based on the systematic Gaussian-3//MP2(full)/6-31G(d) investigations on the PH(3)CH+HX model reactions (X = CH(3), NH(2), OH, F, CH(2)F, CHF(2), CF(3), NHCH(3)/CH(2)NH(2), OCH(3)/CH(2)OH). The novel properties are found to benefit both from the radical and ylidic characters of PH(3)CH, and the type of reagent HX determines the nucleophilic addition reactivity. The present work is the first study on the chemical reactivity of the ylidic radicals.

Computer Simulation↗

Monosilicon-substituted cyanoacetylene: a computational study.

A detailed theoretical investigation of the [H,Si,C(2),N] potential energy surfaces including 28 minimum isomers and 65 interconversion transition states is reported at the Gaussian-3//B3LYP/6-31G(d) level. Generally, the triplet species lie energetically higher than the singlet ones. The former three low-lying isomers are linear HCCNSi 1 (0.00 kcal/mol), branched SiC(H)CN 12 (7.09 kcal/mol), and bent HNCCSi 7 (14.22 kcal/mol), which are separated by rather high barriers from each other and are kinetically very stable with the least conversion barriers of 32.6-70.5 kcal/mol. Two energetically high-lying isomers HCNCSi 3 (42.99 kcal/mol) and SiC(H)NC 13 (36.05 kcal/mol) are also kinetically stable with a barrier of 49.19 and 21.42 kcal/mol, respectively. Additionally, five high-lying isomers, that is, three chainlike isomers, HCCSiN 2 (55.17), HCSiNC 6 (47.80), HSiNCC 11 (78.83), and one three-membered ring isomer HN-cSiCC 19 (51.21), and one four-membered ring isomer cSiCN(H)C 27 (50.6 kcal/mol), are predicted to each have lower conversion barriers of 12-18 kcal/mol and can be considered as meta-stable species. All of the predicted 10 isomers could exist as stable or meta-stable intermediates under suitable conditions. Finally, the structural and bonding analysis indicate that the [H,Si,C(2),N] molecule contains various properties that are of chemical interest (e.g., silylene, SiC triple bonding, and conjugate SiN triple bonding and CC triple bonding, charge-transfer specie, planar aromatic specie, cumulate double bonding). This is the first detailed theoretical study on the potential energy surfaces of the series of hydrogenated Si,C,C,N-containing molecules. The knowledge of the present monohydrogenated SiC(2)N isomerism could provide useful information for more highly hydrogenated or larger Si,C(2),N-containing species.

Acetylene↗

Radical reaction C3H+NO: a mechanistic study.

Although a number of hydrocarbon radicals including the heavier C(3)-radicals C(3)H(3) and C(3)H(5) have been experimentally shown to deplete NO effectively, no theoretical or experimental attempts have been made on the reactivity of the simplest C(3)-radical towards NO. In this article, we report our detailed mechanistic study on the C(3)H+NO reaction at the Gussian-3//B3LYP/6-31G(d) level by constructing the singlet and triplet electronic state [H,C(3),N,O] potential energy surfaces (PESs). The l-C(3)H+NO reaction is shown to barrierlessly form the entrance isomer HCCCNO followed by the direct O-elimination leading to HCCCN+(3)O on triplet PES, or by successive O-transfer, N-insertion, and CN bond-rupture to generate the product (1)HCCN+CO on singlet PES. The possible singlet-triplet intersystem crossings are also discussed. Thus, the novel reaction l-C(3)H+NO can proceed effectively even at low temperatures and is expected to play an important role in both combustion and interstellar processes. For the c-C(3)H+NO reaction, the initially formed H-cCCC-NO can most favorably isomerize to HCCCNO, and further evolution follows that of the l-C(3)H+NO reaction. Quantitatively, the c-C(3)H+NO reaction can take place barrierlessly on singlet PES, yet it faces a small barrier 2.7 kcal/mol on triplet PES. The results will enrich our understanding of the chemistry of the simplest C(3)-radical in both combustion and interstellar processes, which to date have received little attention despite their importance and available abundant studies on its structural and spectroscopic properties.

Computer Simulation↗

Theoretical mechanistic study on the radical-molecule reaction of CHCl2/CCl3 with NO2.

The radical-molecule reaction mechanism of CHCl(2) and CCl(3) with NO(2) have been explored theoretically at the B3LYP/6-311G(d,p) and MC-QCISD (single-point) levels. For the singlet potential energy surface (PES) of CHCl(2) + NO(2) reaction, the association of CHCl(2) with NO(2) was found to be a barrierless carbon-to-nitrogen approach forming an energy-rich adduct a (HCl(2)CNO(2)) followed by isomerization to b(1) (trans-cis-HCl(2)CONO), which can easily interconvert to b(2), b(3), and b(4). Subsequently, the most feasible pathway is the 1,3-chlorine migration associated with N-O1 bond cleavage of b(1) leading to P(1) (CHClO + ClNO). The second competitive pathway is the 1,4-chlorine migration along with N-O1 bond rupture of b(4) giving rise to P(2) (CHClO + ClON). Moreover, some of P(1) and P(2) can further dissociate to give P(6) (CHClO + Cl + NO). The lesser followed competitive channel is the 1,3-H-shift from C to N atom along with N-O1 bond rupture of b(1) to form P(3) (CCl(2)O + HNO). The concerted 1,4-H-shift accompanied by N-O1 bond fission of b(3) to product P(4) (CCl(2)O + HON) is even much less feasible. For the singlet PES of CCl(3) + NO(2) reaction, the only primary product is found to be P(1) (CCl(2)O + ClNO), which can lead to P(2) (CCl(2)O + Cl + NO) via dissociation of ClNO. The obtained major products CHClO and CCl(2)O for CHCl(2) + NO(2) and CCl(3) + NO(2) reactions, respectively, are in good agreement with kinetic detection in experiment. Compared with the singlet pathways, the triplet pathways may have less contributions to both reactions. Because the rate-determining transition state involved in the feasible pathways lie above the reactants R, the title reactions may be important in high-temperature processes. The similarities and discrepancies among the CH(n)Cl(3-n) + NO(2) (n == 0-2) reactions are discussed in terms of the substitution effect. The present study may be helpful for further experimental investigation of the title reactions.

Chloroform↗

Theoretical study of Al(n) and Al(n)O (n = 2-10) clusters.

The stable structures, energies, and electronic properties of neutral, cationic, and anionic clusters of Al(n) (n = 2-10) are studied systematically at the B3LYP/6-311G(2d) level. We find that our optimized structures of Al5(+), Al9(+), Al9(-), Al10, Al10(+), and Al10(-) clusters are more stable than the corresponding ones proposed in previous literature reports. For the studied neutral aluminum clusters, our results show that the stability has an odd/even alternation phenomenon. We also find that the Al3, Al7, Al7(+), and Al7(-) structures are more stable than their neighbors according to their binding energies. For Al7(+) with a special stability, the nucleus-independent chemical shifts and resonance energies are calculated to evaluate its aromaticity. In addition, we present results on hardness, ionization potential, and electron detachment energy. On the basis of the stable structures of the neutral Al(n) (n = 2-10) clusters, the Al(n)O (n = 2-10) clusters are further investigated at the B3LYP/6-311G(2d), and the lowest-energy structures are searched. The structures show that oxygen tends to either be absorbed at the surface of the aluminum clusters or be inserted between Al atoms to form an Al(n-1)OAl motif, of which the Al(n-1) part retains the stable structure of pure aluminum clusters.

Journal Article↗

Theoretical mechanistic study on the radical-molecule reaction of CH2OH with NO2.

The complex singlet potential energy surface for the reaction of CH2OH with NO2, including 14 minimum isomers and 28 transition states, is explored theoretically at the B3LYP/6-311G(d,p) and Gaussian-3 (single-point) levels. The initial association between CH2OH and NO2 is found to be the carbon-to-nitrogen approach forming an adduct HOCH2NO2 (1) with no barrier, followed by C-N bond rupture along with a concerted H-shift leading to product P1 (CH2O + trans-HONO), which is the most abundant. Much less competitively, 1 can undergo the C-O bond formation along with C-N bond rupture to isomer HOCH2ONO (2), which will take subsequent cis-trans conversion and dissociation to P2 (HOCHO + HNO), P3 (CH2O + HNO2), and P4 (CH2O + cis-HONO) with comparable yields. The obtained species CH2O in primary product P1 is in good agreement with kinetic detection in experiment. Because the intermediate and transition state involved in the most favorable pathway all lie blow the reactants, the CH2OH + NO2 reaction is expected to be rapid, as is confirmed by experiment. These calculations indicate that the title reaction proceeds mostly through singlet pathways; less go through triplet pathways. In addition, a mechanistic comparison is made with the reactions CH3 + NO2 and CH3O + NO2. The present results can lead us to deeply understand the mechanism of the title reaction and may be helpful for understanding NO2-combustion chemistry.

Journal Article↗

Theoretical mechanistic study on the radical-radical reaction of ketenyl with nitrogen dioxide.

The radical-radical reaction between the ketenyl radical (HCCO) and nitrogen dioxide (NO(2)) played a very important role in atmospheric and combustion chemistry. Motivated by recent laboratory characterization about the reaction kinetics of ketenyl radical with nitrogen dioxide, in this contribution, we applied the coupled cluster and density functional theory to explore the mechanism of the title reaction. These calculations indicate that the title reaction proceeds mostly through singlet pathways, less go through triplet pathways. It is found that the HCCO + NO(2) reaction initially favors formation of adduct OCCHNO(2) (1) with no barrier. Subsequently, starting from isomer 1, the most feasible pathway is ring closure of 1 to isomer O-cCCHN(O)O (2) followed by CO(2) extrusion to product HCNO + CO(2) (P(1)), which is the major product with predominant yields. Much less competitively, 1 can take the successive 1,3-H- and 1,3-OH-shift interconversion to isomer OCCNOHO (3(a), 3(b), 3(c)) and then to isomer OCOHCNO (4(a), 4(b)), which can finally take a concerted H-shift and C-C bond fission to give HCNO + CO(2) (P(1)). The least competitive pathway is the ring-closure of isomer 3(a) to form isomer O-cCCN(OH)O (5(a), 5(b)) followed by dissociation to HONC + CO(2) (P(2)) through the direct side CO(2) elimination. Because the intermediates and transition states involved in the most favorable channel all lie below the reactants, the title reaction is expected to be rapid, as is confirmed by experiment. Therefore, it can be significant for elimination of nitrogen dioxide pollutants. The present results can lead us to a deep understanding of the mechanism of the title reaction and can be helpful for understanding NO(x)-combustion chemistry.

Journal Article↗

Characterization of solvated electrons in hydrogen cyanide clusters: (HCN)n- (n=3, 4).

Theoretical studies of the solvated electrons (HCN)n- (n=3, 4) reveal a variety of electron trapping possibilities in the (HCN)n (n=3, 4) clusters. Two isomers for (HCN)3- and four isomers for (HCN)4- are obtained at the MP2aug-cc-pVDZ+dBF (diffusive bond functions) level of theory. In view of vertical electron detachment energies (VDEs) at the CCSD(T) level, the excess electron always "prefers" locating in the center of the system, i.e., the isomer with higher coordination number shows larger VDE value. However, the most stable isomers of the solvated electron state (HCN)3- and (HCN)4- are found to be the linear Cinfinitynu and Dinfinityh structures, respectively, but not the fullyl symmetric structures which have the largest VDE values.

Journal Article↗

Nonlinear optical properties of alkalides Li+(calix[4]pyrrole)M- (M = Li, Na, and K): alkali anion atomic number dependence.

A new type of alkalide compound, Li+(calix[4]pyrrole)M- (M = Li, Na, and K), is presented in theory, which may be stable at room temperature. It has been shown by our calculations that the first hyperpolarizability (beta) is considerably large by means of the density functional theory method. The beta values are determined at the B3LYP/6-311++G level (for the alkali atoms the 6-311++G(3df) basis set is employed) as 8.9 x 103, 1.0 x 104, and 2.4 x 104 au for M = Li, Na, and K, respectively. These beta values are much larger than that of electride Li+(calix[4]pyrrole)e- (beta = 7.3 x 103 au) by a factor of 1.2 to 3.4. Comparing to the cryptand calix[4]pyrrole, the beta values of Li+(calix[4]pyrrole)M- are enhanced by 20-60 times. It is revealed, for the first time, that the beta value of alkalide compounds depends on the atomic number of the alkali anion, and it can be enhanced by choosing the akali anions with larger atomic numbers. The alkali anion in the alkalide compound decreases the transition energy and also increases the oscillator strength of the main transition, consequently the beta value is enhanced. This study proposes such a novel way to synthesize and design new NLO materials by using the alkali atom with a larger atomic number to create an anion in alkalide compounds.

Journal Article↗

The dependence of nanostructures on the molecule rigidity of A2(B4)2-type miktoarm block copolymer.

Using the dissipative particle dynamics simulation technique, we have studied the influence of the molecule rigidity on the nanostructures of the A2(B4)2-type miktoarm block copolymers. A typical spherical micellar ordered structure is obtained for a coil-coil miktoarm block copolymer in melt. By introducing a bond angle potential in our model to enhance the molecule rigidity systematically, we find, respectively, a hexagonal cylindrical structure and a parallel ellipsoid in lamellae structure which is discovered for the first time.

Journal Article↗

Role of each residue in catalysis in the active site of pyrimidine nucleoside phosphorylase from Bacillus subtilis: a hybrid QM/MM study.

Pyrimidine nucleoside phosphorylase (PYNP) catalyzes the reversible phosphorolysis of pyrimidines in the nucleotide synthesis salvage pathway. We have built a model of a closed active conformation of the three-dimensional structure of PYNP from Bacillus subtilis. Using docking, molecular dynamics, and hybrid quantum-mechanical/molecular-mechanical methods to study the reaction mechanics between PYNP and a substrate, we identified the role of each residue in the active site during the entire catalytic process. The results indicate that the function of His(82), Arg(169), and Lys(188) is to stabilize the uridine in a high-energy conformation by means of electrostatic interactions and that these residues are involved in catalysis. In addition, the function of Asp(162) is likely to activate Lys(188) for phosphorolytic catalysis through polarization effects.

Amino Acids↗

Gaseous reaction mechanism of C2F radical with water.

The kinetic properties of the carbon-fluorine radicals are little understood except those of CFn (n =1-3). In this article, a detailed mechanistic study was reported on the gas-phase reaction between the simplest pi-bonded C2F radical and water as the first attempt to understand the chemical reactivity of the C2F radical. Various reaction channels are considered. The most kinetically competitive channel is the quasi-direct hydrogen-abstraction route forming P5 HCCF + OH. At the CCSD(T)/6-311+G(2d,2p)//B3LYP/6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels, the overall H-abstraction barriers (4.5, 4.7, and 4.2 kcal/mol) for the C2F + H2O reaction are comparable to the corresponding values (5.5, 3.7, and 5.7 kcal/mol) for the analogous C2H + H2O reaction. This suggests that C2F is a reactive radical like the extensively studied C2H, in contrast to the situation of the CF and CF2 radicals that have much lower reactivity than the corresponding hydrocarbon species. Thus, the C2F radical is expected to play an important role in the combustion processes of the carbon-fluorine chemistry. Furthermore, addition of a second H2O can catalyze the reaction with the H-abstraction barrier significantly reduced to a marginally zero value (0.5 kcal/mol). This is also indicative of the potential relevance of the title reactions in the low-temperature atmospheric chemistry.

Journal Article↗

Theoretical study on structures and stabilities of [H,Ge,C,N].

Theoretical investigations are performed for the first time on the simplest hydrogenated germanium cyanide [H,Ge,C,N], whose analogs [H,C(2),N] and [H,Si,C,N] have been detected in space and laboratory, respectively. The detailed potential energy surfaces in both singlet and triplet states are constructed at the CCSD(T)/6-311+G(3df,2p)//B3LYP/6-31G(d)+ZPVE level, including 18 minimum isomers and 26 interconversion transition states. The former three low-lying and kinetically stabilized isomers are HGeCN (1)1 (0.0 kcal/mol), HGeNC (1)2 (5.1 kcal/mol), and cyclic cCHNGe(1)7 (11.1 kcal/mol). In addition, five isomers HCNGe (1)3 (33.8), HNCGe (1)5 (29.8), cNHCGe (1)8 (37.9), HGeCN (3)1 (30.1), and HNCGe (3)5 (26.5) each have considerable barriers, despite their high energies. Future laboratory characterization and astrophysical detection of the eight [H,Ge,C,N] isomers, especially the former three low-lying species (1)1, (1)2, and (1)7, are highly recommended. The accurate spectroscopic data at the QCISD/6-311G(d,p) level are provided. For some species, the CBS-QB3 calculations are also performed. Wherever possible, comparisons with the analogous [H,C(2),N] and [H,Si,C,N] are made on the structural, energetic, and bonding properties.

Chemical Phenomena↗

Reaction mechanism of the CCN radical with nitric oxide.

To investigate the possibility of the carbyne radical CCN in removal of nitric oxide, a detailed computational study is performed at the Gaussian-3//B3LYP/6-31G(d) level on the CCN + NO reaction by constructing the singlet and triplet electronic state [C(2)N(2)O] potential energy surfaces (PESs). The barrierless formation of the chain-like isomers NCCNO (singlet at -106.5, triplet cis at -48.2 and triplet trans at -47.6 kcal/mol) is the most favorable entrance attack on both singlet and triplet PESs. Subsequently, the singlet NCCNO takes an O-transfer to form the branched intermediate singlet NCC(O)N (-85.6), which can lead to the fragments CN + NCO (-51.2) via the intermediate singlet NCOCN (-120.3). The simpler evolution of the triplet NCCNO is the direct N-O rupture to form the weakly bound complex triplet NCCN...O (-56.2) before the final fragmentation to NCCN + (3)O (-53.5). However, the lower lying products (3)NCN + CO (-105.6) and (3)CNN + CO (-74.6) are kinetically much less competitive. All the involved transition states for generation of CN + NCO and NCCN + (3)O lie much lower than the reactants. Thus, the novel reaction CCN + NO can proceed effectively even at low temperatures and is expected to play a role in both combustion and interstellar processes. Significant differences are found on the singlet PES between the CCN + NO and CH + NO reaction mechanisms.

Journal Article↗

Theoretical study on the mechanism of the CH2F + NO2 reaction.

Despite the importance of the Fluoromethyl radicals in combustion chemistry, very little experimental information on their reactions toward stable molecules is available in the literature. Motivated by recent laboratory characterization about the reaction kinetics of Chloromethyl radicals with NO2, we carried out a detailed potential energy survey on the CH2F + NO2 reaction at the B3LYP/6-311G(d,p) and MC-QCISD (single-point) levels as an attempt toward understanding the CH2F + NO2 reaction mechanism. It is shown that the CH2F radical can react with NO2 to barrierlessly generate adduct a (H2FCNO2), followed by isomerization to b1 (H2FCONO-trans) which can easily interconvert to b2 (H2FCONO-cis). Subsequently, Starting from b (b1, b2), the most feasible pathway is the C--F and N--O1 bonds cleavage along with N--F bond formation of b (b1, b2) leading to P1 (CH2O + FNO), or the direct N--O1 weak-bond fission of b (b1, b2) to give P2 (CH2FO + NO), or the 1,3-H-shift associated with N--O1 bond rupture of b1 to form P3 (CHFO + HNO), all of which may have comparable contribution to the reaction CH2F + NO2. Much less competitively, b2 either take the 1,4-H-shift and O1--N bond cleavage to form product P4 (CHFO + HON) or undergo a concerted H-shift to isomer c2 (HFCONOH), followed by dissociation to P4. Because the rate-determining transition state (TSab1) in the most competitive channels is only 0.3 kcal/mol higher than the reactants in energy, the CH2F + NO2 reaction is expected to be rapid, and may thus be expected to significantly contribute to elimination of nitrogen dioxide pollutants. The similarities and discrepancies among the CH2X + NO2 (X = H, F, and Cl) reactions are discussed in terms of the electronegativity of halogen atom. The present article may assist in future experimental identification of the product distributions for the title reaction, and may be helpful for understanding the halogenated methyl chemistry.

Journal Article↗

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface.

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid-liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. [Figure: see text].

Carbon↗