• Proceedings of the Korean Institute of Surface Engineering Conference
• /
• 2016.11a
• /
• pp.97-97
• /
• 2016

Development of advanced Li-ion battery cells with high durability is critical for safe operation, especially in applications to electric vehicles and portable electronic devices. Understanding fundamental mechanism on the formation of a solid-electrolyte interphase (SEI) layer, which plays a substantial role in the electrochemical stability of the Li-ion battery, in a cathode was rarely reported unlike in an anode. Using first-principles density functional theory (DFT) calculations and ab-initio molecular dynamic (AIMD) simulations we demonstrate atomic-level process on the generation of the SEI layer at the interface of a carbonate-based electrolyte and a spinel $LiMn_2O_4$ cathode. To accomplish the object we calculate the energy band alignment between the work function of the cathode and frontier orbitals of the electrolyte. We figure out that a proton abstraction from the carbonate-based electrolyte is a critical step for the initiation of an SEI layer formation. Our results can provide a design concept for stable Li-ion batteries by optimizing electrolytes to form proper SEI layers.

• Bulletin of the Korean Chemical Society
• /
• v.28 no.6
• /
• pp.959-964
• /
• 2007

The benzene complexes with dimethyl sulfur (DMS) and fluorinated DMS (FDMS) have been investigated using ab initio calculations. The natural bond orbital (NBO) charge population on S atom varies remarkably for different conformations of DMS and FDMS, which determines the possible binding modes for their benzene complexes. The electronegative substituent at the methyl group of DMS causes a significant change in the molecular electrostatic potential around the sulfur atom and changes the interaction mode with aromatic ring. It was found that the sulfur…π interaction mode does not occur in the DMS-benzene complex, while it does in the FDMS-benzene complex. Both B3LYP and MP2 methods provide reliable structures, while the interaction energy obtained by B3LYP is unreliable.

• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.294-298
• /
• 2013

A FeRh alloy is a well-known efficient magnetocaloric material and some experimental and theoretical studies of bulk FeRh have been reported already by several groups. In this study we report first-principles calculations on magnetic properties of different thickness FeRh thin films in order to investigate the possibility to enhance further the magnetocaloric efficiency. We used two methods of a Vienna Ab-initio Simulation Package (VASP) code and SIESTA package. We found that the FeRh thin films have quite different magnetic properties from the bulk when the thickness is thinner than 6-atomic-layers. While bulk FeRh has a G-type antiferromagnetic(AFM) state, thin films which are thinner than 6-atomic-layers have an A-type AFM state or a ferromagnetic (FM) state. We will discuss possibility of magnetic phase transitions of the FeRh thin films in the view point of a magnetocaloric effect. And we found 4-, 5-, 6-layers films with Fe surface and 7-layers film with Rh surface are FM and they have relatively small magnetocrystalline anisotropy (MCA) energy about less than 70 meV. The small MCA energy leads to reduction of the strength of magnetic field in operating a magnetic refrigerator.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2005.07a
• /
• pp.172-173
• /
• 2005

The electronic state of ZnO doped with Y was calculated using the density functional theory. In this study, the program used for the calculation on theoretical structures of ZnO and doped ZnO was Vienna Ab-initio Simulation Package (VASP), which is a sort of pseudo potential method. The detail of electronic structure was obtained by the descrite variational $X\alpha$ (DV-$X\alpha$) method, which is a sort of molecular orbital full potential method. The optimized crystal structures obtained by calculations were compared to the measured structure. The density of state and energy levels of dopant elements was shown and discussed in association with optical properties.

• Bulletin of the Korean Chemical Society
• /
• v.14 no.1
• /
• pp.101-107
• /
• 1993

Ab initio calculations are carried out at the 6-311G$^{**}$ level for the $C_{2v}$ interactions of Be and Li atoms with acetylene molecule. The main contribution to the deep minima on the $^3B_2\;BeC_2H_2\;and\;^2B_2 LiC_2H_2$ potential energy curves is the b_2\;(2p(3b_2)-l{\pi}_g^*(4b_2))$ interaction, the $a_1\;(2s(6a_1)-I{\pi}_u(5a_1))$ interaction playing a relatively minor role. The exo deflection of the C-H bonds is basically favored, as in the $b_2$ interaction, due to steric crowding between the metal and H atoms, but the strong in-phase orbital interaction, or mixing, of the $a_1$ symmetry hydrogen orbital with the $5a'_1,\;6a'_1,\;and\;7a'_1$ orbitals can cause a small endo deflection in the repulsive complexes. The Be complex is more stable than the Li complex due to the double occupancy of the 2s orbital in Be. The stability and structure of the $MC_2H_2$ complexes are in general determined by the occupancy of the singly occupied frontier orbitals.

• Bulletin of the Korean Chemical Society
• /
• v.21 no.10
• /
• pp.953-954
• /
• 2000

