XVIth International Workshop on
Quantum Systems in
Chemistry and Physics
QSCP logo Abstracts (Posters)
 Main | Registration | Abstract Submission | List of Registrants | List of Abstracts | Program 
Accommodation | Location and Access | Previous Workshops | ISTCP-VII | Photos
Oral presentations only. Posters only. Show all.

First-Principles Calculations of Hydrogen Impurities in Graphenes and Carbon Nanotubes
Mohammad Shafiul Alam,1 Nyayu Siti Nurainun,1,2 Fahdzi Muttaqien,1,2 Agung Setiadi1,2 and Mineo Saito1,3
1Division of Mathematical and Physical Sciences, Graduate School of Natural Science and Technology,Kanazawa University, Kakuma Kanazawa 920-1192, Japan
2Department of Computational Science, Bandung Institute of Technology, Indonesia
3Collaborative Research Center for Frontier Simulation Software for Industrial Science, Institute of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguroku Tokyo 153-8505, Japan
Carbon Nanomaterials have attracted much attention since they are candidates for post-silicon materials. Since hydrogen is a common impurity in graphene [1-3] and carbon nanotubes (CNT), it is important to clarify how the hydrogen impurity affects the electronic structure of graphene and carbon nanotubes.
Now we are going to study mono-hydrogen in graphene, CNT (5,5) and CNT (10,0). We perform spin-polarized generalized gradient approximation by using first-principles calculations. In the most stable structure of mono-hydrogen, the hydrogen atom is bonded to one of the carbon atom in the graphene sheet [Fig. 1(a)] and also hydrogen atom is located on the outer side of CNT (5,5) [Fig. 1(b)] and CNT (10,0) [Fig. 1(c)]. The bond lengths between the hydrogen atom and the nearest neighbor carbon atom are 1.14 Å, 1.12 Å and 1.12 Å for graphene, CNT (5,5) and CNT (10,0), respectively. These bond lengths are close to that in a CH4 molecule (1.09 Å), which is typical for sp3 C-H bond lengths. In these three cases [Fig. 1] the spin polarized states are found to be most stable and the magnetic moments are 1 μB.

[1] M. Khazaei et al. CARBON 47, 3306 (2009).
[2] A. Ranjbar et al. Phys. Rev. B. 82, 165446 (2010).
[3] E. J. Duplock et al. Phys. Rev. Lett. 92, 225502-1 (2004).

Figure 1. Spin densities of the mono-hydrogen in the graphene, CNT (5,5) and CNT (10,0).
First principle theory for the material constants
Marco Anelli, Dan Jonsson, Kenneth Ruud
Centre for Theoretical and Computational Chemistry, Department of
Chemistry, University of Tromsoe , 9037 Tromsoe , Norway
The interactions between matter and electromagnetic fields are commonly analyzed using the constitutive relations,
which relate the applied electric and magnetic fields (E,B) to the response fields (D,H) through the material constants.
Therefore, material constants (also known as constitutive tensors) describe the response of the matter to the external fields.
In the case of a static perturbation, material constants are observable/measurable quantities, and they can been defined using multipole theory.
In contrast, when a dynamic field is considered, multipole theory leads to origin-dependent expressions for the material constants, and thus do not correspond to
physically acceptable observable quantities. There is therefore a need for origin-independent definitions of the material constants.
Origin-independent expressions of the material constants have been derived by Raab and de Lange applying a transformed multipole theory. They have also derived expressions for the material constants based on a covariant formulation of the constitutive relations. Both these theories provide expressions of the material constants
which are covariant, origin-independent and preserve the spatial invariance of the wave equations. In the dynamic case, a unique definitions of the material constants still do not exist, and the possibility of deriving unique definitions for these quantities is still unresolved.
We will discuss our recent analysis of the problem, focusing on the pure electric and magnetic material constants, namely the electric permittivity and the inverse permeability. In particular, we will provide arguments in order to distinguish between the two possible definitions of the material constants proposed so far.
An important part of our analysis will be how to relate the material constants to well-known measurable quantities. We propose that the absorption coefficient and the scattering amplitude, as defined in QED, can be expressed in terms of the electric permittivity, consistently with the definition provided by transformed multipole theory.
The potential energy surface of Li2+(X2Σg+) alkali dimer colliding with the Xe atom
S. Saidi1, C. Ghanmi1, F. Hassen2 and H. Berriche*1,3
1Laboratoire de Physique et Chimie des Interfaces, Département de Physique, Faculté des Sciences de Monastir, Avenue de l’Environnement, 5019 Monastir, Tunisia.
2Laboratoire de Physique des Semiconducteurs et des Composants Electroniques, Faculté des Sciences, Avenue de l’environnement, 5019 Monastir, Tunisie
3Physics Department, College of Science, King Khalid University, P. O. Box 9004, Abha, Saudi Arabia

*Corresponding author: hamid.berriche@fsm.rnu.tn, hamidberriche@yahoo.fr
The potential energy surface describing the collision between Li2+(X2Σg+) alkali dimer and the Xenon atom have been calculated for the equilibrium distance of the Li2+(X2Σg+) and for an extensive range of the two Jacobi coordinates, R and γ. We have used an ab initio approach based on non-empirical pseudopotential and a parametrized l-dependent polarization potential [1]. The core-core interactions for Li+Xe is included using the (CCSD(T)) accurate potential of Lozeille et al [2]. This numerical potential is interpolated using the analytical form of Tang and Toennies [3] for a better description of the interactions at intermediate and large distances between Li+ and Xe. This technique has reduced the number of active electrons of Li2+(X2Σg+)-Xe system to only one active electron. The Three-Body interactions are also considered in this calculation and an analytical fitting of the potential energy surface has been realized. To our knowledge, there are neither experimental nor theoretical studies on the collision between the Li2+(X2Σg+) alkali ionic molecule and the Xenon atom. Therefore, our results are presented here for the first time. In addition the analytical potential energy surface will be used to investigate solvation of Li2+ ionic alkali molecule in Xenon small clusters.

[1] Ph. Durand and J. C. Barthelat, theorit. Chim. Acta. 38 (1979) 283.
[2] Lozeille et al, Phys. Chem. Chem. Phys. 4 (2002) 3601)
[3] K. T. Tang and J. P. Toennies, J. Chem. Phys. 80 (1984) 3726.
Solvation of Li2+(X2Σg+) in small Xen clusters: structure and stability
S. Saidia, C. Ghanmia, F. Hassenb and H. Berriche*a,c
aLaboratoire de Physique et Chimie des Interfaces, Faculté des sciences, Université de de Monastir, Avenue de l’Environnement, 5019 Monastir, Tunisie
bLaboratoire de Physique des Semiconducteurs et des Composants Electroniques, Faculté des Sciences, Avenue de l’environnement, 5019 Monastir, Tunisie
cPhysics Department, College of Science, King Khalid University, P. O. B. 9004, Abha, Saudi Arabia

*E-mail address: hamidberriche@yahoo.fr
The structure and geometries of the small Li2+(X2Σg+)Xen clusters taken in different geometries with special symmetry groups, are examined. A one-electron pseudopotential approach is used to generate the potential energy surface of the Li2+-Xe triatomic system for a fixed distance between the two Li alkali atoms (the Li2+ equilibrium distance) and for different angles and distances between the Xenon and the Li2+ center of mass. In this approach, the effect of the Li+ core and the electron-Xe interactions are replaced by effective potentials. This has permitted to reduce the number of active electrons of the triatomic system to only one valence electron. In addition, the Li2+Xe interactions are introduced using an analytical form to describe the collision between the alkali ionic dimer and the Xe atom. However, a Lennard-Jones form is used for the Xe-Xe interactions and then included in the total potential energy surface describing the Li2+-Xen cluster. The optimized and stable geometries are performed using the Monte-Carlo Bassin Hoping method. This study has shown that the optimal structures correspond to those where Xenon atoms are aggregated on one or both ends of Li2+ dimer.

Key words: Pseudopotential, Microsolvation of clusters, Monte-Carlo, Structure,
Examining the Limits of Physical Theory:
Analytical Principles and Logical Implications
Erkki J. Brändas
Department of Physical and Analytical Chemistry
Quantum Chemistr
Uppsala University
SE-751 20 Uppsala, Sweden
Our aim is to examine the characteristics and the rationale for developing an analytic foundation for rigorous extensions of quantum mechanics beyond its long-established domain in physics, chemistry and biology. Paradoxical and inconsistent issues related to the various attempts to apply microscopic organization to derive scientific laws in the macro-world are considered. The theoretical framework is augmented with quantum logical principles via a reformulation of Gödel’s theorems. We arrange the assemblage of the mathematical ideas as follows. First we give a detailed examination of the second order differential equation with respect to specific boundary conditions and associated spectral expansions, followed by a description in terms of general complex symmetric representations exemplified and derived from dilation analytic transformations. Associated dynamical time scales are represented and investigated via the corresponding Dunford formula. Relevant applications are reviewed and compared with conventional scattering theory analysing the directive performance and stability of the method. The manifestation and generation of triangular Jordan block entities are derived and further investigated in thermally excited scattering environments of open dissipative systems. Illustrative applications to condensed- and soft condensed matter are provided, and a surprising treatment is given to the Einstein laws of relativity. As a conclusion we emphasize the computational and model building advantages of a conceptual extension of quantum mechanics to rigorously incorporate universal complex resonance structures, their life times and associated localization properties. We also prove the appearance of non-conventional time evolution including the emergence of Jordan blocks in the propagator, which leads to the origin of so-called coherent dissipative structures derived via uniquely defined spatio-temporal neumatic units. This organization yields specific information bearing transformations, cf. the Gödel encoding system, which might connect developmental and building matters with functional issues within a biological framework at the same time providing back-ground dependent features of both special and general relativity theory.
Ab initio calculations of methane dimer interaction energies and
molecular dynamics simulation of fluid methane
Arvin Huang-Te Li and Sheng D. Chao
Institute of Applied Mechanics, National Taiwan University, Taipei 106, Taiwan ROC.

