1
|
Zhang C, Lipparini F, Stopkowicz S, Gauss J, Cheng L. Cholesky Decomposition-Based Implementation of Relativistic Two-Component Coupled-Cluster Methods for Medium-Sized Molecules. J Chem Theory Comput 2024; 20:787-798. [PMID: 38198515 DOI: 10.1021/acs.jctc.3c01236] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
A Cholesky decomposition (CD)-based implementation of relativistic two-component coupled-cluster (CC) and equation-of-motion CC (EOM-CC) methods using an exact two-component Hamiltonian augmented with atomic-mean-field spin-orbit integrals (the X2CAMF scheme) is reported. The present CD-based implementation of X2CAMF-CC and EOM-CC methods employs atomic-orbital-based algorithms to avoid the construction of two-electron integrals and intermediates involving three and four virtual indices. Our CD-based implementation extends the applicability of X2CAMF-CC and EOM-CC methods to medium-sized molecules with the possibility to correlate around 1000 spinors. Benchmark calculations for uranium-containing small molecules were performed to assess the dependence of the CC results on the Cholesky threshold. A Cholesky threshold of 10-4 is shown to be sufficient to maintain chemical accuracy. Example calculations to illustrate the capability of the CD-based relativistic CC methods are reported for the bond-dissociation energy of the uranium hexafluoride molecule, UF6, with up to quadruple-ζ basis sets, and the lowest excitation energy in the solvated uranyl ion [UO22+(H2O)12].
Collapse
Affiliation(s)
- Chaoqun Zhang
- Department of Chemistry, the Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa I-56124, Italy
| | - Stella Stopkowicz
- Fachrichtung Chemie, Universität des Saarlandes, Saarbrücken D-66123, Germany
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, Oslo N-0315, Norway
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, Mainz D-55128, Germany
| | - Lan Cheng
- Department of Chemistry, the Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
2
|
Chamoli S, Surjuse K, Jangid B, Nayak MK, Dutta AK. A reduced cost four-component relativistic coupled cluster method based on natural spinors. J Chem Phys 2022; 156:204120. [PMID: 35649878 DOI: 10.1063/5.0085932] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
We present the theory, implementation, and benchmark results for a frozen natural spinors based reduced cost four-component relativistic coupled cluster method. The natural spinors are obtained by diagonalizing the one-body reduced density matrix from a relativistic second-order Møller-Plesset calculation based on a four-component Dirac-Coulomb Hamiltonian. The correlation energy in the coupled cluster method converges more rapidly with respect to the size of the virtual space in the frozen natural spinor basis than that observed in the standard canonical spinors obtained from the Dirac-Hartree-Fock calculation. The convergence of properties is not smooth in the frozen natural spinor basis. However, the inclusion of the perturbative correction smoothens the convergence of the properties with respect to the size of the virtual space in the frozen natural spinor basis and greatly reduces the truncation errors in both energy and property calculations. The accuracy of the frozen natural spinor based coupled cluster methods can be controlled by a single threshold and is a black box to use.
Collapse
Affiliation(s)
- Somesh Chamoli
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Kshitijkumar Surjuse
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Bhavnesh Jangid
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
3
|
Helmich-Paris B. A trust-region augmented Hessian implementation for state-specific and state-averaged CASSCF wave functions. J Chem Phys 2022; 156:204104. [DOI: 10.1063/5.0090447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we present a one-step second-order converger for state-specific (SS) and state-averaged (SA) complete active space self-consistent field (CASSCF) wave functions. Robust convergence is achieved through step restrictions using a trust-region augmented Hessian (TRAH) algorithm. To avoid numerical instabilities, an exponential parameterization of variational configuration parameters is employed, which works with a nonredundant orthogonal complement basis. This is a common approach for SS-CASSCF and is extended to SA-CASSCF wave functions in this work. Our implementation is integral direct and based on intermediates that are formulated in either the sparse atomic-orbital or small active molecular-orbital basis. Thus, it benefits from a combination with efficient integral decomposition techniques, such as the resolution-of-the-identity or the chain-of-spheres for exchange approximations. This facilitates calculations on large molecules, such as a Ni(II) complex with 231 atoms and 5154 basis functions. The runtime performance of TRAH-CASSCF is competitive with the other state-of-the-art implementations of approximate and full second-order algorithms. In comparison with a sophisticated first-order converger, TRAH-CASSCF calculations usually take more iterations to reach convergence and, thus, have longer runtimes. However, TRAH-CASSCF calculations still converge reliably to a true minimum even if the first-order algorithm fails.
Collapse
Affiliation(s)
- Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
4
|
Helmich-Paris B, de Souza B, Neese F, Izsák R. An improved chain of spheres for exchange algorithm. J Chem Phys 2021; 155:104109. [PMID: 34525816 DOI: 10.1063/5.0058766] [Citation(s) in RCA: 66] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
In the present work, we describe a more accurate and efficient variant of the chain-of-spheres algorithm (COSX) for exchange matrix computations. Higher accuracy for the numerical integration is obtained with new grids that were developed using global optimization techniques. With our new default grids, the average absolute energy errors are much lower than 0.1 kcal/mol, which is desirable to achieve "chemical accuracy." Although the size of the new grids is increased by roughly a factor of 2.5, the excellent efficiency of the original COSX implementation is still further improved in most cases. The evaluation of the analytic electrostatic potential integrals was significantly accelerated by a new implementation of rolled-out versions of the Dupuis-Rys-King and Head-Gordon-Pople algorithms. Compared to our earlier implementation, a twofold speedup is obtained for the frequently used triple-ζ basis sets, while up to a 16-fold speedup is observed for quadruple-ζ basis sets. These large gains are a consequence of both the more efficient integral evaluation and the intermediate exchange matrix computation in a partially contracted basis when generally contracted shells occur. With our new RIJCOSX implementation, we facilitate accurate self-consistent field (SCF) binding energy calculations on a large supra-molecular complex composed of 320 atoms. The binding-energy errors with respect to the fully analytic results are well below 0.1 kcal/mol for the cc-pV(T/Q)Z basis sets and even smaller than for RIJ with fully analytic exchange. At the same time, our RIJCOSX SCF calculation even with the cc-pVQZ basis and the finest grid is 21 times faster than the fully analytic calculation.
Collapse
Affiliation(s)
- Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
5
|
Pototschnig JV, Papadopoulos A, Lyakh DI, Repisky M, Halbert L, Severo Pereira Gomes A, Jensen HJA, Visscher L. Implementation of Relativistic Coupled Cluster Theory for Massively Parallel GPU-Accelerated Computing Architectures. J Chem Theory Comput 2021; 17:5509-5529. [PMID: 34370471 PMCID: PMC8444343 DOI: 10.1021/acs.jctc.1c00260] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this paper, we report reimplementation of the core algorithms of relativistic coupled cluster theory aimed at modern heterogeneous high-performance computational infrastructures. The code is designed for parallel execution on many compute nodes with optional GPU coprocessing, accomplished via the new ExaTENSOR back end. The resulting ExaCorr module is primarily intended for calculations of molecules with one or more heavy elements, as relativistic effects on the electronic structure are included from the outset. In the current work, we thereby focus on exact two-component methods and demonstrate the accuracy and performance of the software. The module can be used as a stand-alone program requiring a set of molecular orbital coefficients as the starting point, but it is also interfaced to the DIRAC program that can be used to generate these. We therefore also briefly discuss an improvement of the parallel computing aspects of the relativistic self-consistent field algorithm of the DIRAC program.
Collapse
Affiliation(s)
- Johann V Pototschnig
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, de Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Anastasios Papadopoulos
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, de Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Dmitry I Lyakh
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Michal Repisky
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Loïc Halbert
- Universite de Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers, Atomes et Molecules, F-59000 Lille, France
| | - André Severo Pereira Gomes
- Universite de Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers, Atomes et Molecules, F-59000 Lille, France
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, de Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
6
|
Szabó PB, Csóka J, Kállay M, Nagy PR. Linear-Scaling Open-Shell MP2 Approach: Algorithm, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2021; 17:2886-2905. [PMID: 33819030 PMCID: PMC8154337 DOI: 10.1021/acs.jctc.1c00093] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A linear-scaling
local second-order Møller–Plesset
(MP2) method is presented for high-spin open-shell molecules based
on restricted open-shell (RO) reference functions. The open-shell
local MP2 (LMP2) approach inherits the iteration- and redundancy-free
formulation and the completely integral-direct, OpenMP-parallel, and
memory and disk use economic algorithms of our closed-shell LMP2 implementation.
By utilizing restricted local molecular orbitals for the demanding
integral transformation step and by introducing a novel long-range
spin-polarization approximation, the computational cost of RO-LMP2
approaches that of closed-shell LMP2. Extensive benchmarks were performed
for reactions of radicals, ionization potentials, as well as spin-state
splittings of carbenes and transition-metal complexes. Compared to
the conventional MP2 reference for systems of up to 175 atoms, local
errors of at most 0.1 kcal/mol were found, which are well below the
intrinsic accuracy of MP2. RO-LMP2 computations are presented for
challenging protein models of up to 601 atoms and 11 000 basis
functions, which involve either spin states of a complexed iron ion
or a highly delocalized singly occupied orbital. The corresponding
runtimes of 9–15 h obtained with a single, many-core CPU demonstrate
that MP2, as well as spin-scaled MP2 and double-hybrid density functional
methods, become widely accessible for open-shell systems of unprecedented
size and complexity.
Collapse
Affiliation(s)
- P Bernát Szabó
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csóka
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
7
|
Senjean B, Sen S, Repisky M, Knizia G, Visscher L. Generalization of Intrinsic Orbitals to Kramers-Paired Quaternion Spinors, Molecular Fragments, and Valence Virtual Spinors. J Chem Theory Comput 2021; 17:1337-1354. [PMID: 33555866 DOI: 10.1021/acs.jctc.0c00964] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Localization of molecular orbitals finds its importance in the representation of chemical bonding (and antibonding) and in the local correlation treatments beyond mean-field approximation. In this paper, we generalize the intrinsic atomic and bonding orbitals [G. Knizia, J. Chem. Theory Comput. 2013, 9, 11, 4834-4843] to relativistic applications using complex and quaternion spinors, as well as to molecular fragments instead of atomic fragments only. By performing a singular value decomposition, we show how localized valence virtual orbitals can be expressed on this intrinsic minimal basis. We demonstrate our method on systems of increasing complexity, starting from simple cases such as benzene, acrylic acid, and ferrocene molecules, and then demonstrate the use of molecular fragments and inclusion of relativistic effects for complexes containing heavy elements such as tellurium, iridium, and astatine. The aforementioned scheme is implemented into a standalone program interfaced with several different quantum chemistry packages.
Collapse
Affiliation(s)
- Bruno Senjean
- Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands
| | - Souloke Sen
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Michal Repisky
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Gerald Knizia
- Department of Chemistry, The Pennsylvania State University, University Park, State College, Pennsylvania 16802, United States
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
8
|
Glasbrenner M, Graf D, Ochsenfeld C. Efficient Reduced-Scaling Second-Order Møller-Plesset Perturbation Theory with Cholesky-Decomposed Densities and an Attenuated Coulomb Metric. J Chem Theory Comput 2020; 16:6856-6868. [PMID: 33074664 DOI: 10.1021/acs.jctc.0c00600] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a novel, highly efficient method for the computation of second-order Møller-Plesset perturbation theory (MP2) correlation energies, which uses the resolution of the identity (RI) approximation and local molecular orbitals obtained from a Cholesky decomposition of pseudodensity matrices (CDD), as in the RI-CDD-MP2 method developed previously in our group [Maurer, S. A.; Clin, L.; Ochsenfeld, C. J. Chem. Phys. 2014, 140, 224112]. In addition, we introduce an attenuated Coulomb metric and subsequently redesign the RI-CDD-MP2 method in order to exploit the resulting sparsity in the three-center integrals. Coulomb and exchange energy contributions are computed separately using specialized algorithms. A simple, yet effective integral screening protocol based on Schwarz estimates is used for the MP2 exchange energy. The Coulomb energy computation and the preceding transformations of the three-center integrals are accelerated using a modified version of the natural blocking approach [Jung, Y.; Head-Gordon, M. Phys. Chem. Chem. Phys. 2006, 8, 2831-2840]. Effective subquadratic scaling for a wide range of molecule sizes is demonstrated in test calculations in conjunction with a low prefactor. The method is shown to enable cost-efficient MP2 calculations on large molecular systems with several thousand basis functions.
Collapse
Affiliation(s)
- Michael Glasbrenner
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstrasse 7, 81377 Munich, Germany
| | - Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstrasse 7, 81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstrasse 7, 81377 Munich, Germany
| |
Collapse
|
9
|
Saue T, Bast R, Gomes ASP, Jensen HJA, Visscher L, Aucar IA, Di Remigio R, Dyall KG, Eliav E, Fasshauer E, Fleig T, Halbert L, Hedegård ED, Helmich-Paris B, Iliaš M, Jacob CR, Knecht S, Laerdahl JK, Vidal ML, Nayak MK, Olejniczak M, Olsen JMH, Pernpointner M, Senjean B, Shee A, Sunaga A, van Stralen JNP. The DIRAC code for relativistic molecular calculations. J Chem Phys 2020; 152:204104. [PMID: 32486677 DOI: 10.1063/5.0004844] [Citation(s) in RCA: 118] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.
Collapse
Affiliation(s)
- Trond Saue
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Radovan Bast
- Department of Information Technology, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - André Severo Pereira Gomes
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Ignacio Agustín Aucar
- Instituto de Modelado e Innovación Tecnológica, CONICET, and Departamento de Física-Facultad de Ciencias Exactas y Naturales, UNNE, Avda. Libertad 5460, W3404AAS Corrientes, Argentina
| | - Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Kenneth G Dyall
- Dirac Solutions, 10527 NW Lost Park Drive, Portland, Oregon 97229, USA
| | - Ephraim Eliav
- School of Chemistry, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
| | - Elke Fasshauer
- Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus, Denmark
| | - Timo Fleig
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Loïc Halbert
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Erik Donovan Hedegård
- Division of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Miroslav Iliaš
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, 974 01 Banská Bystrica, Slovakia
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Stefan Knecht
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Jon K Laerdahl
- Department of Microbiology, Oslo University Hospital, Oslo, Norway
| | - Marta L Vidal
- Department of Chemistry, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
| | - Małgorzata Olejniczak
- Centre of New Technologies, University of Warsaw, S. Banacha 2c, 02-097 Warsaw, Poland
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | | | - Bruno Senjean
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Ayaki Sunaga
- Department of Chemistry, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji-city, Tokyo 192-0397, Japan
| | - Joost N P van Stralen
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| |
Collapse
|