1
|
Zheng J, Frisch MJ. Multiple-time scale integration method based on an interpolated potential energy surface for ab initio path integral molecular dynamics. J Chem Phys 2024; 160:144111. [PMID: 38597307 DOI: 10.1063/5.0196634] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 03/25/2024] [Indexed: 04/11/2024] Open
Abstract
A new multiple-time scale integration method is presented that propagates ab initio path integral molecular dynamics (PIMD). This method uses a large time step to generate an approximate geometrical configuration whose energy and gradient are evaluated at the level of an ab initio method, and then, a more precise integration scheme, e.g., the Bulirsch-Stoer method or velocity Verlet integration with a smaller time step, is used to integrate from the previous step using the computationally efficient interpolated potential energy surface constructed from two consecutive points. This method makes the integration of PIMD more efficient and accurate compared with the velocity Verlet integration. A Nosé-Hoover chain thermostat combined with this new multiple-time scale method has good energy conservation even with a large time step, which is usually challenging in velocity Verlet integration for PIMD due to the very small chain mass when a large number of beads are used. The new method is used to calculate infrared spectra and free energy profiles to demonstrate its accuracy and capabilities.
Collapse
Affiliation(s)
- Jingjing Zheng
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Michael J Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| |
Collapse
|
2
|
Christensen EG, Lutz KT, Spencer RJ, Steele RP. Persistence of a Delocalized Radical in Larger Clusters of Hydrated Copper(II) Hydroxide, CuOH +(H 2O) 3-7. J Phys Chem A 2023; 127:6647-6659. [PMID: 37587877 DOI: 10.1021/acs.jpca.3c03718] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2023]
Abstract
The structures, vibrational spectra, and electronic properties of copper hydroxide hydrates CuOH+(H2O)3-7 were investigated with quantum chemistry computations. As a follow-up to a previous analysis of CuOH+(H2O)0-2, this investigation examined the progression as the square-planar metal coordination environment was filled and as solvation shells expanded. Four-, five-, and six-coordinate structures were found to be low-energy isomers. The delocalized radical character, which was discovered in the small clusters, was found to persist upon continued hydration, although the hydrogen-bonded water network in the larger clusters was found to play a more significant role in accommodating this spin. Partial charges indicated that the electronic structure includes more Cu2+···OH- character than was observed in smaller clusters, but this structure remains decidedly mixed with Cu+···OH· configurations and yields roughly half-oxidation of the water network in the absence of any electrochemical potential. Computed vibrational spectra for n = 3 showed congruence with spectra from recent predissociation spectroscopy experiments, provided that the role of the D2 tag was taken into account. Spectra for n = 4-7 were predicted to exhibit features that are reflective of both the mixed electronic character and proton-/hydrogen-shuttling motifs within the hydrogen-bonded water network.
Collapse
Affiliation(s)
- Elizabeth G Christensen
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Kevin T Lutz
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Ryan J Spencer
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| |
Collapse
|
3
|
Schiltz C, Rappoport D, Mandelshtam VA. Implementation of the self-consistent phonons method with ab initio potentials (AI-SCP). J Chem Phys 2023; 158:2890485. [PMID: 37184023 DOI: 10.1063/5.0146682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/24/2023] [Indexed: 05/16/2023] Open
Abstract
The self-consistent phonon (SCP) method allows one to include anharmonic effects when treating a many-body quantum system at thermal equilibrium. The system is then described by an effective temperature-dependent harmonic Hamiltonian, which can be used to estimate its various dynamic and static properties. In this paper, we combine SCP with ab initio (AI) potential energy evaluation in which case the numerical bottleneck of AI-SCP is the evaluation of Gaussian averages of the AI potential energy and its derivatives. These averages are computed efficiently by the quasi-Monte Carlo method utilizing low-discrepancy sequences leading to a fast convergence with respect to the number, S, of the AI energy evaluations. Moreover, a further substantial (an-order-of-magnitude) improvement in efficiency is achieved once a numerically cheap approximation of the AI potential is available. This is based on using a perturbation theory-like (the two-grid) approach in which it is the average of the difference between the AI and the approximate potential that is computed. The corresponding codes and scripts are provided.
Collapse
Affiliation(s)
- Colin Schiltz
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | | |
Collapse
|
4
|
Mauger N, Plé T, Lagardère L, Huppert S, Piquemal JP. Improving Condensed-Phase Water Dynamics with Explicit Nuclear Quantum Effects: The Polarizable Q-AMOEBA Force Field. J Phys Chem B 2022; 126:8813-8826. [PMID: 36270033 DOI: 10.1021/acs.jpcb.2c04454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We introduce a new parametrization of the AMOEBA polarizable force field for water denoted Q-AMOEBA, for use in simulations that explicitly account for nuclear quantum effects (NQEs). This study is made possible thanks to the recently introduced adaptive Quantum Thermal Bath (adQTB) simulation technique which computational cost is comparable to classical molecular dynamics. The flexible Q-AMOEBA model conserves the initial AMOEBA functional form, with an intermolecular potential including an atomic multipole description of electrostatic interactions (up to quadrupole), a polarization contribution based on the Thole interaction model and a buffered 14-7 potential to model van der Waals interactions. It has been obtained by using a ForceBalance fitting strategy including high-level quantum chemistry reference energies and selected condensed-phase properties targets. The final Q-AMOEBA model is shown to accurately reproduce both gas-phase and condensed-phase properties, notably improving the original AMOEBA water model. This development allows the fine study of NQEs on water liquid phase properties such as the average H-O-H angle compared to its gas-phase equilibrium value, isotope effects, and so on. Q-AMOEBA also provides improved infrared spectroscopy prediction capabilities compared to AMOEBA03. Overall, we show that the impact of NQEs depends on the underlying model functional form and on the associated strength of hydrogen bonds. Since adQTB simulations can be performed at near classical computational cost using the Tinker-HP package, Q-AMOEBA can be extended to organic molecules, proteins, and nucleic acids opening the possibility for the large-scale study of the importance of NQEs in biophysics.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Thomas Plé
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| |
Collapse
|
5
|
Li C, Voth GA. Using Machine Learning to Greatly Accelerate Path Integral Ab Initio Molecular Dynamics. J Chem Theory Comput 2022; 18:599-604. [PMID: 34982562 PMCID: PMC8864787 DOI: 10.1021/acs.jctc.1c01085] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. The situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.
Collapse
|
6
|
Xue Y, Wang JN, Hu W, Zheng J, Li Y, Pan X, Mo Y, Shao Y, Wang L, Mei Y. Affordable Ab Initio Path Integral for Thermodynamic Properties via Molecular Dynamics Simulations Using Semiempirical Reference Potential. J Phys Chem A 2021; 125:10677-10685. [PMID: 34894680 PMCID: PMC9108008 DOI: 10.1021/acs.jpca.1c07727] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Path integral molecular dynamics (PIMD) is becoming a routinely applied method for incorporating the nuclear quantum effect in computer simulations. However, direct PIMD simulations at an ab initio level of theory are formidably expensive. Using the protonated 1,8-bis(dimethylamino)naphthalene molecule as an example, we show in this work that the computational expense for the intramolecular proton transfer between the two nitrogen atoms can be remarkably reduced by implementing the idea of reference-potential methods. The simulation time can be easily extended to a scale of nanoseconds while maintaining the accuracy on an ab initio level of theory for thermodynamic properties. In addition, postprocessing can be carried out in parallel on massive computer nodes. A 545-fold reduction in the total CPU time can be achieved in this way as compared to a direct PIMD simulation at the same ab initio level of theory.
Collapse
Affiliation(s)
- Yuanfei Xue
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Jia-Ning Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Wenxin Hu
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structure, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China,NYU–ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China,NYU–ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
7
|
Mauger N, Plé T, Lagardère L, Bonella S, Mangaud É, Piquemal JP, Huppert S. Nuclear Quantum Effects in Liquid Water at Near Classical Computational Cost Using the Adaptive Quantum Thermal Bath. J Phys Chem Lett 2021; 12:8285-8291. [PMID: 34427440 DOI: 10.1021/acs.jpclett.1c01722] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQEs) in molecular simulations: the adaptive quantum thermal bath (adQTB). In this method, zero-point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation-dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to that of classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Thomas Plé
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Sara Bonella
- CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne, Batochimie, Avenue Forel 2, 1015 Lausanne, Switzerland
| | - Étienne Mangaud
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
- Institut Universitaire de France, 75005 Paris, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Simon Huppert
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| |
Collapse
|
8
|
Brieuc F, Schran C, Uhl F, Forbert H, Marx D. Converged quantum simulations of reactive solutes in superfluid helium: The Bochum perspective. J Chem Phys 2020; 152:210901. [DOI: 10.1063/5.0008309] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Fabien Brieuc
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Felix Uhl
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
9
|
Elton DC, Fritz M, Fernández-Serra M. Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost. Phys Chem Chem Phys 2018; 21:409-417. [PMID: 30534683 DOI: 10.1039/c8cp06077k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge-Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call "monomer PIMD", correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.
Collapse
Affiliation(s)
- Daniel C Elton
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA.
| | | | | |
Collapse
|
10
|
Schran C, Brieuc F, Marx D. Converged Colored Noise Path Integral Molecular Dynamics Study of the Zundel Cation Down to Ultralow Temperatures at Coupled Cluster Accuracy. J Chem Theory Comput 2018; 14:5068-5078. [PMID: 30217111 DOI: 10.1021/acs.jctc.8b00705] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
For a long time, performing converged path integral simulations at ultralow but finite temperatures of a few Kelvin has been a nearly impossible task. However, recent developments in advanced colored noise thermostatting schemes for path integral simulations, namely, the Path Integral Generalized Langevin Equation Thermostat (PIGLET) and the Path Integral Quantum Thermal Bath (PIQTB), have been able to greatly reduce the computational cost of these simulations, thus making the ultralow temperature regime accessible in practice. In this work, we investigate the influence of these two thermostatting schemes on the description of hydrogen-bonded systems at temperatures down to a few Kelvin as encountered, for example, in helium nanodroplet isolation or tagging photodissociation spectroscopy experiments. For this purpose, we analyze the prototypical hydrogen bond in the Zundel cation (H5O2+) as a function of both oxygen-oxygen distance and temperature in order to elucidate how the anisotropic quantum delocalization and, thus, the shape of the shared proton adapts depending on the donor-acceptor distance. The underlying electronic structure of the Zundel cation is described in terms of Behler's Neural Network Potentials of essentially converged Coupled Cluster accuracy, CCSD(T*)-F12a/AVTZ. In addition, the performances of the PIQTB and PIGLET methods for energetic, structural, and quantum delocalization properties are assessed and directly compared. Overall, our results emphasize the validity and practical usefulness of these two modern thermostatting approaches for path integral simulations of hydrogen-bonded systems even at ultralow temperatures.
Collapse
Affiliation(s)
- Christoph Schran
- Lehrstuhl für Theoretische Chemie , Ruhr-Universität Bochum , 44780 Bochum , Germany
| | - Fabien Brieuc
- Lehrstuhl für Theoretische Chemie , Ruhr-Universität Bochum , 44780 Bochum , Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie , Ruhr-Universität Bochum , 44780 Bochum , Germany
| |
Collapse
|
11
|
Stock G, Hamm P. A non-equilibrium approach to allosteric communication. Philos Trans R Soc Lond B Biol Sci 2018; 373:20170187. [PMID: 29735740 PMCID: PMC5941181 DOI: 10.1098/rstb.2017.0187] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/16/2018] [Indexed: 12/16/2022] Open
Abstract
While the theory of protein folding is well developed, including concepts such as rugged energy landscape, folding funnel, etc., the same degree of understanding has not been reached for the description of the dynamics of allosteric transitions in proteins. This is not only due to the small size of the structural change upon ligand binding to an allosteric site, but also due to challenges in designing experiments that directly observe such an allosteric transition. On the basis of recent pump-probe-type experiments (Buchli et al. 2013 Proc. Natl Acad. Sci. USA110, 11 725-11 730. (doi:10.1073/pnas.1306323110)) and non-equilibrium molecular dynamics simulations (Buchenberg et al. 2017 Proc. Natl Acad. Sci. USA114, E6804-E6811. (doi:10.1073/pnas.1707694114)) studying an photoswitchable PDZ2 domain as model for an allosteric transition, we outline in this perspective how such a description of allosteric communication might look. That is, calculating the dynamical content of both experiment and simulation (which agree remarkably well with each other), we find that allosteric communication shares some properties with downhill folding, except that it is an 'order-order' transition. Discussing the multiscale and hierarchical features of the dynamics, the validity of linear response theory as well as the meaning of 'allosteric pathways', we conclude that non-equilibrium experiments and simulations are a promising way to study dynamical aspects of allostery.This article is part of a discussion meeting issue 'Allostery and molecular machines'.
Collapse
Affiliation(s)
- Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, Freiburg, Germany
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| |
Collapse
|
12
|
|
13
|
Mouhat F, Sorella S, Vuilleumier R, Saitta AM, Casula M. Fully Quantum Description of the Zundel Ion: Combining Variational Quantum Monte Carlo with Path Integral Langevin Dynamics. J Chem Theory Comput 2017; 13:2400-2417. [DOI: 10.1021/acs.jctc.7b00017] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Félix Mouhat
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Sandro Sorella
- International
School for Advanced Studies (SISSA), Via Bonomea 26, 34136 Trieste, Italy
- INFM Democritos
National Simulation Center, 34151 Trieste, Italy
| | - Rodolphe Vuilleumier
- PASTEUR,
Département de chimie, École normale supérieure,
UPMC Univ Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ Paris 06, École normale supérieure, CNRS, Processus d’activation
sélective par transfert d’énergie uni-électronique
ou radiatif (PASTEUR), 75005 Paris, France
| | - Antonino Marco Saitta
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Michele Casula
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| |
Collapse
|