Intermolecular interaction potentials of the methane dimers have been calculated for 12 symmetric conformations using the second-order Møller-Plesset (MP2) perturbation theory and the coupled-cluster with single and double and perturbative triple excitations (CCSD(T)) theory. With increasing basis size, a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy (~0.01 kcal/mol). Only the BSSE corrected results systematically converge to the destined potential curve with increasing basis size. The binding energy calculated and the equilibrium bond length using the CCSD(T) method are close to the results at the basis set limit. For molecular dynamics (MD) simulation, a 4-site potential model with sites located at the hydrogen atoms was used to fit the ab initio potential data. MD simulations using the ab initio PES show good agreement on both the atom-wise radial distribution functions and the self-diffusion coefficients over a wide range of experimental conditions.
Ultrafast electronic motion in hydrogen molecular ion induced by a high power intense laser
H. Mineo1 , Y. Teranishi2 , S.D. Chao1, and S.H. Lin3,4
1 Institute of Applied Mechanics, National Taiwan University, Taipei 106, Taiwan ROC
2 Department of Physics, National Chiao Tung University, Hsinchu 300, Taiwan
3 Institute of Atomic and Molecular Sciences, Academia Sinica, P. O. Box 23-166, Taipei 106, Taiwan ROC
4 Institute of Applied Chemistry, Institute of Molecular Science, Chiao-Tung University, Hsin-Chu, Taiwan ROC
In this report we show how to control the electronic localization in a molecular ion on attosecond time scale by high intense laser based on two different excitation mechanisms, i.e., one takes place during ionization, and the other one sequentially takes place after the ionization [1]. The electronic excited states of hydrogen molecular ion are created during the ionization by taking the configuration interaction (CI) mixing of neutral molecule into account. Photoionization process is calculated by the generalized Keldysh-Faisal-Reiss (KFR) theory [2]. We detect ultrafast oscillatory electronic motion between two atoms in a hydrogen molecular ion due to the creation of excited states during the ionization.

[1] H. Mineo, Y. Teranishi, S.D. Chao, and S.H. Lin, CPL 499, 45 (2010).
[2] L. V. Keldysh, Sov. Phys. JETP 20, 1307 (1965); F. H. M. Faisal, J. Phys. B6, L89 (1973); H. R. Reiss, Phys. Rev. A22, 1786 (1980).
Valence XPS and Raman Spectral analysis of chitosan film modified by Kr+ ion beam bombardments by quantum chemical calculations
K. Endo1, H. Shinomiya1, T. Ida2, S. Shimada2,K. Takahashi1, Y. Suzuki3, H. Yajima1
1Center for Colloid and Interface Science, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601 Japan,
2Laboratory of Theoretical Chemistry, Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa 920-1192, Japan,
3Advanced Development and Supporting Center, RIKEN, 2-1 Hirosawa, Wako, Saitama, 351-0198, Japan
Valence X-ray photoelectron and Raman spectra of chitosan film modified by Kr+ ion beam bombardments were analyzed from quantum chemical calculations. Experimental Raman spectra of the Kr+ ion bombarded film were found-out to be due to four component contributions of chitosan (Chito), diamond-like carbon (DLC), graphite carbon (GC) and amorphous carbon (AC). By considering the four components contribution, we performed depth profile assignments of nm and μm orders for the chitosan film in valence XPS and Raman experiments, respectively from MO-DFT hybrid calculations in GAUSSIAN 09 using the model molecules of the four components. Carbonizations of the film by Kr+ irradiation were obtained as Chito: DLC: AC : GC = 2:1:0.5:0.375 in the μm order from Raman shift spectral analysis, while they were determined as Chito: DLC: AC : GC = 2:1:1:2 in the nm order from valence XPS analysis.

[1] K. Takahashi, R. Shizume, K. Uchida, H. Yajima, J. Biorheol, 21,64(2009).
[2] S. Danielache, M. Mizuno, S. Shimada, K. Endo, T. Ida, K. Takaoka, E.Z. Kurmaev, Polym. J. 37,21-29(2005).
[3] K. Tamura, K. Endo, Y. Takagi, K. Kato, D. Matsumoto, T. Ida, M. Mizuno, Y. Suzuki, K. Takahashi, K. Uchida, H. Yajima, J. Surf. Anal. 14, 344-347(2008).
[4] Gaussian 09, Revision B.01, http://www.gaussian.com.
Valence XPS, IR,and C13 NMR Spectral analysisof 6 polymers by quantum chemical calculations
Kazunaka Endo1, Tomonori Ida2, Shingo Shimada2
1Research Center for Colloid and Interface Science, Tokyo University of Science 1-3, Kagurazaka, Shinjuku-ku 162-8601, Japan,
2Laboratory of Theoretical Chemistry, Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa 920-1192, Japan
Valence XPS, IR, and C13 NMR spectra of 6 polymers (PE, PS, PMMA, PET, Nylon6, PVC) have been analyzed by quantum chemical calculations using the polymer model oligomers from B3LYP/6-31G** basis calculations in GAUSSIAN 09. Simulated valence XPS, IR, and C13 NMR spectra of PS as an example are in good accordance with the experimental results.

[1] Gaussian 09, Revision B.01, http://www.gaussian.com.
[2] K. Endo, Y. Kaneda, H. Okada, D. P. Chong, P. Duffy, J. Phys. Chem.,100, 19455 (1996).
[3] Scott A.P., Radom L., J. Phys. Chem., 100, 16502 (1996).
Second Hyperpolarizabilities of Open-Shell Singlet Extended Metal Atom Chains (EMACs)
Hitoshi Fukui,1 Yudai Inoue,1 Yasuteru Shigeta,1 Benoît Champagne2 and Masayoshi Nakano1
1Graduate School of Engineering Science, Osaka University, Japan
2Laboratoire de Chimie Théorique, Facultés Universitaires Notre-Dame de la Paix (FUNDP), Belgium
We have theoretically investigated open-shell singlet molecules as a novel class of nonlinear optical (NLO) systems, and have revealed their structure−property relationship: singlet diradical systems with intermediate diradical characters (y) tend to exhibit larger second hyperpolarizabilities (γ) than pure diradical and closed-shell systems [1,2]. In addition to organic compounds, polynuclear transition metal complexes are expected to show singlet multiradical nature due to the unpaired electrons in the metal d-orbitals [3]. In particular, open-shell singlet extended metal atom chains (S-EMACs), which have metal−metal direct bonds between transition metal atoms are promising NLO candidates because of the extended d−d conjugation and the singlet multiradical characters originating from the interactions between the d-orbitals, which can create dσ, dπ and dδ electron conjugation. In the present study, we theoretically investigate the multiradical characters of the dσ, dπ and dδ orbitals as well as the γ values of S-EMACs of chromium(II) without ligands. We also analyze the contributions of these orbital electrons to γ [γ(dX), where X = σ, π or δ] in order to reveal relationships between the orbital symmetry, the multiradical character and the γ value.

[1] M. Nakano et al., J. Phys. Chem. A 109, 885 (2005).
[2] M. Nakano et al., J. Chem. Phys., 125, 074113 (2006); Phys. Rev. Lett., 99, 033011 (2007).
[3] H. Fukui et al., J. Phys. Chem. Lett. 10.1021/jz2007897 (2011).
Laser Electron-Gamma-Nuclear Spectroscopy of Atoms and Multicharged Ions and NEET Effects in Heavy Nuclei: Relativistic Energy Approach
Alexander V. Glushkov,1,2 Olga Yu. Khetselius1 Svetlana Malinovskaya1 and Andrey A. Svinarenko1
1Odessa University, Odessa-9, SE, Ukraine

