• Bulletin of the Korean Chemical Society
• /
• v.28 no.12
• /
• pp.2454-2458
• /
• 2007

Exploring the basic concepts for the design of CO2-philic molecules is important due to the possibility for “green” chemistry in supercritical CO2 as substitute solvent systems. The Lewis acid-base interactions and C?H…O weak hydrogen bonding were suggested as two key factors for the solubility of CO2-philic molecules. We have performed high level quantum mechanical calculations for the van der Waals complexes of CO2 with trimethylphosphate and trimethylphosphine oxide, which have long been used for metal extractants in supercritical CO2 fluid. Structures and energies were calculated using the MP2/6-31+G(d) and recently developed multilevel methods. These studies indicate that the Lewis acid-base interactions have larger impact on the stability of structure than the C?H…O weak hydrogen bonding. The weak hydrogen bonds in trimethylphosphine oxide have an important role to the large supercritical CO2 solubility when a metal is bound to the oxygen atom of the P=O group. Trimethylphosphate has many Lewis acid-base interaction sites so that it can be dissolved into supercritical CO2 easily even when it has metal ion on the oxygen atom of the P=O group, which is indispensable for a good extractant.

• Bulletin of the Korean Chemical Society
• /
• v.23 no.2
• /
• pp.241-244
• /
• 2002

Using ab initio method, we have studied the structural stabilities, the electronic structures and properties between the two isomers with $C_{2h}$ and $C_{2v}$ symmetry of AlN four-membered-ring single precursors $[Me_2AlNHR]_2$ (R = Me, $^iPr$, and $^iBu$). In the viewpoint of bond lengths in optimized structures, the N-C bonds are considerably affected by the change of the R groups bonded to nitrogen, but the bonding characters of the Al-N and Al-C bonds are little affected. Also the structural stabilities between the two isomers with $C_{2h}$ and $C_{2v}$ symmetry by using Hartree-Fock (HF) and the second order Moeller-Pleset (MP2) calculations agree well with the experimental results for the relative stability of bis(dimethyl- m-isopropylamido-aluminum) (BDPA) and bis(dimethyl- m-t-butylamido-aluminum) (BDBA), while the semiempirical AM1 and PM3 calculations for BDPA were reverse. Thus, our results may aid in designing an optimum precursor for a given process by explaining the experimental results through the elimination of the R groups bonded to nitrogen.

• Bulletin of the Korean Chemical Society
• /
• v.16 no.5
• /
• pp.410-416
• /
• 1995

The proton transfer energies of gas phase glycine and alanine and those of hydrated glycine and alanine were calculated both with Hartree-Fock and $M{\Phi}ller-Plesset$ ab initio molecular orbital (MO) calculations with 6-31G** basis set. The transition states of the proton transfer of gas phase glycine was also investigated. For zwitterions, both for glycine and alanine, the water bound to -NH3+ site stabilize the complex more compared with the water bound to -CO2-. The proton transfer energy, ΔEpt, of glycine, alanine, mono-hydrated glycine, mono-hydrated alanine, di-hydrated glycine and di-hydrated alanine were obtained as 30.78 (MP2: 22.57), 31.43, 23.99 (MP2: 17.00), 24.98, 22.87, and 25.63 kcal/mol, respectively. The activation energy for proton transfer from neutral (Nt) glycine to zwitterion (Zw) glycine, Ea, was obtained as 16.13 kcal/mol and that for reverse process, Ear, was obtained as 0.85 kcal/mol. Since the transition state of the proton transfer of gas phase glycine locate near the glycine zwitterion on the potential energy surface and the shape of the potential well of the zwitterion is shallow, the zwitterion easily changed to neutral glycine through the proton transfer.

• Bulletin of the Korean Chemical Society
• /
• v.30 no.12
• /
• pp.2973-2978
• /
• 2009

The decomposition reaction mechanism of $CH_3CF_2O$ radical formed from hydroflurocarbon, $CH_3CHF_2$ (HFC-152a) in the atmosphere has been investigated using ab-initio quantum mechanical methods. The geometries of the reactant, products and transition states involved in the decomposition pathways have been optimized and characterized at DFT-B3LYP and MP2 levels of theories using 6-311++G(d,p) basis set. Calculations have been carried out to observe the effect of basis sets on the optimized geometries of species involved. Single point energy calculations have been performed at QCISD(T) and CCSD(T) level of theories. Out of the two prominent decomposition channels considered viz., C-C bond scission and F-elimination, C-C bond scission is found to be the dominant path involving a barrier height of 12.3 kcal/mol whereas the F-elimination path involves that of a 28.0 kcal/mol. Using transition-state theory, rate constant for the most dominant decomposition pathway viz., C-C bond scission is calculated at 298 K and found to be 1.3 ${\times}$ 10$^4s{-1}$. Transition states are searched on the potential energy surfaces involving both decomposition channels and each of the transition states are characterized. The existence of transition states on the corresponding potential energy surface are ascertained by performing Intrinsic Reaction Coordinate (IRC) calculation.

• Nuclear Engineering and Technology
• /
• v.53 no.11
• /
• pp.3543-3562
• /
• 2021

