1
|
Di Grande S, Barone V. Toward Accurate Quantum Chemical Methods for Molecules of Increasing Dimension: The New Family of Pisa Composite Schemes. J Phys Chem A 2024; 128:4886-4900. [PMID: 38847454 DOI: 10.1021/acs.jpca.4c01673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
The new versions of the Pisa composite scheme introduced in the present paper are based on the careful selection of different quantum chemical models for energies, geometries, and vibrational frequencies, with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, the computation of accurate electronic energies has been further improved introducing more reliable complete basis set extrapolations and estimation of core-valence correlation, together with improved basis sets for third-row atoms. Furthermore, the reduced-cost frozen natural orbital (FNO) model has been introduced and validated for large molecules. Accurate molecular structures can be obtained avoiding complete basis set extrapolation and evaluating core-valence correlation at the MP2 level. Unfortunately, analytical gradients are not available for the FNO version of the model. Therefore, for large molecules, an accurate reduced-cost alternative is offered by evaluation of valence contributions with a double-hybrid functional in conjunction with the same MP2 contribution for core-valence correlation or by means of a one-parameter approximation. The same double-hybrid functional and basis set are employed to evaluate zero-point energies and partition functions. After the validation of the new models for small systems, a panel of molecular bricks of life has been used to analyze their performances for problems of current fundamental or technological interest. The fully black-box implementation of the computational workflow paves the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by experimentally oriented researchers.
Collapse
Affiliation(s)
- Silvia Di Grande
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
| | | |
Collapse
|
2
|
Mendolicchio M, Barone V. Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions. J Chem Theory Comput 2024; 20:2842-2857. [PMID: 38556752 DOI: 10.1021/acs.jctc.4c00023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/02/2024]
Abstract
The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.
Collapse
|
3
|
Lazzari F, Mendolicchio M, Barone V. Accurate Geometries of Large Molecules by Integration of the Pisa Composite Scheme and the Templating Synthon Approach. J Phys Chem A 2024; 128:1385-1395. [PMID: 38347709 DOI: 10.1021/acs.jpca.3c08382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
An effective yet reliable computational workflow is proposed, which permits the computation of accurate geometrical structures for large flexible molecules at an affordable cost thanks to the integration of machine learning tools and DFT models together with reduced scaling computations of vibrational averaging effects. After validation of the different components of the overall strategy, a panel of molecules of biological interest have been analyzed. The results confirm that very accurate geometrical parameters can be obtained at reasonable cost for molecules including up to about 50 atoms, which are the largest ones for which comparison with high-resolution rotational spectra is possible. Since the whole computational workflow can be followed employing standard electronic structure codes, accurate results for large-sized molecules can be obtained at DFT cost also by nonspecialists.
Collapse
Affiliation(s)
- Federico Lazzari
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Mendolicchio
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | |
Collapse
|
4
|
Barone V. Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase. Phys Chem Chem Phys 2024; 26:5802-5821. [PMID: 38099409 DOI: 10.1039/d3cp05169b] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2024]
Abstract
Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| |
Collapse
|
5
|
Barone V, Crisci L, Di Grande S. Accurate Thermochemical and Kinetic Parameters at Affordable Cost by Means of the Pisa Composite Scheme (PCS). J Chem Theory Comput 2023; 19:7273-7286. [PMID: 37774410 PMCID: PMC10601482 DOI: 10.1021/acs.jctc.3c00817] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Indexed: 10/01/2023]
Abstract
A new strategy for the computation at an affordable cost of geometrical structures, thermochemical parameters, and rate constants for medium-sized molecules in the gas phase is proposed. The most distinctive features of the new model are the systematic use of cc-pVnZ-F12 basis sets, the addition of MP2 core-valence correlation in geometry optimizations by a double-hybrid functional, the separate extrapolation of MP2 and post-MP2 contributions, and the inclusion of anharmonic contributions in zero-point energies and thermodynamic functions. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme outperforms the most well-known model chemistries without the need for any empirical parameter. Additional tests show that the computed zero-point energies and thermal contributions can be confidently used for obtaining accurate thermochemical and kinetic parameters. Since the whole computational workflow is translated in a black-box procedure, which can be followed with standard electronic structure codes, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56125 Pisa, Italy
| | - Luigi Crisci
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56125 Pisa, Italy
| | - Silvia Di Grande
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56125 Pisa, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
| |
Collapse
|
6
|
Song S, Vuckovic S, Kim Y, Yu H, Sim E, Burke K. Extending density functional theory with near chemical accuracy beyond pure water. Nat Commun 2023; 14:799. [PMID: 36781855 PMCID: PMC9925738 DOI: 10.1038/s41467-023-36094-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Accepted: 01/13/2023] [Indexed: 02/15/2023] Open
Abstract
Density functional simulations of condensed phase water are typically inaccurate, due to the inaccuracies of approximate functionals. A recent breakthrough showed that the SCAN approximation can yield chemical accuracy for pure water in all its phases, but only when its density is corrected. This is a crucial step toward first-principles biosimulations. However, weak dispersion forces are ubiquitous and play a key role in noncovalent interactions among biomolecules, but are not included in the new approach. Moreover, naïve inclusion of dispersion in HF-SCAN ruins its high accuracy for pure water. Here we show that systematic application of the principles of density-corrected DFT yields a functional (HF-r2SCAN-DC4) which recovers and not only improves over HF-SCAN for pure water, but also captures vital noncovalent interactions in biomolecules, making it suitable for simulations of solutions.
Collapse
Affiliation(s)
- Suhwan Song
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea ,grid.266093.80000 0001 0668 7243Department of Chemistry, University of California, Irvine, CA 92697 USA
| | - Stefan Vuckovic
- grid.472716.10000 0004 1758 7362Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy ,grid.12380.380000 0004 1754 9227Departments of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Youngsam Kim
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea
| | - Hayoung Yu
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722, Korea.
| | - Kieron Burke
- grid.266093.80000 0001 0668 7243Department of Chemistry, University of California, Irvine, CA 92697 USA ,grid.266093.80000 0001 0668 7243Departments of Physics & Astronomy, University of California, Irvine, CA 92697 USA
| |
Collapse
|
7
|
Grupa U, Liebl K, Zacharias M. Orientation Dependence of DNA Blunt-End Stacking Studied by Free-Energy Simulations. J Phys Chem B 2021; 125:13850-13857. [PMID: 34928161 DOI: 10.1021/acs.jpcb.1c07829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
DNA blunt ends can associate mediated by stacking interactions between the terminal base pairs that form blunt ends. The blunt end association plays a role in DNA repair and recombination processes and can also be of importance for the design of DNA-based nano-materials. Its function depends on the sequence and on the geometric arrangement that leads to stable interaction. For a stacked state, the relative orientation (twisting) of the base pairs is important. Molecular dynamics and advanced sampling simulations were used to calculate free energy change associated with twist changes of the stacked blunt-end base pairs. The calculations reproduce blunt stacking arrangements found in crystal structures of DNA oligonucleotides as free energy minima. To elucidate the physical origin of the stabilization of certain angular arrangements, the interactions between backbone atoms in the blunt-end stack were switched off in additional free energy calculations. It allows us to decipher the contributions to stacking stabilization due to the nucleobases and the backbone and to analyze the sequence dependence of the angular stacking preferences. Good qualitative agreement was also found for the comparison with quantum mechanical calculations. The results may help in the design of novel DNA-based materials.
Collapse
Affiliation(s)
- Ulrich Grupa
- Center of Functional Protein Assemblies and Physics Department, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| | - Korbinian Liebl
- Center of Functional Protein Assemblies and Physics Department, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| | - Martin Zacharias
- Center of Functional Protein Assemblies and Physics Department, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| |
Collapse
|
8
|
Mráziková K, Šponer J, Mlýnský V, Auffinger P, Kruse H. Short-Range Imbalances in the AMBER Lennard-Jones Potential for (Deoxy)Ribose···Nucleobase Lone-Pair···π Contacts in Nucleic Acids. J Chem Inf Model 2021; 61:5644-5657. [PMID: 34738826 DOI: 10.1021/acs.jcim.1c01047] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The lone-pair···π (lp···π) (deoxy)ribose···nucleobase stacking is a recurring interaction in Z-DNA and RNAs that is characterized by sub-van der Waals lp···π contacts (<3.0 Å). It is a part of the structural signature of CpG Z-step motifs in Z-DNA and r(UNCG) tetraloops that are known to behave poorly in molecular dynamics (MD) simulations. Although the exact origin of the MD simulation issues remains unclear, a significant part of the problem might be due to an imbalanced description of nonbonded interactions, including the characteristic lp···π stacking. To gain insights into the links between lp···π stacking and MD, we present an in-depth comparison between accurate large-basis-set double-hybrid Kohn-Sham density functional theory calculations DSD-BLYP-D3/ma-def2-QZVPP (DHDF-D3) and data obtained with the nonbonded potential of the AMBER force field (AFF) for NpN Z-steps (N = G, A, C, and U). Among other differences, we found that the AFF overestimates the DHDF-D3 lp···π distances by ∼0.1-0.2 Å, while the deviation between the DHDF-D3 and AFF descriptions sharply increases in the short-range region of the interaction. Based on atom-in-molecule polarizabilities and symmetry-adapted perturbation theory analysis, we inferred that the DHDF-D3 versus AFF differences partly originate in identical nucleobase carbon atom Lennard-Jones (LJ) parameters despite the presence/absence of connected electron-withdrawing groups that lead to different effective volumes or vdW radii. Thus, to precisely model the very short CpG lp···π contact distances, we recommend revision of the nucleobase atom LJ parameters. Additionally, we suggest that the large discrepancy between DHDF-D3 and AFF short-range repulsive part of the interaction energy potential may significantly contribute to the poor performances of MD simulations of nucleic acid systems containing Z-steps. Understanding where, and if possible why, the point-charge-type effective potentials reach their limits is vital for developing next-generation FFs and for addressing specific issues in contemporary MD simulations.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Šlechtitelů 241/27, 783 71 Olomouc-Holice, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Šlechtitelů 241/27, 783 71 Olomouc-Holice, Czech Republic
| | - Pascal Auffinger
- Architecture and Reactivity of RNA, University of Strasbourg, Institute of Molecular and Cellular Biology of the CNRS, 67084 Strasbourg, France
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
9
|
Kaur S, Grover P, Wetmore SD, Sharma P. Role of Stacking Interactions in the Stability of Primitive Genetics: A Quantum Chemical View. J Chem Inf Model 2021; 61:4321-4330. [PMID: 34491053 DOI: 10.1021/acs.jcim.1c00276] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The origin of genetic material on earth is an age-old, entangled mystery that lacks a unanimous explanation. Recent studies have suggested that noncanonical bases such as barbituric acid (BA), melamine (MM), cyanuric acid (CA), and 2,4,6-triaminopyrimidine (TAP) may have undergone molecular selection within the "prebiotic soup" to spontaneously form supramolecular assemblies, which then covalently assembled into an RNA-like polymer (preRNA). However, information on the role of intrinsic interactions of these candidate heterocycles in their molecular selection as the components of preRNA, and the subsequent transition from preRNA to RNA, is currently missing in the literature. To fill this gap in our knowledge on the origin and evolution of primitive genetics, the present work employs density functional theory (B3LYP-D3) to evaluate and compare the stacking propensities of dimers containing prebiotic noncanonical (BA, MM, CA, and TAP) and/or canonical RNA bases (A, C, G, and U). Our detailed analysis of the variation in stacking strength with respect to four characteristic geometrical parameters between the monomers [i.e., the vertical distance, the angle of rotation, and (two) displacements in the x and y directions] reveals that stacking between nonidentical bases is preferred over identical bases for both prebiotic-prebiotic and canonical-canonical dimers. This not only underscores the similarity between the fundamental chemical properties of preRNA and RNA constituents but also supports the likelihood of the evolution of modern (RNA) genetics from primitive (preRNA) genetics. Furthermore, greater average stacking stabilization of canonical dimers than that of dimers containing one canonical and one preRNA nucleobase (by ∼5 kJ mol-1) or dimers solely containing preRNA nucleobases (by ∼12 kJ mol-1) indicates that enhanced stacking is an important factor that may have spurred the evolution of preRNA to an intermediate informational polymer to RNA. More importantly, our study identifies the central roles of CA, BA, and TAP in stacking stabilization within the preRNA and of BA in stacking interactions within the intermediate polymers and suggests that these heterocycles may have played distinct roles in various stages during the evolution from preRNA to RNA. Overall, our results highlight the significance of stacking interactions in the selection of nucleobase components of preRNA.
Collapse
Affiliation(s)
- Sarabjeet Kaur
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Payal Grover
- Department of Chemistry, Dayanand Anglo-Vedic (DAV) College, Sector 10, Chandigarh 16011, India
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| |
Collapse
|
10
|
Gyevi-Nagy L, Kállay M, Nagy PR. Accurate Reduced-Cost CCSD(T) Energies: Parallel Implementation, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2021; 17:860-878. [PMID: 33400527 PMCID: PMC7884001 DOI: 10.1021/acs.jctc.0c01077] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Indexed: 11/28/2022]
Abstract
The accurate and systematically improvable frozen natural orbital (FNO) and natural auxiliary function (NAF) cost-reducing approaches are combined with our recent coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] implementations. Both of the closed- and open-shell FNO-CCSD(T) codes benefit from OpenMP parallelism, completely or partially integral-direct density-fitting algorithms, checkpointing, and hand-optimized, memory- and operation count effective implementations exploiting all permutational symmetries. The closed-shell CCSD(T) code requires negligible disk I/O and network bandwidth, is MPI/OpenMP parallel, and exhibits outstanding peak performance utilization of 50-70% up to hundreds of cores. Conservative FNO and NAF truncation thresholds benchmarked for challenging reaction, atomization, and ionization energies of both closed- and open-shell species are shown to maintain 1 kJ/mol accuracy against canonical CCSD(T) for systems of 31-43 atoms even with large basis sets. The cost reduction of up to an order of magnitude achieved extends the reach of FNO-CCSD(T) to systems of 50-75 atoms (up to 2124 atomic orbitals) with triple- and quadruple-ζ basis sets, which is unprecedented without local approximations. Consequently, a considerably larger portion of the chemical compound space can now be covered by the practically "gold standard" quality FNO-CCSD(T) method using affordable resources and about a week of wall time. Large-scale applications are presented for organocatalytic and transition-metal reactions as well as noncovalent interactions. Possible applications for benchmarking local CCSD(T) methods, as well as for the accuracy assessment or parametrization of less complete models, for example, density functional approximations or machine learning potentials, are also outlined.
Collapse
Affiliation(s)
- László Gyevi-Nagy
- 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
|
11
|
Kruse H, Szabla R, Šponer J. Surprisingly broad applicability of the cc-pVnZ-F12 basis set for ground and excited states. J Chem Phys 2020; 152:214104. [PMID: 32505162 DOI: 10.1063/5.0006871] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
Excellent convergence properties for the (aug-)cc-pVnZ-F12 basis set family, purpose-made for explicitly correlated calculations, are demonstrated with conventional wave function methods and Kohn-Sham density functional theory for various ground and excited-state calculations. Among the ground-state properties studied are dipole moments, covalent bond lengths, and interaction and reaction energies. For excited states, we looked at vertical excitation energies, UV absorption, and excited-state absorption spectra. Convergence is compared against the basis sets cc-pVnZ, def2-nVD, aug-pcseg-n, and nZaPa-NR. It is established that the cc-pVnZ-F12 family consistently yields results of n + 1 quality and better. Especially, the cc-pVDZ-F12 basis set is found to be a basis set of good cost vs performance trade-off.
Collapse
Affiliation(s)
- Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61265 Brno, Czech Republic
| | - Rafał Szabla
- EaStCHEM, School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61265 Brno, Czech Republic
| |
Collapse
|