2ISAN, Russian Academy of Sciences, Troitsk, Moscow reg., Russia
In the resonant process of nuclear excitation by electron transition (NEET) or electron capture (NEEC) an electron is captured into a bound atomic shell with the simultaneous excitation of the nucleus. The excited nucleus can then decay radiatively or by internal conversion. In the latter case, a resonant inelastic electron scattering on the nucleus occurs. Here we present consistent, relativistic approach to calculation of the probabilities of the different cooperative laser electron-gamma-nuclear processes in atoms, ions, nuclei and resonant NEET (NEEC) processes in heavy nuclei, based on the relativistic density functional (DF) formalism and energy approach (S-matrix formalism of Gell-Mann and Low) [2]. Decay and excitation probability is linked with the imaginary part of energy of the excited state for the “electron shell- nucleus-photons” system. For radiate decays it is manifested as effect of retarding in interaction and self-action and calculated within QED- DFT theory [2]. We firstly present data about intensities of the electron satellites in gamma-spectra of nuclei in the neutral (low lying transitions) and multicharged O-and F-like ions for isotopes 57Fe, 133Cs, 171Yb and discover a new effect of the giant increasing electron satellites intensities under transition from the neutral atoms to multicharged ions. We present new, more accurate data about NEET probabilities in the nuclei of 189Os, 197Au (with comparison with theoretical data by Tkalya and experimantal data of Argonne Nat.Lab. and Japan Synchrotron Centre [3]) and firstly for nuclei of 193Ir, 235U, 268Mt.
[1]. L.N.Ivanov, V.S.Letokhov, JETP. 93, 396 (1987); A.V.Glushkov, L.Ivanov, Phys.Lett.A 170, 33 (1992); A.V.Glushkov, L.N.Ivanov, V.S.Letokhov, Preprint of ISAN N AS-4, Troitsk, (1992); E.V.Tkalya, Phys.Rev.A.75, 022509 (2007); T.J. Burvenichm J.Evers, C.H.Keitel, Phys. Rev. C. 74, 044601 (2007) ; A.Shahbaz, C.Muller, A.Staudt, T.J.Burvenich, C.H.Keitel, Phys.Rev.Lett.98, 263901 (2007).
[2]. A.Glushkov et al, J.Phys.CS. 11, 188 (2005); 35, 425 (2005); Int.J.Quant.Chem. 104, 512, 562 (2005); 99, 889, 936 (2004); Europ.Phys.Journ. 160, 195 (2009); Phys.Scr.T135, 014022 (2009).
[3]. S.Kishimoto, Y.Yoda, Y.Kobayashi etal, Phys.Rev.C74, 031301 (2006); I. Ahmad, R.Dunfird, H.Esbensen etal, Phys.Rev.C61, 051304 (2000).
Theoretical Study on the Redox Reaction of Azurin in Water Solvent
Masashi Iwayama, Hiroaki Saito, Kazutomo Kawaguchi and Hidemi Nagao
Division of Mathematical and Physical Science, Graduate School of Natural Science and Technology, Kanazawa University, Japan
Metalloproteins including a transition-metal ion in the active site have important characteristics such as the metabolism and the intracellular signaling in vivo. In metalloproteins, Azurin is a kind of blue copper protein containing a copper atom in the active site, and has been known to show the oxidation-reduction reaction by charge transfer between the active sites of proteins in the solution. In this study, we carry out molecular dynamics simulations and the quantum mechanics calculations of the oxidized Azurin (Az(Ⅱ)) and reduced Azurin (Az(Ⅰ)) in water solvent to reveal the redox character of Azurin in water solvent. According to the Born-Haber Cycle model, the free energy change of the solute molecule in solvent can be expressed as a summation of the free energy change of solute in gas phase and the difference of solvation free energies of the solute. From these results, we estimate of the redox potential of Azurin and compare with the experimental value.
Calculation of Magnetic Properties and Spectroscopic Parameters of Manganese Clusters with Density Functional Theory Methods
Keita Kanda, Toru Saito, Yasutaka Kitagawa, Takashi Kawakamai, Shushuke Yamanaka, Kizashi Yamaguchi and Mitsutaka Okumura
Graduate School of Science, Osaka University, Japan
Many manganese complexes bridged by μ-oxo ligands are synthesized as model of the oxygen-evolving complex (OEC) in Photosystem II. Catalytic cycle and oxidation states of Mn cluster in OEC is still unrevealed. Spectroscopic parameters for those model complex, especially hyperfine coupling constants (HFCs) measured by EPR suggest deep insight for electronic structure.
Many calculation of HFCs have been carried out with ab initio DFT methods, however, there is very few investigation of functional dependence. Today’s Density Functional Theory approximates exchange-correlation terms, which may cause inaccurate spin density nearby nuclei. We have done benchmark calculations to check the accuracy of various functionals for HFC of several manganese-μ-oxo systems. In addition, the most reliable functional for estimations of s-d exchange interactions was identified.
Ab initio and density functional calculation of calcium binding sphingomyelin lipid molecules: A pin holder model approach
Hiroyuki Kawabe1 and Kimikazu Sugimori2
1Department of Social Work, Faculty of Social Work, Kinjo University, Japan

2Department of Physical Therapy, Faculty of Health Sciences, Kinjo University, Japan
One of the phospholipids, sphingomyelin (SM, N-acyl-sphingosine-1-phosphorylcholine) is the most abundant component of mammalian membranes in brain and nervous tissues. It plays an important role for apoptosis, aging, signal transduction with cations. Recently, Yappert and co-workers have shown that human lens sphingomyelin and its hydrogenated derivative, dihydrosphingomyelin (DHSM) are interacted with Ca2+ ions to develop human cataracts [1,2]. Previously we have investigated conformational differences between an isolated SM / DHSM molecule and Ca2+ -coordinated form by using density functional theory (DFT) which B3LYP functional and 6-31G(d,p) double-zeta split-valence basis set is applied for geometry optimization and normal mode analysis. As a result, one of conformers of SMs has hydrogen bonding between hydroxyl group and phosphate group, whereas another conformer has hydrogen bonding between hydroxyl and phosphate amide group [3]. Moreover, 31P-Nuclear Magnetic Resonance (NMR) shielding constants of the obtained conformers are investigated by using ab initio and DFT with gauge invariant atomic orbitals (NMR-GIAO) calculations [4].
In this study, we apply previous results to dimer and trimer model of lipid molecules with Ca2+ and water molecules by quantum chemical calculations. We assume that tails of the lipids arbitrary fixed specific length in plasma membrane, like pin holder. Except for fixed tail groups, all of molecular structure of head group and solvated Ca2+ are fully optimized. We analyze potential energy surface (PES) with the fixed parameters that are length between tails of lipid and coordinate of Ca2+ approaching from outer region to phosphate.

[1] M. Rujoi, D. Borchman, D. B. Dupré, and M. C. Yappert, Biophys. J. 82 (2002) 3096-3104.
[2] M. C. Yappert and D. Borchman, Chem. Phys. Lipids 129 (2004) 1-20.
[3] K. Sugimori, H. Kawabe, H. Nagao and K. Nishikawa, Int. J. Quantum Chem. 108, 2935-2942 (2008). doi:10.1002/qua.21858
[4] K. Sugimori, H. Kawabe, H. Nagao, and K. Nishikawa, Int. J. Quantum Chem., 109 (2009) 3685-3693. doi:10.1002/qua.22404
Structural analysis of ligand binding mechanism of Hsp90
Kazutomo Kawaguchi, Acep Purqon, Hiroaki Saito, Hidemi Nagao
Institute of Science and Engineering, Kanazawa University, Kanazawa 920-1192 Japan
Hsp90 is one of the molecular chaperones which protect their client proteins in the folding processes. Hsp90 usually forms a dimer and undergoes large conformational changes during its functional cycle. Binding and hydrolysis of ATP to ADP are required for these conformational changes. We performed molecular dynamics simulation of Hsp90 with ADP to elucidate molecular mechanisms of ligand binding process. The results of structural analysis of ligand binding process of Hsp90 will be presented and discussed through recent calculations.
Blackbody, synchrotron radiation, bremsstrahlung and plasmon analyzed by Tsallis nonextensive entropy
Jun Kawai, Abbas Alshehabi, Hiroyuki Iwasaki, Koretaka Yuge, Ágnes Nagy*
Kyoto University, Department of Materials Science and Engineering,
Sakyo-ku, Kyoto 606-8501, Japan
*Debrecen University, H-4010 Debrecen, Hungary
Manne Siegbahn once wrote that “the curve (bremsstrahlung from X-ray tube) rises abruptly on the short wave-length side. In this respect the curve differs fundamentally from the curve for the radiation of a ‘black’ body, though otherwise the curves have an external resemblance”[1]. Blackbody and bremsstrahlung or synchrotron radiation are similar in their overall spectral shape but a minor difference exists [2-4]. The Planck’s blackbody equation is expressed through the Tsallis nonextensive entropy [5]: N(E)=1/[(1+(q-1)x)1/(q-1)-1], where N is the number of photons, x=E/kT, and q the Tsallis parameter [6,7]. When q->1 , the above equation becomes Plank’s blackbody equation. We have fitted Eqn.(1) to a synchrotron radiation (relativistic bremsstrahlung) and when q=1.05 the agreement is best and satisfactory as shown in Fig. 1. The parameter q >1 means the long-range correlation due to nonextensitivity. The blackbody is ideal gas of photons without interaction, but the photons are weakly interacting each other for the synchrotron radiation.
Similarly the plasmon energy loss peaks when Si is irradiated by e.g., 1500 eV electron beam, decays exponentially from the 1st loss to the higher orders. The deviation from the exponential decay is also expressed by the Tsallis q parameter, which indicate the intrinsic and extrinsic paths of the plasmon energy loss processes.
exp(-E/kT)appears in physical chemistry, and the deviation from exponential could be interpreted by the long range interaction through the Tsallis parameter q.

[1] M. Siegbahn, “The Spectroscopy of X-Rays”, transl. G. A. Lindsay, Oxford Univ. Press, London, 1925, p. 206.
[2] J. Kawai, H. Ishii, Spectrochim. Acta, Part B, 60, 1586 (2005).
[3] J. Kawai, H. Ishii, Radiation Phys. Chem., 75, 1716 (2006).
[4] T. Tanigaki, J. Kawai, X-Ray Spectrom., 36, 321 (2007).
[5] C. Tsallis, Introduction to Nonextensive Statistical Mechanics, Springer (2009).
[6] Q. A. Wang, A. Le Mehaute, Phys. Lett. A, 242, 301 (1998).
[7] S. Martnez, F. Pennini, A. Plastino, C. J. Tessone, Physica A, 309, 85 (2002).

Figure 1. Comparison of synchrotron radiation (dotted)and blackbody radiation (solid line) with q=1.05.

Figure 2. Plasmon energy loss spectra.

Calculation of magnetic constants D in Zero-Field Splitting by ab initio methods
Takashi Kawakami, Keiji Kinoshita, Akira Ito, Yasutaka Kitagawa, Shusuke Yamanaka, Kizashi Yamaguchi, Mitsutaka Okumura
Graduate School of Science, Osaka University, Japan
For theoretical treatments of ZFS parameters, previously Pederson and co-workers have developed a method calculating the ZFS parameters with the DFT methods. In the series of our studies, we developed new ab initio MO program package ("Q" by Ryo Takeda in our group). The Pederson's treatments are also included in this program. Thus, several ab initio MO methods can be applied to calculating zero-field splitting parameters caused by mainly spin-orbit coupling. We examine the behavior and tendency of the method by applying it to some basic molecules such as small molecules (ex. carbene), pure organic molecular magnets, single-molecular magnets (SMM), etc. It shows the method has good accuracy. We suggest applicable region of the method.