The methods and performance of a 3D pin-by-pin neutronics code based on the 2D/1D decoupling method are presented. The code was newly developed as an effort to achieve enhanced accuracy and high calculation performance that are sufficient for the use in practical nuclear design analyses. From the 3D diffusion-based finite difference method (FDM) formulation, decoupled planar formulations are established by treating pre-determined axial leakage as a source term. The decoupled axial problems are formulated with the radial leakage source term. To accelerate the pin-by-pin calculation, the two-level coarse mesh finite difference (CMFD) formulation, which consists of the multigroup node-wise CMFD and the two-group assembly-wise CMFD is implemented. To enhance the accuracy, both the discontinuity factor method and the super-homogenization (SPH) factor method are examined for pin-wise cross-section homogenization. The parallelization is achieved with the OpenMP package. The accuracy and performance of the pin-by-pin calculations are assessed with the VERA and APR1400 benchmark problems. It is demonstrated that pin-by-pin 2D/1D alternating calculations within the two-level 3D CMFD framework yield accurate solutions in about 30 s for the typical commercial core problems, on a parallel platform employing 32 threads.

• Journal of the Korean Chemical Society
• /
• v.61 no.6
• /
• pp.328-338
• /
• 2017

The density functional theory(DFT) and ab initio calculations have been applied to investigate hydrogen interaction of $H_2O_3(H_2O)_n$ clusters(n=1-5). The structures, IR spectra, and H-bonding energies are calculated at various levels of theory. The $trans-H_2O_3$ monomer is predicted to be thermodynamically more stable than cis form at the CCSD(T)/cc-pVTZ level of theory. For clusters, the geometries are optimized at the MP2/cc-pVTZ level of theory. The binding energy of $H_2O_3-H_2O$ cluster is predicted to be -6.39 kcal/mol at the CCSD(T)//MP2/cc-pVTZ level of theory after zero-point vibrational energy (ZPVE) and basis set superposition error (BSSE) correction. This result implies that $H_2O_3$ is a stronger proton donor(acid) than either $H_2O$ or $H_2O_2$. The average binding energies per $H_2O$ are predicted to be 8.25 kcal/mol for n=2, 7.22 kcal/mol for n=3, 8.50 kcal/mol for n=4, and 8.16 kcal/mol for n=5.

• Bulletin of the Korean Chemical Society
• /
• v.34 no.7
• /
• pp.2051-2055
• /
• 2013

Existing theoretical calculations predict that infrared spectra of the two most fundamental reactive carbo-ions, methyl cation $CH{_3}^+$ with $D_{3h}$ symmetry and protonated methyl cation $CH{_5}^+$ with $C_s(I)$, $C_s(II)$, and $C_{2v}$ symmetries, appear together in the 7-${\mu}m$ region corresponding to the C-H bending modes. Vibrational band profiles of $CH{_3}^+$ and $CH{_5}^+$ have been compared by ab initio calculation methods that use the basis sets of MP2/aug-cc-pVTZ and CCSD(T)/cc-pVTZ. Our results indicate that the bands of rotation-vibration transitions of $CH{_3}^+$ and $CH{_5}^+$ should overlap not only in the 3-${\mu}m$ region corresponding to the C-H stretching modes but also in the 7-${\mu}m$ region corresponding to the C-H bending modes. Five band intensities of $CH{_5}^+$ among fifteen vibrational transitions between 6 and 8 ${\mu}m$ region are stronger than those of the ${\nu}_2$ and ${\nu}_4$ bands in $CH{_3}^+$. Ultimate near degeneracy of the two bending vibrations ${\nu}_2$ and ${\nu}_4$ of $CH{_3}^+$along with the stronger intensities of $CH{_5}^+$ in the three hydrogen scrambling structures may cause extreme complications in the analysis of the high-resolution carbo-ion spectra in the 7-${\mu}m$ region.

• Bulletin of the Korean Chemical Society
• /
• v.34 no.10
• /
• pp.3022-3026
• /
• 2013

The adsorption of carbon dioxide on graphene sheets was theoretically investigated using density functional theory (DFT) and MP2 calculations. Geometric parameters and adsorption energies were computed at various levels of theory. The $CO_2$ chemisorption energies on graphene-$C_{40}$ assuming high pressure are predicted to be 71.2-72.1 kcal/mol for the lactone systems depending on various C-O orientations at the UCAM-B3LYP level of theory. Physisorption energies of $CO_2$ on graphene were predicted to be 2.1 and 3.3 kcal/mol, respectively, at the single-point $UMP2/6-31G^{**}$ level of theory for perpendicular and parallel orientations.

• Bulletin of the Korean Chemical Society
• /
• v.33 no.3
• /
• pp.917-924
• /
• 2012

The conformational study of glycerol has been carried out using the M06-2X/cc-pVTZ level of theory in the gas phase and the SMD M06-2X/cc-pVTZ level of theory in water in order to understand its conformational preferences and solvation effects. Most of the preferred conformers of glycerol have two $C_5$ hydrogen bonds in the gas phase, as found by the analysis of calorimetric data. It has been known that the solvation drove the hydrogen bonds of glycerol to be weaker and its potential surface to be fatter and that glycerol exists as an ensemble of many feasible local minima in water. The calculated populations of glycerol in the gas phase and in water are consistent with the observed values, which are better than the previously calculated ones at the G2(MP2), CBS-QB3, and SM5.42 HF/6-31G(d) levels of theory.

• Bulletin of the Korean Chemical Society
• /
• v.33 no.11
• /
• pp.3797-3804
• /
• 2012

We present calculations for valine (Val) - $(H_2O)_n$ (n = 0-5) to examine the effects of microsolvating water on the relative stability of the zwitterionic vs. canonical forms of Val. We calculate the structures, energies and Gibbs free energies of the conformers at B3LYP/6-311++G(d,p), wB97XD/6-311++G(d,p) and MP2/aug-cc-pvdz level of theory. We find that five water molecules are needed to stabilize the zwitterionic form of Val. By calculating the barriers of the canonical ${\leftrightarrow}$ zwitterionic pathways of Val - $(H_2O)_5$ conformers, we suggest that both forms of Val - $(H_2O)_5$ may be observed in low temperature gas phase.