The equilibrium geometrical structures of silver atom clusters at their electronic ground states have been theo-retically determined by using the nonrelativistic semiempirical INDO/1 method. The clusters investigated are Agn, Agn+, and Agn- (n = 7 , 8, 9). In order to find the most stable structure, i.e., the global minimum in energy hypersurface, geometry optimization and energy calculation processes have been repeatedly performed for all the possible graphical models by changing the bond parameters (resonance integral values). The heptamers are pentagonal bipyramidal-Ag7(D5h), Ag7+ (D5h), Ag7- (D5h); the octamers are pentagonal bipyramidal with one atom capped-Ag8(D2d), Ag8+ (Cs), Ag8- (D2d); the nonamers are pentagonal bipyramidal with two atoms capped -Ag9(C2v), Ag9+ (C2v), Ag9- (C2v). Our structures are in good agreement with those by ab initio calculations ex-cept for the anionic Ag9- cluster. And it is noted that the INDO/1 method can accurately predict the Ag cluster geometries when a proper set of bond parameters is used.

• 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.

• Journal of the Korean Chemical Society
• /
• v.67 no.5
• /
• pp.333-338
• /
• 2023

In this study, the H₂O reaction with SiO clusters was investigated using ab initio Monte Carlo simulations and density functional theory calculations. Three chemistry models, PBE1/DGDZVP (Model 1), PBE1/DGDZVP (Si atom), and aug-cc-pVDZ (O and H atoms), (Model 2) and PBE1/aug-cc-pVDZ (Model 3), were used. The average bond lengths, as well as the relative and reaction energies, were calculated using Models 1, 2, and 3. The average bond lengths of Si-O and O-H are 1.67-1.75 Å and 0.96-0.97 Å, respectively, using Models 1, 2, and 3. The most stable structures were formed by the H transfer from an H₂O molecule except for Si₃O₃-H₂O-1 cluster. The Si₃O₃ cluster with H₂O exhibited the lowest reaction energy. In addition, the Bader charge distributions of the Si_nO_n and (SiO)_n-H₂O clusters with n = 1-7 were calculated using Model 1. We determined that the reaction sites between H₂O and the SiO clusters possessed the highest fraction of electrons.

• Bulletin of the Korean Chemical Society
• /
• v.35 no.5
• /
• pp.1422-1432
• /
• 2014

Convergent all-electron multi-reference configuration-interaction (MRCI) calculations are performed for a lithium dimer with Kaufmann's Rydberg basis functions. A comparison of the results of these calculations with those of the effective core potential/core polarization potential (ECP/CPP) method and experimental data reveals the deficiency of the all-electron ab initio method. The deficiency is related to the mere 51.9% attainment of electron correlation for the ground state. The percent attainment of electron correlation for the first excited state is slightly better than that for the ground state, preventing us from obtaining better agreements with experimental data by means of increasing the size of basis sets. The Kaufmann basis functions are then used with the ECP/CPP method to obtain the accurate convergent potential energy curves for the $^1\prod_u$ states correlated to Li(2p) + Li(2p) and Li(2s) + Li(n = 2, 3, 4). Quantum defect curves (QDCs) calculated for both the $X^2\sum_g$ and 1 $^2\prod_u$ states of the $Li{_2}^+$ ion and the Lu-Fano plot reveal a strong series-series interaction between the two $2snp{\pi}$ and $2pnp{\pi}$ Rydberg series. The QDCs are then used to resolve assignment problems in the literature. The reassignments, performed by Jedrzejewski-Szemek et al., of the dissociation product of the D $^1\prod$ state from (2s+3d) to (2s+3p) and that of the 6 $^1\prod_u$ from (2s+4d) to (2s+4p) are found to be incorrect. It may be more natural to assign their $2snp{\pi}$ Rydberg series as a $2snd{\pi}$ series. The state, assigned as 5p $^1\prod_u$ by Ross et al. and 4d $^1\prod$ by Jedrzejewski-Szemek et al., is assigned as the 7 $^1\prod_u$ state, correlated to the Li(2s) + Li(4f) limit.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.185.2-185.2
• /
• 2014

We have studied the atomic and electronic structure of graphene nanoribbons (GNRs) on a hexagonal boron nitride (h-BN) sheet with intercalated atoms using first-principles calculations. The h-BN sheet is an insulator with the band gap about 6 eV and then it may a good candidate as a supporting dielectric substrate for graphene-based nanodevices. Especially, the h-BN sheet has the similar bond structure as graphene with a slightly longer lattice constant. For the computation, we use the Vienna ab initio simulation package (VASP). The generalized gradient approximation (GGA) in the form of the PBE-type parameterization is employed. The ions are described via the projector augmented wave potentials, and the cutoff energy for the plane-wave basis is set to 400 eV. To include weak van der Waals (vdW) interactions, we adopt the Grimme's DFT-D2 vdW correction based on a semi-empirical GGA-type theory. Our calculations reveal that the localized states appear at the zigzag edge of the GNR on the h-BN sheet due to the flat band of the zigzag edge at the Fermi level and the localized states rapidly decay into the bulk. The open-edged graphene with a large corrugation allows some space between graphene and h-BN sheet. Therefore, atoms or molecules can be intercalated between them. We have considered various types of atoms for intercalation. The atoms are initially placed at the edge of the GNR or inserted in between GNR and h-BN sheet to find the effect of intercalated atoms on the atomic and electronic structure of graphene. We find that the impurity atoms at the edge of GNR are more stable than in between GNR and h-BN sheet for all cases considered. The nickel atom has the lowest energy difference of ~0.2 eV, which means that it is relatively easy to intercalate the Ni atom in this structure. Finally, the magnetic properties of intercalated atoms between GNR and h-BN sheet are investigated.