In addition, Neese has also developed another calculation scheme with coupled-perturbed equation. This method is involved in their program package "ORCA". Thus, we also evaluated ZFS parameters by Neese's methods.
Free energy of cell-penetrating peptide in lipid bilayer membrane: coarse-grained simulation
Shuhei Kawamoto1, Takeshi Miyakawa2, Masako Takasu2, Ryota Morikawa2, Tatsuki Oda1, Hiroaki Saito1, Shiroh Futaki3 and Hidemi Nagao1
1Graduate School of Natural Science and Technology, Kanazawa University
2School of Life Sciences, Tokyo University of Pharmacy and Life Sciences
3Institute for Chemical Research, Kyoto University
Cell-penetrating peptides (CPPs) can permeate plasma membrane [1], and are useful for delivery of drugs and transfection of DNA. Free energy ΔG to insert a molecule into lipid bilayer membrane is important to investigate the permeability of the molecule through the membrane. ΔG of small molecules, such as water and oxygen, are obtained by all-atom simulations [2] with cavity insertion method. However, large molecules, such as CPP, are difficult to insert into the cavity in lipid bilayer membrane.
We estimate ΔG of CPP by steered MD method with our coarse-grained model [3]. When there is large attractive potential between CPP and lipid heads, the CPP makes inverted micelle in the lipid bilayer membrane, and the difference of ΔG between water region and lipid bilayer region is more than 66 kBT. It suggests that the spontaneous permeation of CPPs needs large driving force to overcome the barrier.

[1] I. Nakase, T. Takeuchi, G. Tanaka, S. Futaki, Adv. Drug Deliv., 60, 598 (2008)
[2] K. Shinoda, W. Shinoda, T. Baba, M. Mikami, J. Chem. Phys., 121, 9648 (2004)
[3] S. Kawamoto, et al., J. Chem. Phys., 134 095103 (2011)
Nuclear-Relativistic Many-Body Perturbation Theory to Parity Nonconservation Effect in Heavy Atoms and Nuclei
Olga Yu. Khetselius
Department of Mathematics, Odessa OSENU University, Odessa-9, Ukraine
During the past two decades, the nuclear and atomic-optical experiments to detect parity non-conservation (PNC) have progressed to the point where PNC amplitudes can be measured with accuracy on the level of a few percents in certain heavy atoms and significantly worse in some nuclei (Mossbauer spectroscopy). Nowadays the PNC in atoms has a potential to probe a new physics beyond the standard model. Promising idea (Forston) is to apply the techniques of laser cooling and ion trapping to measurement of the PNC in 6s2S1/2-5d2D3/2 transition of the Ba+. In our paper we systematically apply the nuclear-relativistic many-body perturbation theory formalism [1] to precise studying PNC effect in heavy atoms with account for nuclear, correlation and QED corrections. There are determined the PNC radiative amplitudes for a set of nuclei (atoms): 133Cs, 137Ba+, 173Yb with account of the exchange-correlation, Breit, weak e-e interactions, QED and nuclear (magnetic moment distribution, finite size, neutron skin) corrections, nuclear-spin dependent corrections due to anapole moment, Z-boson [(AnVe) current] exchange, hyperfine-Z exchange [(VnAe) current]. The weak charge is found for 133Cs, 205Tl and 173Yb and comparison with Standard Model is done. Using the experimental value (EPNC/b)= 39mV/cm (Tsigutkin et al, 2009) and our calculated amplitude value 9.707*10-10 ieaBone could find for 173Yb (Z=70, N=103) the weak charge value QW=-92.31 (the SM gives QW=-95.44). The received data are compared with known earlier and recent results [1,2]. The role of the nuclear effects contribution (core-polarization contributions, which are induced by valent protons of a nucleus), spatial distribution of magnetization in a nucleus (the Bohr-Weisskopf effect), neutron skin correction and the non-accounted high order QED corrections are analyzed.
[1] O.Khetselius, Phys.Scripta T34, 014023 (2009); Int.J.Quant.Chem. 109, 3330 (2009); A.V.Glushkov, O.Khetselius etal, Frontiers in Quantum Systems in Chem. and Phys. (Springer) 18, 505 (2008).
[2] W. Johnson, M.S.Safronova, U.I.Safronova, Phys.Rev. A69, 062106 (2003); V.V.Flambaum, J.S.Ginges, Phys.Rev. A72, 052115 (2005).
Computational study of conformational preferences in intermediates and transition states of the hydrolysis of dimethyl phosphate
Makoto Kita,1 Haruki Nakamura1 and Yu Takano1
1Institute for Protein Research, Osaka University, Japan
Phosphate diester is a basic structure in DNA and RNA. Hydrolysis mechanism of phosphate diester is important for understanding decomposition reaction of nucleic acid.
In this study, we have explored the hydrolysis reaction pathway of dimethyl phosphate, which is the simplest phosphate diester, with two nucleophiles (hydroxide and water), using the density functional theory (B3LYP/6-311++G(2d,2p)) with PCM (polarizable-continuum model). Since it was reported that the conformations of intermediates and transition states influence the reaction mechanism of transacylation of methyl acetate with methoxide [1], we focused on the conformational preferences in the hydrolysis reaction. The computed reaction pathway was confirmed by IRC (intrinsic reaction coordinate) calculation.
In hydrolysis of dimethyl phosphate with hydroxide, a concerted reaction pathway was obtained in the gas phase, whereas a stepwise reaction pathway was obtained in the aqueous phase. In the gas phase, the conformation of the transition state was different from that of the earlier study [2], but the reaction mechanism was similar to each other. The conformations hardly influence the energetics of the gas-phase hydrolysis of dimethyl phosphate. On the other hand, in the aqueous phase, the reaction pathway and the activation and reaction energies were very different from those of the earlier study [2], due to the differences in the computational procedure.
We have also explored the reaction mechanism of the hydrolysis of dimethyl phosphate with water. In the gas phase, the most stable conformation of the reactant (conf1) changed to the other conformation (conf2). The reaction proceeded with conf2, and the stepwise reaction pathway was obtained. This is because the activation energy of the hydrolysis of conf2 is lower than that of conf1, and because the activation energy of conformational change from conf1 to conf2 is as low as 2.2 kcal/mol. In the aqueous phase, the proton transfer, the first step of hydrolysis, did not occur, implying that the hydrolysis occurs by a proton shuttle mechanism.

[1] Takano, Y.; Houk, N. K. J. Phys. Chem. A 2004, 108(52), 11740-11751.
[2] Ribeiro, J. M. A.; Ramos, J. M.; Fernandes, A. P. J. Chem. Theory Comput. 2010, 6(8), 2281-2292.
The band 12 issue of norbornane: a comparison between symmetry adapted cluster expansion configuration interaction (SAC-CI) and the third order algebraic diagrammatic construction scheme [ADC(3)]
S. Knippenberg1 and B. Hajgató
1 Service de Chimie des Matériaux Nouveaux, Université de Mons, Place du Parc 20, B-7000 Mons, Belgium

2 Eenheid Algemene Chemie, Vrije Universiteit Brussel (VUB), Faculteit Wetenschappen, Pleinlaan 2, 1050 Brussels, Belgium
In line with a recent study of the electronic structure of the cage compound norbornane [1, 2], symmetry adapted cluster expansion configuration interaction (SAC-CI) general R calculations have been performed and compared with results obtained by the third order algebraic diagrammatic construction scheme [ADC(3)]. Comparison has been made with previously performed electron momentum spectroscopy (EMS) and ultra violet photo-electron measurements. The region around 25 eV (band 12), characterized by an elaborated band in the EMS spectrum which is missing in previous Green's function and ADC calculations, is investigated. This study is completed with outer-valence Green's function (OVGF) and SAC-CI/SD-R calculations. Since ADC(3) only includes 2h-1p shake-up states, while SAC-CI general-R also includes higher order ones, the agreement between both methods assures that the higher order shake-up states do not play an important role in the ionization spectra of norbornane. The band-12 issue of norbornane is therefore still open, while a tentative description in terms of ultrafast nuclear dynamical effects and autoionization processes has become more plausible.

[1] S. Knippenberg, K. L. Nixon, M. J. Brunger, T. Maddern, L. Campbell, N. Trout, F. Wang, W. R. Newell, M. S. Deleuze, J.-P. François, D. A. Winkler, J. Chem. Phys. 121 (2004), 10525.
[2] S. Knippenberg, M. S. Deleuze, T. J. Cleij, J.-P. François, L. S. Cederbaum, J. H. D. Eland, J. Phys. Chem. A 109 (2005), 4267.
Analysis for constructing protein nano-fiber including metal ions
Yu Komatsu,1 Shuhei Kawamoto,2 Takeshi Miyakawa,1 Ryota Morikawa,1 Masako Takasu,1 Satoshi Akanuma1 and Akihiko Yamagishi1
1School of Life Sciences, Tokyo University of Pharmacy and Life Sciences, Japan
2Graduate School of Natural Science & Technology, Kanazawa University, Japan
As materials used in nanotechnology, proteins can be applied to many fields. Our aim is to fabricate a nano-fiber using arbitrary proteins as maintaining their native structures. By inducing mutations in the adaptor proteins, hydrophobic and electrostatic interaction can be modified. By means of mixing the two proteins, the fiber may be organized spontaneously. As the materials of nano-fiber formation, we are testing three proteins: LARFH, Sulerythrin and IPMDH.
In this research, we analyzed the three proteins which are candidates for constructing a nano-fiber and investigated conditions that show high stability using molecular dynamics (MD). Sulerythrin contains two pairs of Zn2+ and Fe2+. We perform quantum calculation to investigate how these ions affect the stability as a fiber by quantum computing. Moreover, we perform umbrella sampling simulations to obtain PMF (Potential of Mean Force), from which binding energy Δ Gbind is calculated.

[1]S. Akanuma et al., J. Biochem, 147, 371 (2010)
[2]J. A. Lemkul et al., J. Phys. Chem. B, 114, 1652 (2010)
General Coalescence Conditions for the Exact Wave Functions
Yusaku I. Kurokawa, Hiroyuki Nakashima, and Hiroshi Nakatsuji
Quantum Chemistry Research Institute, JST, CREST, Kyodai Katsura Venture Plaza 106, Goryo Oohara 1-36, Nishikyo-ku, Kyoto 615-8245, Japan
  We derived necessary conditions for exact wave functions of the time-independent Schrödinger equation that must be satisfied at a coalescence (or cusp) region. Some of such conditions are already known as the Kato's cusp conditions (CC) [1] and Rassolov and Chipman's CC [2]. In this study, we have generalized them to be relations among higher order derivatives of wave functions. Furthermore, we have extended them to be applicable not only to the Coulombic system but also to any systems where the interaction between two particles is represented in power series of inter-particle distance. We named these conditions to be General Coalescence Conditions (GCCs). Any wave functions must satisfy GCCs; otherwise they never become exact. The Kato's CC and Rassolov and Chipman's CC are included in GCCs as special cases.
  We applied the Free Complement (FC) wave functions [3] (which is nearly exact) of the hydrogen atom in the ground and excited states, the harmonic oscillator, a system with interacting potential of V = r, and the helium atom to GCCs and confirmed that the FC wave functions satisfy GCCs.

 [1] T. Kato, Commun. Pure Appl. Math. 10, 151 (1957).
 [2] V. A. Rassolov and D. M. Chipman, J. Chem. Phys. 104 (24), 9908 (1996).
 [3] H. Nakatsuji, Phys. Rev. Lett. 93 (3), 030403 (2004); H. Nakatsuji, Phys. Rev. A 72 (6), 062110 (2005).
Laser Electron-Gamma-Nuclear Spectroscopy of Diatomic and Multiatomic Molecules
Svetlana V. Malinovskaya,1 and Andrey A. Svinarenko2
1Department of Chemistry, Odessa University, Odessa-9, Ukraine

2Department of Mathematics, Odessa University, Odessa-9, Ukraine
An important class of problems connected with modelling the cooperative laser-electron-gamma-nuclear phenomena in diatomic and multiatomic molecules [1,2] is now of a great interest. It includes a calculation of the probabilities and energies of the mixed gamma-nuclear and optical quantum transitions in molecules, intensities of the complicated gamma-transitions due to the changing of the molecular excited states because of the gamma nuclear transition. We present a consistent, quantum approach to calculation of the probabilities of the different cooperative laser electron-gamma-nuclear processes in different molecules [2] (including the set of electron or vibration-rotational satellites of the gamma-nuclear spectrum). The calculation results for the electron-gamma-nuclear transition probabilities in the diatomic (the nucleus 127I with E=203keV in molecule of H127I) and vibration-nuclear transition probabilities for some three-atomic XY2, five-atomic XY4 molecules are given. In particular, we present the results of calculation for the vibration-nuclear transition probabilities in a case of the emission and absorption spectrum of nucleus 188Os (E=155 keV) in the molecule of OsO4 and nucleus 191Ir (E=82 keV) in the molecule IrO4.
[1]. V.S.Letokhov, V.Minogin, JETP. 69, 1369 (1975); A.Glushkov, L.Ivanov, Phys.Lett.A 170, 33 (1992); A.Glushkov, L.N.Ivanov, V.S.Letokhov, Preprint of ISAN N AS-4, Moscow-Troitsk, (1991);
[2]. A.Glushkov, S.Malinovskaya, O.Khetselius, Europ.Phys.Journ. T160, 195 (2009); Molec.Phys.(UK) 108, 1257 (2008); Int.J. Quant. Chem. 104, 496 (2005); Frontiers in Quantum Systems in Chem. and Phys. (Springer) 18, 523 (2008).
From Koopmans to accurate estimates of molecular core ionization energies
through parameterized allometric scaling
Jean Maruani and Christiane Bonnelle
Laboratoire de Chimie Physique - Matière et Rayonnement
11, rue Pierre et Marie Curie - 75005 Paris, France
Spin Catalysis of Dioxygen Activation by Enzymes
Boris Minaev,1,2 Valentina Minaeva1 and Hans Agren2
1Department of Chemistry, Bogdan Khmelnitskij National University, Cherkasy, Ukraine
2Department of Theoretical Chemistry, Institute of Biotechnology, Royal Institute of Technology, Stockholm, Sweden
In stable organic substances all electron spins are paired and molecules have the singlet ground state. In order to activate chemical transformations and cleave the chemical bond one has to produce spin uncoupling. Interaction with collision partner in bimolecular reactions of diamagnetic species usually leads to activation barrier produced by avoided crossing of the closed shell reactant state and the doubly-triplet excited singlet state. In order to lower the barrier one can add new non-paired electrons to the reacting system or induce spin flip for enhancement of the inter-system crossing to the singly-excited triplet state. This can be done by exchange interaction with paramagnetic transition metal, or by internal (spin-orbit coupling) and external magnetic fields. This is the key idea behind spin-catalysis [1-7].
Activation barriers in chemical reactions are determined by the exchange repulsion between molecules in the singlet ground states. The barriers are often getting much lower when, at least, two spins are unpaired, thus when the singlet-triplet transition occurs. Chemical reactivity is often coded by the triplet excited state of the molecule, in which two spins are unpaired (because exchange interaction destabilizes chemical bond in this case). Enzymatic reactions in live matter are so efficient because they often involve strong spin uncoupling induced by transition metals. This is especially important for O2 (dioxygen) production by photosynthesis of green plants and for dioxygen consumption by respiration of mammals.
Molecular oxygen is a vital elixir beneath the sun. Lavoisier had established that dioxygen is essential for aerobic life and is the oxidant for combustion of organic fuels. But the main puzzle of molecular oxygen is still unclear in modern biochemistry: the O2 molecule has a triplet ground state and how it can react with stable organic substances at 36.7 C to produce water and CO2 if the processes are spin-forbidden? These transformations are equivalent to combustion in the net energetic and material balance. Combustion is a radical chain process. In order to burn a fuel one needs to produce a spark (to create first radical). The radical can interact with triplet O2 and create new radical (the brunch chain reaction). In mammals the organic food is not burned by radical chain process, since the temperature is lower than in the fire. Though the oxidation of food produces enough energy for life, the reactions are spin-forbidden and can be activated and controlled by oxidase enzymes. This should involve particular spin catalysis: the strictly controlled T-S transition in the enzyme active site with low activation energy, which does not produce any radicals, but only short-lived biradical intermediate. The T-S transition means the subsequent spin flip in the biradical.
We present DFT calculations of the electronic structure, zero-field splitting (ZFS) and the phosphorescence spectra of biopolymers that contain chromophores being excited to the triplet state (porphyrin systems, chlorophylls and bacteriochlorophylls, aromatic amino acids, polypeptides and proteins, riboflavin, FAD, FADH2, nicotinamide, NAD+, NADH and active centers of glucose oxidase, copper amine oxidase, hemoglobin, horse-reddish peroxidase and cytochrome P450). This is used for better understanding of the role of spectra in structural analysis of biopolymers and also for predictions of spin catalysis models in electron transfer and oxygen activation processes in these biosystems. The spin transitions are important for enzymatic reactions with and without paramagnetic transition metal involvement. Both types of spin catalysis are studied by DFT calculations of reaction models.
Our experience in fine structure calculations and spin-selective photoprocesses in ZFS triplet states helps us to consider the influence of external magnetic field on the T-S transitions in dark enzymatic reactions. Analysis of possible magnetic field effects on bioprocesses, including birds navigation in the Earth field, is also presented.

1. B.F. Minaev. Electronic mechanisms of molecular oxygen activation. Rus. Chem. Rev., 76 (2007) 998-1023.
2. B.F. Minaev, V.A. Minaeva. Spin-dependent binding of dioxygen to heme. Ukrainica Bioorganica Acta,2 (2008) 56-64.
3. B.F. Minaev. Solvent induced emission of singlet oxygen. J. Mol. Struct. (Theochem), 183 (1989) 207-214.
4. B.F. Minaev. Intensities of spin-forbidden transitions in molecular oxygen and selective heavy-atom effects. Int. J. Quant. Chem., 17 (1980) 367-374.
5. R. Prabhakar, P. Siegbahn, B.F. Minaev, H. Agren. Activation of dioxygen by glucose oxidase. J. Phys. Chem., B, 106 (2002) 3742-3750.
6. B.F. Minaev. Spin effects in reductive activation of O2 by oxidases. RIKEN Rev., 44 (2002) 147-149.
7. B.F. Minaev, H. Agren. Spin catalysis phenomena. Int. J. Quant. Chem., 57 (1996) 510-525.
Intermonomer Interaction Effect on the Electromagnetically Induced Transparency on Molecular Aggregate Model
Takuya Minami and Masayoshi Nakano
Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
The intermonomer interaction effect on the electromagnetically induced transparency (EIT) for a dipole-coupled dimer model is investigated. The absorption properties are evaluated based on the imaginary part of the dynamic polarizability α using the quantum master equation approach [1]. In order to investigate the intermonomer interaction effect on the EIT, several orientations and intermonomer distances are examined. We have found that the EIT can be observed even in the presence of near-degenerete excitation states induced by the intermonomer interaction through an adjustment of the incident field frequency.
[1] M. Nakano et al., J. Phys. Chem. A 105, 22, 5473 (2001); T. Minami et al., J. Phys. Chem. C 114, 13, 6067-6076 (2010).
The potentials of the atoms around Mg2+ in the H-ras GTP complex and in the H-ras GDP complex
Takeshi Miyakawa,1 Ryota Morikawa1, Masako Takasu1, Kimikazu Sugimori2, Kazutomo Kawaguchi3, Hiroaki Saito3 and Hidemi Nagao3
1School of Life Sciences, Tokyo University of Pharmacy and Life Sciences, Japan

2Department of Physical Therapy, Kinjo University, Japan

3Institute of Science and Engineering, Kanazawa University, Japan
In the H-ras GTP complex and the H-ras GDP complex, the coordination bonds between Mg2+ and oxygen atoms are very important. We derived the potentials of atoms around Mg2+ in H-ras GTP complex and in the H-ras GDP complex by quantum chemical calculations in order to modify the AMBER fields ff03.
Because it is not clear if SER17 and THR35 have OH group around Mg2+, we consider two candidate subsystems. In one system, SER17 and THR35 have OH group, although in the other system, SER17 and THR35 have O- instead of OH group. The details will be shown at the presentation.
Validation of quantum chemical methods for geometrical optimizations of sulfonamide derivatives
Akifumi Oda,1,2 Yu Takano2 and Ohgi Takahashi1
1Faculty of Pharmaceutical Sciences, Tohoku Pharmaceutical University, Japan
2Institute for Protein Research, Osaka University, Japan
Sulfonamide is one of the most important chemical groups in drug design because sulfonamide derivatives are stable in living cells and water-soluble. In this study, we assessed the validity of quantum chemical methods and basis sets for the geometrical parameters of various sulfonamides compared to crystallographic data. Introducing f-type polarization functions into basis sets improved the geometry optimizations using Hartree-Fock, MP2, and B3LYP, indicating that f-type polarization functions play an important role in the description of chemical bonds in sulfonamide derivatives.
Theoretical investigation of the hetero-junction effect in polymer stabilized precious metal clsuters
Mitsutaka Okumura,1,2 Yasutaka Kitagawa,1 Takashi Kawakamai,1and Shushuke Yamanaka1
1Department of Chemistry, Graduate School of Science, Osaka University, Japan

2Core Research for Evaluational Science and Technology, Japan
In recent years, gold is attracting industrial and scientific interests for its catalytic activity in such reactions as propylene partial oxidation, odor decomposition H2O2 direct production and CO oxidation at low temperatures, especially when Au is deposited as nanoparticles on selected metal oxides. Highly dispersed gold catalysts exhibit unique catalytic features in low-temperature CO oxidation. In this reaction, the use of oxides of several 3d transition metals and the hydroxides of alkaline earth metals as supports leads to high activities even at a temperature as low as 203 K.
Lately, Tsukuda has found that the Au nanoclusters stabilized by poly(N-Vinyl-2-pyrrolidone) [PVP; (C6H9ON)n] , abbreviated as Au:PVP, can oxidize p-hydroxybenzyl alcohol selectively into the corresponding aldehyde in water without degradation.
This observation indicates that Au cluster can exhibit high catalytic activity without any metal oxide supports. Therefore, the interaction (hetero-junction) between Au cluster and PVP seems to be an important factor to oxidize benzylic alcohols to the corresponding aldehyde over the Au:PVP catalyst.
These results suggest that the hetero-junctions, such as metal cluster-support junctions in heterogeneous catalysts, metal-ligands junctions in metal complex, active sites in enzymes, metal cluster-stabilizer junction in the polymer stabilized clusters and etc., are the key factors to promote and/or modify the catalytic reaction over several types of catalysts. Therefore, we have carried out hybrid density functional theory (DFT) calculations on the hetero-junction effect between Au clusters and PVP. The present theoretical study has been undertaken to explain the hetero-junction effect between PVP molecules and Au clusters as a first step for understanding the hetero-junction effect of the catalytic reactions over the polymer stabilized Au nanoclusters.
Origin of the variety of the Cu2S2 core structure of the CuA sites: a density functional theory study
Orio Okuyama1, Yasuteru Shigeta1,2, Haruki Nakamura1, and Yu Takano1

1Institute for Protein Research, Osaka University, Japan
2Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Japan
The CuA site, a proximate electron-transfer intermediate from an electron source, is contained both in cytochrome c oxidase (C cO), the terminal electron acceptor in aerobic respiration, and in nitrous oxide reductase (N2OR), the terminal acceptor in anaerobic respiration. The CuA site is a binuclear copper site with two bridging cystenyl thiolate groups. Each copper ion is coordinated by a histidine side chain and a methionine side chain or a carbonyl oxygen of the backbone.
Spectroscopic studies have suggested the σu* ground state of the oxidized CuA site is due to the direct Cu–Cu interaction, facilitating a rapid electron transfer over long distances with low driving forces. We have shown that the ground state of the Cu2S2 core is the πu state even with a shorter Cu–Cu distance in the oxidized state, and that a combination of electrostatic and orbital interactions is required to stabilized the σu* state of the Cu2S2 core rather than the πu state [1].
High resolution X-ray crystal structures of the CuA sites in several proteins have so far been reported. The three-dimensional structures of the Cu2S2 core of the CuA sites show various Cu–Cu or S–S distances, 2.35–2.59 Å for the Cu–Cu distance and 3.76–4.24 Å for the S–S distance. In order to elucidate the origin of the variety of the Cu2S2 core structure of several CuA sites, we have examined electronic structures of the Cu2S2 core of the CuA sites, using the density functional theory (DFT) with Gaussian09 program.
We built models from 10 types of C cO (PDB ID: 1V54, 1V55, 2CUA, 2GSM, 2DYR, 2EIJ, 3ABK, 3ABM, 3AG2, 3AG3), 3 types of N2OR (PDB ID: 1FWX, 2IWF, 2IWK), and an engineered CuA azurin (PDB ID: 1CC3). DFT calculations were performed using the exchange-correlation functional, M06, because the M06 functional reproduces the energy difference of the Cu2S2 core between the σu* and πu states at the CCSD(T) level of theory. Wachters+f basis set was used for copper ions and the other atoms were treated with 6-311G++(df,pd) basis set.

[1] Y. Takano, Y. Shigeta, K. Koizumi, H. Nakamura, Int. J. Quantum Chem. in press.
Structure and dynamics of glutathione and glutathione-transferaseT2-2: a molecular dynamics study
Yuriko Omae, Hiroaki Saito, Kazutomo Kawaguchi and Hidemi Nagao
Division of Mathematical and Physical Science, Graduate School of Natural Science and Technology, Kanazawa University
Glutathione is an antioxidant concerned with the detoxification metabolism in vivo and has a function to defecate toxic compounds from cells. Glutathione transferase has been known to reduce the reactivity of toxic compounds by catalyzing their conjugation with glutathione. The function of glutathione and glutathione transferase has been through to implicate in detoxification. Thus, the understanding of the binding character between the glutathione and glutathione transferase should be important in the realization of catalytic reaction mechanism and contribute the development of antibiotic drug. In this study, we carry out the molecular dynamics simulation of glutathione transferase T2-2 in the presence and absence of glutathione in water solvent to analyze the dynamical structure of binding site of the ligand-protein complexes. We calculate the binding energy of the glutathione by these simulations. The specific interactions, which contribute to the binding stability of the ligand molecule, between the glutathione and glutathione transferase T2-2 are shown in detail by these analyses.
Determination of Atomic Charges in Molecules and Ions
Valentin Oshchapovsky
Lviv State University of Vital Activity Safety
35 Kleparivska Str., Lviv, 79007, Ukraine
A new method of calculation of the lattice energy of binary ionic crystals of MX type was developed [1]. It enabled to deduce a new universal formula for the lattice energy calculation taking into account only ion radii values without introducing any additional arbitrary factors. It is pointed out that the exactness of Ulat calculation depends upon the ionnity bond degree [1].
The above-mentioned allowed deducting an equation for a priori calculation of the length of interatomic distances in crystals and gaseous molecules assuming that all the bonds are of pure ionic type [2].
Furthermore, the equation for the calculation of ion radii with an arbitrary effective charge was deduced [3]. The ion radii in the binary non-polar molecules of halogens, chalcogens etc. were calculated: e.g., RIV(F+) = 0.255 Å, RIV(Br+) = 0.918 Å, RIV(At+) = 1.147 Å .
On the basis of these values, during the use of the previously received equation for a priori determination of interatomic distances [2] there were calculated R12 values for the large group of molecules of different type: halogens, interhalogenides, chalcogens, nitrogens and their combinations. It allowed to solve a reverse side of the problem i.e., to estimate the atomic charges of the large group of binary gaseous molecules according to the value of internuclear distance: Î2 (O1.92+O1.92–), ÑÎ (C1.76+C1.76–), CO2 (C3.68+O1.84–), N2 (N2.69+N2.69–), NO (N1.87+O1.87–), NO2 (N3.83+O1.915–).
Considering the experimental R12 values there were also calculated the values of atomic charges in the ionized molecules, e.g. Î2+ (O2.87+O1.87–) and Î2 (O0.97+O1.97–). The radius of He+ ion was calculated. The possibility of estimation of ion co-ordination, bond multiplicity as well as some correction of ion radii R sizes is shown.

1. Oshchapovskii V.V. Russian J. of General Chemistry, 2008, V.78, No. 4, p.532-542.
2. Oshchapovskii V.V. Russian J. of Inorganic Chemistry, 2010, V.55, No. 3, p. 401-409.
3. Oshchapovskii V.V. Russian J. of Inorganic Chemistry (in press).

Electronic structure calculations on vacancy-manganese center in nanodiamond systems
Takao Otsuka,1 Yoshitaka Tateyama,2,3 Masahito Morita,4 Makoto Taiji1
1RIKEN Quantitative Biology Center (QBiC), Kobe, Hyogo, Japan
2International Center for Materials Nanoarchitectonics (MANA), National Institute for Materials Science (NIMS), Tsukuba, Ibaraki, Japan
3CREST, Japan Science and Technology Agency (JST), Kawaguchi, Saitama, Japan
4WPI Immunology Frontier Research Center (iFReC), Osaka University, Suita, Osaka, Japan
We present our theoretical study on vacancy-manganese center in nanodiamond systems by electronic structure calculations. We investigate the stability of spin states in some vacancy-manganese complexes in nanodiamond by using cluster models calculations. From the analysis of electronic structures by density functional theory (DFT) and Hartree-Fock with the perturbative correction calculation (with unrestricted/restricted open-shell), we found that the high spin state of [V-Mn-V] complex is most probable among different spin states and defect types. The results on the surface oxidation systems will be also reported.
Calculation of Energy Deposition by Swift Ions in Biomolecules: Glycine to DNA
John R. Sabin,1,2 Stephan P. A. Sauer,{3 and Jens Oddershede1,2
1Department of Physics and Chemistry, University of Southern Denmark, Odense Denmark
2Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida, USA
3Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
The effects of energy transfer from swift ion radiation to biomolecules are best described by the stopping cross section of the target molecule for the projectile ion. In turn, the mean excitation energy of the target is the determining factor in the stopping cross section. Using polarization propagator methodology, the mean excitation energies of components of several biomolecular systems, ranging from amino acids to nucleotides have been calculated, and are reported here. The calculated mean excitation energies could then be used to determine the stopping cross sections of the various biomolecular systems.

Free energy profile of water across cholesterol containing lipid bilayer
Hiroaki Saito1,3 and Wataru Shinoda2,3
1 Faculty of Mathematics and Physics, Institute of Science and Engineering, Kanazawa University, Kanazawa 920-1192 Japan
2 Nanosystem Research Institute (NRI), Research Institute of Advanced Industrial Science and Technology (AIST), Central 2, 1-1-1 Umezono, Tsukuba 305-8568, Japan
3 CREST, Japan Science and Technology Agency (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
The free energy profiles of a water molecule across two different lipid bilayers of dipalmitoyl phosphatidylcholine (DPPC) and palmitoyl sphingomyelin (PSM) in the absence and presence of cholesterol (0-50 mol%) have been studied by molecular dynamics simulations to elucidate the molecular mechanism of the reduction in water leakage across the membranes by the addition of cholesterol. An enhanced free energy barrier was observed in these membranes with increased cholesterol concentration, and this was explained by the reduced cavity density around the cholesterol in the hydrophobic membrane core. There was an increase of trans conformers in the hydrophobic lipid chains adjacent to the cholesterol, which determined the cavity density. At low cholesterol concentrations the PSM bilayer exhibited a higher free energy barrier than the DPPC bilayer for water permeation, while at greater than 30 mol% of cholesterol the difference became minor. This tendency for the PSM and DPPC bilayers to resemble each other at higher cholesterol concentrations is similar to commonly observed trends in several structural properties, such as order parameters, cross-sectional area per molecule, and cavity density profiles in the hydrophobic regions of bilayer membranes. These results demonstrate that DPPC and PSM bilayers with high cholesterol contents possess similar physical properties, which suggests that the solubility of cholesterol in these lipid bilayers has importance for an understanding of multicomponent lipid membranes with cholesterol.
Modifications for the electron and the polarization propagator formalisms
Masaaki Saitow1, Tomonori Ida2, and Yuji Mochizuki1
1 Department of Chemistry and Research Center for Smart Molecules, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
2 Department of Chemistry, Kanazawa University, Kakuma-machi, Kanazawa 920-1192, Japan
We have derived some improvements for the ab initio Green's function formalisms including the electron and the polarization propagators. For the former case, which describes one-electron detachment process, a spin adapted formalism was derived by considering spin-coupling of the configuration space explicitly. And this approach gives improved description for the open-shell species over the conventional method, although the working equation become rather cumbersome because of the structure of the spin functions. And for the polarization propagator theory, two distinct modifications for the second order framework are shown. The first approach is specifically effective for the Rydberg excitation, in which drastic rearrangement of the charges occur, and recognizable as a branch of denominator shift approach. The second approach improves the differential correlation in terms of the ground and excited states and gives good results for the strongly degenerated states. Both methods retain the polynomial scaling to O(N^5) .

Evaluation of multicenter integrals of Slater type orbitals and Coulomb-Yukawa like electric field gradient potentials using their one-range addition theorems
Faculty of Science, Department of physics, Trakya university, Edirne, Turkey
Using one-range addition theorems established in Ref.[1] for Slater type orbitals and Coulomb-Yukawa like interaction potentials with the help of complete orthonormal sets of ψ^α-exponential type orbitals (α=1,0,-1,-2,…) [2], the calculations are performed for multicenter electric field gradient integrals. The results of computations show good rate of convergence and numerical stability. The convergences of the series is tested by calculating multicenter electric field gradient integrals for the arbitrary values of parameters of potentials and location of orbitals.

[1] I. I. Guseinov, Bull. Chem. Soc. Jpn., 78 (2005) 611
[2] I. I. Guseinov, Int. J. Quantum Chem. ,90 (2002) 114
Cluster Embedding Method with Non-Orthogonal Wave Functions: DFT Kohn-Sham Embedding Equations and Vacant States
Emma K. Shidlovskaya1,2
1Information Systems Management Institute, Riga, Latvia

2Institute of Chemical Physics, University of Latvia, Riga, Latvia
When we theoretically study processes in large or infinite electron systems we have to treat the whole quantum system as two subsystems: small fragment of the system (cluster) and the remaining part of it. Problem "cluster in the field of the rest of system" is successfully solved in the frameworks of embedded molecular cluster (EMC) model [1] with orthogonal wave functions. Unfortunately, standard realization of EMC model leads to well-pronounced boundary effects [2].

To overcome limitations of the standard EMC model, we have treated cluster embedding problem in the frameworks of one-electron approximation with non-orthogonal wave functions. Equations for the cluster are obtained varying total energy of the system expressed in terms of non-orthogonal one-electron wave functions. Using these cluster embedding equations we have developed modified cluster embedding scheme. We have demonstrated that application of this scheme may radically reduce boundary effects in EMC model [2].

Possibility to generalize our embedding approach on the case of DFT Kohn-Sham one-electron equations is studied. We demonstrate that our variation procedure is compatible with Kohn-Sham method. Cluster embedding equations remain the same [3] if instead of Fock operator we use one-electron Kohn-Sham Hamiltonian.

For further applications of our cluster embedding method we should overcome limitations of one-electron approximation. It may be done by configurations interaction (CI) or perturbation theory (PT) methods. For this purpose we need occupied and vacant cluster states of the same localization. We have established that our initial embedding equations [2] give localized in the cluster region occupied states and delocalized vacant ones [4]. To get the same localization degree for the both occupied and vacant states, modified equations are proposed [4].

Modified cluster embedding equations [4] lead to correct description of electron transitions from occupied states to vacant ones. Treatment of electron correlation effects by CI or PT and proper description of excited states both by CI or DFT become possible. As the result, our embedding scheme now may be applied for quantum-chemical simulation of various phenomena.

[1] L.N.Kantorovich, J. Phys. C: Solid State Phys. 21, 5041 (1988).
[2] E.K. Shidlovskaya, Int. J. Quantum Chem. 89, 349 (2002).
[3] E.K. Shidlovskaya, in Topics in Chemistry and Material Science, Volume 2. Theoretical Aspects of Catalysis, eds. G. Vayssilov, T. Mineva, Heron Press, Sofia, 2009, pp. 11-18.
[4] E.K. Shidlovskaya, Computer Modelling and New Technologies 10, No 4, 17 (2006).

Real-Time TDDFT Simulation for Coherent Phonon Generation
Yasushi Shinohara1, Kazuhiro Yabana1,2, Jun-Ichi Iwata2, Tomohito Otobe3 and George F. Bertsch4
1Graduate School for Pure and Applied Sciences, University of Tsukuba, Japan

2Center for Computational Sciences, University of Tsukuba, Japan

3Advanced Photon Research Center, Japan Atomic Energy Agency, Japan

4Department of Physics, University of Washington, USA
We have been developing a first-principles description for quantum dynamics of electrons and ions induced by an intense and ultrashort laser pulse. In our method, electron dynamics is described by the time-dependent density-functional theory (TDDFT) solving the time-dependent Kohn-Sham equation in real-time, while the ion motion is described by classical mechanics with the force evaluated by the Ehrenfest theorem. We consider a dynamics under a spatially-uniform electric field which includes the polarization field as well as the applied electric field. Our scheme is capable of describing dielectric function in the linear response regime [1]. It has also been applied to optical breakdown in nonlinear regime [2].
In my presentation, we report an application of our framework for coherent phonon generation in bulk Si [3]. Coherent phonon is a macroscopic lattice vibration induced by the ultrashort laser pulse whose duration is shorter than the vibration period.
The figure below shows the calculated force acting on optical phonon coordinate at photon frequencies of 2.25, 2.5 and 2.75 eV. In our calculation, direct band gap of bulk Si is 2.4 eV. When the laser frequency is below the direct band gap (2.25 eV, red solid line), the impulsive force is only seen during the irradiation of the laser pulse (0-16 fs). When the laser frequency is above the direct band gap (2.75 eV, blue dotted line), the force persists even after the laser pulse ends. These features are consistent with two mechanisms of coherent phonon generation which have been proposed in phenomenological studies, the impulsive stimulated Raman scattering (ISRS) and the diplacive excitation of coherent phonon (DECP). In the ISRS, the force is generated from virtual electron-hole excitations. In the DECP, the force is generated from real electron-hole excitations which persist even after the laser pulse ends. We thus conclude that, in our framework based on TDDFT, both mechanisms of ISRS and DECP are included in a unified way.
We also apply the method to semimetal, Sb. Preliminary results will be presented as well.

[1] G.F. Bertsch, J.-I. Iwata, A. Rubio and K. Yabana, Phys. Rev. B 62, 7998 (2000)
[2] T. Otobe, M. Yamagiwa, J.-I. Iwata, K. Yabana, T. Nakatsukasa and G.F. Bertsch, Phys. Rev. B 77,165104 (2008)
[3] Y. Shinohara, K. Yabana, Y. Kawashita, J.-I. Iwata, T. Otobe and G.F. Bertsch, Phys. Rev. B 82, 155110 (2010)

Figure. The force on the optical phonon coordinate for three laser frequencies: 2.25 eV (red solid); 2.5 eV (green dashed); and 2.75 eV (blue dotted).
Theoretical study of isotope-induced additivity of chemical shift in benzene
Kimikazu Sugimori,1 Hiroyuki Kawabe2, and Hideto Shimahara3
1Department of Physical Therapy, Faculty of Health Sciences, Kinjo University, Japan

2Department of Social Work, Faculty of Social Work, Kinjo University, Japan

3Center for Nano Materials and Technology, Japan Advanced Institute of Science and Technology, Japan
The isotope effect induced by deuterium substituted species is observed in molecular properties, such as geometry, kinetics, and electronic state of the molecules through nuclear-electron interaction. Theoretical considerations and experimental alignments have been studied by ab initio molecular orbital (MO), density functional (DF) theory, and other empirical strategies. The Born-Oppenheimer (BO) approximation with nuclear vibrational wavefunction can treat isotope effect because nuclear mass effect account for the average distance of vibrational motion [1]. We have been reported for diatomic molecules [2].
In this study, we apply Morse anharmonic oscillator model to calculate average internuclear distance of isotopologues and isotopomers of methan-like molecules and benzene derivatives. NMR shielding constants are calculated again on the average distance by using Gauge-Independent Atomic Orbital (GIAO) method with B3LYP/cc-pVTZ and CCSD/cc-pVTZ. The calculated isotope shifts will be discussed with relation to additivity, charge distribution, and experimental values.

[1] A. C. de Dios and C. J. Jameson. Annu. Rep. NMR Spectrosc. 29 (1994) 1-69.
[2] K. Sugimori and H. Kawabe, Int. J. Quantum Chem. 110 (2010) 2989-2995. doi: 10.1002/qua.22917
Development of First-Principles Maxwell+TDDFT Multi-Scale Simulator for
Propagation of High-Intensity Laser Pulse
Takeru Sugiyama1,Yasushi Shinohara1,Tomohito Otobe2,Kazuhiro Yabana1,3 and George F. Bertsch4
1Graduate School for Pure and Applied Sciences, University of Tsukuba, Japan
2Center for Computational Sciences, University of Tsukuba, Japan
3Advanced Photon Research Center, Japan Atomic Energy Agency, Japan
4Department of Physics, University of Washington, USA
Interaction between light and matter is described by the Schroedinger and Maxwell equations. The Schroedinger equation describes electron dynamics while the Maxwell equation describes propagation of electromagnetic fields. For ordinary weak light-wave, one can apply the perturbation theory for the Schroedinger equation which decouples two equations with the dielectric function. However, for intensive laser pulses, one cannot separate them because of the nonlinear electron responses to the strong electric field of the laser pulse. Previously, we developed a framework in TDDFT to describe electron dynamics under spatially-uniform time-varying electric field solving the time-dependent Kohn-Sham equation in real-time [1,2]. We now extend it to a first-principles simulator calculating simultaneously the coupled nonlinear dynamics of electrons and electromagnetic field. Since the length-scale is much different between the laser wavelength (μm) and the electron dynamics (nm), we employ two different spatial grids and express the vector potential in the macroscopic grids and the Kohn-Sham orbitals in the microscopic grids. As a preliminary demonstration, we will show our calculation for the one dimensional propagation of electromagnetic field incident on bulk Si.

[1]T. Otobe et al. Phys.Rev.B77,165104(2008)
[2]Y. Shinohara et al. Phys.Rev.B82,155110(2010)

Figure. Laser pulse irradiated on bulk Si surface. The intensity and the frequency of the laser pulse is set to I=5×1012W/cm2 and ℏw=1.55eV (below calculated band-gap, 2.4eV), respectively. The upper panels show propagation of electromagnetic field. The lower panels show the ground-state electron density (left) and the density change of electrons from that in the ground state at the surface (middle and right).

L. Sandoval*, I. Urdaneta and A. Palma
*Facultad de Ciencias de la Computacion, Benemerita Universidad Autonoma
de Puebla. Puebla, Pue. 72570, Mexico

Instituto de Fisica, Benemerita Universidad Autonoma de Puebla. Puebla, Pue. 72570, Mexico
In this work we develop a harmonic approximation of the Franck-Condon Factors (FCF) for any anharmonic potential of diatomic molecules.

The method is based on the Taylor series expansion of the potential and the second quantization formalism. Well known recurrence relations are used to calculate the FCF, which are incorporated in the hamiltonian. Keeping then only the quadratic terms, by the Bogoliubov-Tyablikov transformation the hamiltonian is found to be equivalent to an harmonic oscillator. The derivation is an alternative route to the Iterative Bogoliubov transformation (IBT), widely used to treat anharmonic potentials in a non-perturbative way. The FCF obtained by our method for the Morse potential are compared with numerical techniques like the RKR method and the Morse itself. Our results are in agreement with these numerical methods, being our technique entirely analytical and much simpler to use.
Quantum Chemistry on Quantum Computers
L. Veis and J. Pittner
J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague, Czech Republic
Quantum computers are appealing for their ability to solve some tasks much faster than their classical counterparts, e.g. efficiently factore integers. Quantum chemistry could in principle benefit from them as well, for example by an efficient solution of many-body Hamiltonian eigenvalue problem [1]. As was shown in the seminal work by Aspuru-Guzik et. al. [2], quantum computers, if available, would be able to perform the full configuration interaction (FCI) energy calculations with only a polynomial scaling, in contrast to conventional computers where FCI scales exponentially.

We have developed a code for simulation of quantum computers and implemented our version of the quantum full configuration interaction (QFCI) method which uses the iterative phase estimation algorithm. This approach reduces demands on the total number of quantum bits (qubits) as only one is needed in the read-out part of the quantum register and the whole algorithm proceeds in an iterative manner.

We have tested its performance and applicability for non-relativistic as well as relativistic CI energy calculations. Non-relativistic QFCI calculations of the four lowest lying electronic states of methylene molecule (CH2), which exhibit a multireference character were performed [3]. It has been shown that with a suitably chosen initial state of the quantum register, one is able to achieve the probability amplification regime of the iterative phase estimation even for nearly dissociated molecule. Relativistic Kramers-restricted CI calculations employing the QFCI algorithm have been applied to the spin-orbit coupling in the SbH molecule [4]. We have also designed the quantum circuits for the simplest proof-of-principle physical realizations of relativistic quantum chemical computations on quantum computers.

[1] Abrams, D. S.; Lloyd, S. Phys.Rev.Lett. 1999, 83, 5162–5165.
[2] Aspuru-Guzik, A.; Dutoi, A. D.; Love, P. J.; Head-Gordon, M. Science 2005, 309, 1704–1707.
[3] Veis, L.; Pittner, J. J. Chem. Phys. 2010, 133, 194106.
[4] Veis, L; Višnák, J.; Fleig, T.; Knecht, S.; Saue, T.; Visscher, L.; Pittner, J. in preparation
Antidot Structure Dependences of Open-shell Characters and Aromaticities for Hexagonal Graphene Nanoflakes
Kyohei Yoneda,1 Yudai Inoue,1 Tomoya Inui,1 Yasuteru Shigeta,1 Takashi Kubo,2 Benoît Champagne3 and Masayoshi Nakano1
1Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Japan
2Department of Chemistry, Graduate School of Science, Osaka University, Japan
3Laboratoire de Chimie Théorique (LCT), Facultés Universitaires Notre-Dame de la Paix (FUNDP), BELGIUM
The open-shell singlet state of graphene nanoflakes is at the origin of their unique physico-chemical properties. In this study, we theoretically investigate the open-shell characters and aromaticities of hexagonal graphene nanoflakes with different sizes of antidot structures, using the long-range corrected spin-unrestricted density functional theory, LC-UBLYP, method. It is found that the open-shell character exhibits an oscillatory behavior with increasing the size of the antidot structure, and the nucleus-independent chemical shift (NICS) – an index of aromaticity – also depends on the open-shell character. These antidot structure dependences are rationalized in terms of the variations in the HOMO–LUMO energy gaps, which is explained by the molecular orbital correlation diagram.
Representation of structure-property relationships in polymorphic systems
Koretaka Yuge
Department of Materials Science and Engineering, Kyoto University, Kyoto 606-8501, Japan
While structure-property relationships in crystals are fundamental prerequisite to design and optimize desirable materials, quantitative assessment for the relationships has been just beginning. Cluster expansion (CE) is one of the most successful approach based on ab-initio calculation to predict configurational scalar as well as tensor properties such as internal energy, bandgap, density of states, and elasticity. However, application of the current CE is essentially limited to the single lattice, which makes it really difficult to applying to the polymorphic systems.
Here we introduce variable lattice CE (VLCE)[1] enabling to treat multiple structures in a single Hamiltonian, which overcomes the limitation in CE. We demonstrate derivation, concept, and interpretation of the proposed VLCE, and show application to 2-dimensional system to search ground-state structure in polymorphs.

[1] K. Yuge, J. Phys.: Condens. Matter 22, 125402 (2010).
Use of one-range addition theorems in evaluation of multicenter electric field integrals of Slater type orbitals and Yukawa-like interaction potential
Nimet ZAIM
Department of Physics, Faculty of Sciences, Trakya University, Edirne, Turkey
The calculation is performed for multicenter electric field integrals containing Slater type orbitals and Yukawa-like interaction potential using their one-range addition theorems established in Ref.[1] with the help of complete orthonormal sets of ψα- exponantial type orbitals (α=1,0,-1,-2,…) [2]. The results of computer calculations are presented. The convergences of the series is tested by calculating concrete cases for the arbitrary values of quantum numbers, orbital parameters and internuclear distances.

[1] I. I. Guseinov, Bull. Chem. Soc. Jpn., 78 (2005) 611
[2] I. I. Guseinov, Int. J. Quantum Chem. ,90 (2002) 114
Mn-doped Thiolated Au25 Nanoclusters: Atomic Confguration, Magnetic Properties, and A Possible High-performance Spin Filter
M. Zhou,1 Y. Q. Cai,1 M. G. Zeng,1 C. Zhang,1,2 and Y. P. Feng1
1Department of Physics, National University of Singapore, 2 Science Drive 3, Singapore, 117542
2Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore, 117543
We report an ab inito investigation on the ground-state atomic configuration, electronic structures, magnetic and spin-dependent transport properties of Mn-doped Au25 nanoclusters protected by thiolate. It is found that the most stable dopant sites are near surfaces, rather than the center position of the nanoparticles. Transport calculations show that high-performance spin filters can be achieved by sandwiching these doped clusters between two nonmagnetic Au electrodes. The nearly perfect spin filtering ori g! inates from localized magnetic moments of these clusters that are well protected by ligands from the presence of electrodes

 Main | Registration | Abstract Submission | List of Registrants | List of Abstracts | Program 
Accommodation | Location and Access | Previous Workshops | ISTCP-VII | Photos