1
|
Genoni A, Martín Pendás Á. Critical assessment of the x-ray restrained wave function approach: Advantages, drawbacks, and perspectives for density functional theory and periodic ab initio calculations. J Chem Phys 2024; 160:234108. [PMID: 38899684 DOI: 10.1063/5.0213247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Accepted: 06/03/2024] [Indexed: 06/21/2024] Open
Abstract
The x-ray restrained wave function (XRW) method is a quantum crystallographic technique to extract wave functions compatible with experimental x-ray diffraction data. The approach looks for wave functions that minimize the energies of the investigated systems and also reproduce sets of x-ray structure factors. Given the strict relationship between x-ray structure factors and electron distributions, the strategy practically allows determining wave functions that correspond to given (usually experimental) electron densities. In this work, the capabilities of the XRW approach were further tested. The aim was to evaluate whether the XRW technique could serve as a tool for suggesting new exchange-correlation functionals for density functional theory or refining existing ones. Additionally, the ability of the method to address the influences of the crystalline environment was also assessed. The outcomes of XRW computations were thus compared to those of traditional gas-phase, embedding quantum mechanics/molecular mechanics, and fully periodic calculations. The results revealed that, irrespective of the initial conditions, the XRW computations practically yield a consensus electron density, in contrast to the currently employed density functional approximations (DFAs), which tend to give a too large range of electron distributions. This is encouraging in view of exploiting the XRW technique to develop improved functionals. Conversely, the calculations also emphasized that the XRW method is limited in its ability to effectively address the influences of the crystalline environment. This underscores the need for a periodic XRW technique, which would allow further untangling the shortcomings of DFAs from those inherent to the XRW approach.
Collapse
Affiliation(s)
- Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, 57078 Metz, France
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, Avenida Julian Clavería 8, 33006 Oviedo, Spain
| |
Collapse
|
2
|
Macetti G, Genoni A. Introduction of a weighting scheme for the X-ray restrained wavefunction approach: advantages and drawbacks. Acta Crystallogr A Found Adv 2023; 79:25-40. [PMID: 36601761 DOI: 10.1107/s2053273322010221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2022] [Accepted: 10/23/2022] [Indexed: 11/22/2022] Open
Abstract
In a quite recent study [Genoni et al. (2017). IUCrJ, 4, 136-146], it was observed that the X-ray restrained wavefunction (XRW) approach allows a more efficient and larger capture of electron correlation effects on the electron density if high-angle reflections are not considered in the calculations. This is due to the occurrence of two concomitant effects when one uses theoretical X-ray diffraction data corresponding to a single-molecule electron density in a large unit cell: (i) the high-angle reflections are generally much more numerous than the low- and medium-angle ones, and (ii) they are already very well described at unrestrained level. Nevertheless, since high-angle data also contain important information that should not be disregarded, it is not advisable to neglect them completely. For this reason, based on the results of the previous investigation, this work introduces a weighting scheme for XRW calculations to up-weight the contribution of low- and medium-angle reflections, and, at the same time, to reasonably down-weight the importance of the high-angle data. The proposed strategy was tested through XRW computations with both theoretical and experimental structure-factor amplitudes. The tests have shown that the new weighting scheme works optimally if it is applied with theoretically generated X-ray diffraction data, while it is not advantageous when traditional experimental X-ray diffraction data (even of very high resolution) are employed. This also led to the conclusion that the use of a specific external parameter λJ for each resolution range might not be a suitable strategy to adopt in XRW calculations exploiting experimental X-ray data as restraints.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques, 1 Boulevard Arago, Metz, F-57078, France
| | - Alessandro Genoni
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques, 1 Boulevard Arago, Metz, F-57078, France
| |
Collapse
|
3
|
Genoni A. On the termination of the X-ray constrained wavefunction procedure: reformulation of the method for an unequivocal determination of λ. Acta Crystallogr A Found Adv 2022; 78:302-308. [DOI: 10.1107/s2053273322003746] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Accepted: 04/05/2022] [Indexed: 11/10/2022] Open
Abstract
The X-ray constrained/restrained wavefunction (XCW/XRW) approach of quantum crystallography is revisited by introducing the stationary condition of the Jayatilaka functional with respect to the Lagrange multiplier λ. The theoretical derivation has unequivocally shown that the right value of λ is a maximum stationary point of the functional to optimize, thus enabling the solution of the longstanding problem of establishing the point at which to halt the XCW/XRW procedure. Based on the new finding, a reformulation of the X-ray constrained wavefunction algorithm is proposed and its implementation is envisaged. In addition to relying on more solid mathematical grounds, the new variant of the method will be intrinsically more physically meaningful, allowing a straightforward evaluation of the highest level of confidence with which the experimental X-ray diffraction data can be possibly reproduced.
Collapse
|
4
|
Bürgi HB, Genoni A. Remarks on X-ray constrained/restrained wavefunction fitting. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2022; 78:298-304. [PMID: 35695103 DOI: 10.1107/s2052520622004164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2022] [Accepted: 04/19/2022] [Indexed: 06/15/2023]
Abstract
X-ray constrained/restrained wavefunctions (XCWs/XRWs) result from a combination of theory and experiment and are therefore affected by experimental errors and model uncertainties. The present XCW/XRW procedure does not take this into account, thus limiting the meaning and significance of the obtained wavefunctions.
Collapse
Affiliation(s)
- Hans Beat Bürgi
- Department of Chemistry, Biochemistry and Pharmacy, University of Berne, Freiestr. 12, Bern, CH-3012, Switzerland
| | - Alessandro Genoni
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques, 1 Boulevard Arago, Metz, 57050, France
| |
Collapse
|
5
|
Milašinović V, Molčanov K, Krawczuk A, Bogdanov NE, Zakharov BA, Boldyreva EV, Jelsch C, Kojić-Prodić B. Charge density studies of multicentre two-electron bonding of an anion radical at non-ambient temperature and pressure. IUCRJ 2021; 8:644-654. [PMID: 34258012 PMCID: PMC8256703 DOI: 10.1107/s2052252521005273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Accepted: 05/18/2021] [Indexed: 06/13/2023]
Abstract
The variation of charge density of two-electron multicentre bonding (pancake bonding) between semi-quinone radicals with pressure and temperature was studied on a salt of 5,6-di-chloro-2,3-di-cyano-semi-quinone radical anion (DDQ) with 4-cyano-N-methyl-pyridinium cation (4-CN) using the Transferable Aspheric Atom Model (TAAM) refinement. The pancake-bonded radical dimers are stacked by non-bonding π-interactions. With rising pressure, the covalent character of interactions between radicals increases, and above 2.55 GPa, the electron density indicates multicentric covalent interactions throughout the stack. The experimental charge densities were verified and corroborated by periodic DFT computations. The TAAM approach has been tested and validated for atomic resolution data measured at ambient pressure; this work shows this approach can also be applied to diffraction data obtained at pressures up to several gigapascals.
Collapse
Affiliation(s)
- Valentina Milašinović
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 10000, Croatia
| | - Krešimir Molčanov
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 10000, Croatia
| | - Anna Krawczuk
- Institute of Inorganic Chemistry, University of Göttingen, Tammannstrasse 4, Göttingen 37077, Germany
- Faculty of Chemistry, Jagiellonian University in Krakow, Gronostajowa 2, Krakow 30-387, Poland
| | - Nikita E. Bogdanov
- Boreskov Institute of Catalysis, SB RAS, Lavrentiev Avenue 5, Novosibirsk 630090, Russian Federation
- Novosibirsk State University, Pirogova Street 2, Novosibirsk 630090 Russian Federation
| | - Boris A. Zakharov
- Boreskov Institute of Catalysis, SB RAS, Lavrentiev Avenue 5, Novosibirsk 630090, Russian Federation
- Novosibirsk State University, Pirogova Street 2, Novosibirsk 630090 Russian Federation
| | - Elena V. Boldyreva
- Boreskov Institute of Catalysis, SB RAS, Lavrentiev Avenue 5, Novosibirsk 630090, Russian Federation
- Novosibirsk State University, Pirogova Street 2, Novosibirsk 630090 Russian Federation
| | - Christian Jelsch
- CRM2, CNRS, UMR 7036, Université de Lorraine, BP 70239 Nancy, France
| | - Biserka Kojić-Prodić
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 10000, Croatia
| |
Collapse
|
6
|
Podhorský M, Bučinský L, Jayatilaka D, Grabowsky S. HgH 2 meets relativistic quantum crystallography. How to teach relativity to a non-relativistic wavefunction. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2021; 77:54-66. [PMID: 33399131 DOI: 10.1107/s2053273320014837] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Accepted: 11/09/2020] [Indexed: 11/10/2022]
Abstract
The capability of X-ray constrained wavefunction (XCW) fitting to introduce relativistic effects into a non-relativistic wavefunction is tested. It is quantified how much of the reference relativistic effects can be absorbed in the non-relativistic XCW calculation when fitted against relativistic structure factors of a model HgH2 molecule. Scaling of the structure-factor sets to improve the agreement statistics is found to introduce a significant systematic error into the XCW fitting of relativistic effects.
Collapse
Affiliation(s)
- Michal Podhorský
- Institute of Physical Chemistry and Chemical Physics FCHPT, Slovak University of Technology, Radlinskeho 9, Bratislava SK-812 37, Slovakia
| | - Lukáš Bučinský
- Institute of Physical Chemistry and Chemical Physics FCHPT, Slovak University of Technology, Radlinskeho 9, Bratislava SK-812 37, Slovakia
| | - Dylan Jayatilaka
- School of Molecular Sciences, The University of Western Australia, 35 Stirling Hwy, Perth WA 6009, Australia
| | - Simon Grabowsky
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland
| |
Collapse
|
7
|
Affiliation(s)
- Piero Macchi
- Department, Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Milano, Italy
- Center for Nano Science and Technology CNST@polimi, Italian Institute of Technology, Milano, Italy
| |
Collapse
|
8
|
Abstract
In this review article, we report on the recent progresses in the field of quantum crystallography that has witnessed a massive increase of production coupled with a broadening of the scope in the last decade. It is shown that the early thoughts about extracting quantum mechanical information from crystallographic experiments are becoming reality, although a century after prediction. While in the past the focus was mainly on electron density and related quantities, the attention is now shifting toward determination of wavefunction from experiments, which enables an exhaustive determination of the quantum mechanical functions and properties of a system. Nonetheless, methods based on electron density modelling have evolved and are nowadays able to reconstruct tiny polarizations of core electrons, coupling charge and spin models, or determining the quantum behaviour at extreme conditions. Far from being routine, these experimental and computational results should be regarded with special attention by scientists for the wealth of information on a system that they actually contain.
Collapse
|
9
|
Post-Hartree-Fock methods for Hirshfeld atom refinement: are they necessary? Investigation of a strongly hydrogen-bonded molecular crystal. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.127934] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
10
|
Ernst M, Genoni A, Macchi P. Analysis of crystal field effects and interactions using X-ray restrained ELMOs. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.127975] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
11
|
Macetti G, Wieduwilt EK, Assfeld X, Genoni A. Localized Molecular Orbital-Based Embedding Scheme for Correlated Methods. J Chem Theory Comput 2020; 16:3578-3596. [PMID: 32369363 DOI: 10.1021/acs.jctc.0c00084] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Embedding strategies currently provide the best compromise between accuracy and computational cost in modeling chemical properties and processes of large and complex systems. In this framework, different methods have been proposed all over the years, from the very popular QM/MM approaches to the more recent and very promising density matrix and density functional embedding techniques. Here, we present a further development of the quantum mechanics/extremely localized molecular orbital technique (QM/ELMO) method, a recently proposed multiscale embedding strategy in which the chemically active region of the investigated system is treated at a fully quantum mechanical level, while the rest is described by frozen extremely localized molecular orbitals previously transferred from proper libraries or tailor-made model molecules. In particular, in this work we discuss and assess in detail the extension of the QM/ELMO approach to density functional theory and post-Hartree-Fock techniques by evaluating its performances when it is used to describe chemical reactions, bond dissociations, and intermolecular interactions. The preliminary test calculations have shown that, in the investigated cases, the new embedding strategy enables the results of the corresponding fully quantum mechanical computations to be reproduced within chemical accuracy in almost all the cases but with a significantly reduced computational cost, especially when correlated post-Hartree-Fock strategies are used to describe the quantum mechanical subsystem. In light of the obtained results, we already envisage the future application of the new correlated QM/ELMO techniques to the investigation of more challenging problems, such as the modeling of enzyme catalysis, the study of excited states of biomolecules, and the refinement of macromolecular X-ray crystal structures.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Erna K Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Xavier Assfeld
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, Boulevard des Aiguilletes, BP 70239, F-54506 Vandoeuvre-lès-Nancy, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
12
|
Franchini D, Forni A, Genoni A, Pieraccini S, Gandini E, Sironi M. The Origin of the σ-Hole in Halogen Atoms: a Valence Bond Perspective. ChemistryOpen 2020; 9:445-450. [PMID: 32269900 PMCID: PMC7136650 DOI: 10.1002/open.202000062] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Indexed: 11/26/2022] Open
Abstract
A detailed Valence Bond-Spin Coupled analysis of a series of halogenated molecules is here reported, allowing to get a rigorous ab initio demonstration of the qualitative models previously proposed to explain the origin of halogen bonding. The concepts of σ-hole and negative belt observed around the halogen atoms in the electrostatic potential maps are here interpreted by analysis of the relevant Spin Coupled orbitals.
Collapse
Affiliation(s)
- Davide Franchini
- Department of ChemistryUniversità degli Studi di Milano and INSTM RUvia Golgi 1920133MilanoItaly
| | | | - Alessandro Genoni
- Université de Lorraine and CNRSLaboratoire de Physique et Chimie Théoriques (LPCT)UMR CNRS 70191 Boulevard AragoMetz57078France
| | - Stefano Pieraccini
- Department of ChemistryUniversità degli Studi di Milano and INSTM RUvia Golgi 1920133MilanoItaly
- SCITEC-CNRvia Golgi 1920133MilanoItaly
| | - Enrico Gandini
- Department of ChemistryUniversità degli Studi di Milano and INSTM RUvia Golgi 1920133MilanoItaly
| | - Maurizio Sironi
- Department of ChemistryUniversità degli Studi di Milano and INSTM RUvia Golgi 1920133MilanoItaly
- SCITEC-CNRvia Golgi 1920133MilanoItaly
| |
Collapse
|
13
|
Genoni A. On the use of the Obara–Saika recurrence relations for the calculation of structure factors in quantum crystallography. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2020; 76:172-179. [DOI: 10.1107/s205327332000042x] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/07/2019] [Accepted: 01/14/2020] [Indexed: 11/10/2022]
Abstract
Modern methods of quantum crystallography are techniques firmly rooted in quantum chemistry and, as in many quantum chemical strategies, electron densities are expressed as two-centre expansions that involve basis functions centred on atomic nuclei. Therefore, the computation of the necessary structure factors requires the evaluation of Fourier transform integrals of basis function products. Since these functions are usually Cartesian Gaussians, in this communication it is shown that the Fourier integrals can be efficiently calculated by exploiting an extension of the Obara–Saika recurrence formulas, which are successfully used by quantum chemists in the computation of molecular integrals. Implementation and future perspectives of the technique are also discussed.
Collapse
|
14
|
Grabowsky S, Genoni A, Thomas SP, Jayatilaka D. The Advent of Quantum Crystallography: Form and Structure Factors from Quantum Mechanics for Advanced Structure Refinement and Wavefunction Fitting. 21ST CENTURY CHALLENGES IN CHEMICAL CRYSTALLOGRAPHY II 2020. [DOI: 10.1007/430_2020_62] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
15
|
Genoni A, Macetti G, Franchini D, Pieraccini S, Sironi M. X-ray constrained spin-coupled technique: theoretical details and further assessment of the method. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2019; 75:778-797. [PMID: 31692454 DOI: 10.1107/s2053273319011021] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/24/2019] [Accepted: 08/07/2019] [Indexed: 11/11/2022]
Abstract
One of the well-established methods of modern quantum crystallography is undoubtedly the X-ray constrained wavefunction (XCW) approach, a technique that enables the determination of wavefunctions which not only minimize the energy of the system under examination, but also reproduce experimental X-ray diffraction data within the limit of the experimental errors. Initially proposed in the framework of the Hartree–Fock method, the strategy has been gradually extended to other techniques of quantum chemistry, but always remaining limited to a single-determinant ansatz for the wavefunction to extract. This limitation has been recently overcome through the development of the novel X-ray constrained spin-coupled (XCSC) approach [Genoni et al. (2018). Chem. Eur. J.
24, 15507–15511] which merges the XCW philosophy with the traditional spin-coupled strategy of valence bond theory. The main advantage of this new technique is the possibility of extracting traditional chemical descriptors (e.g. resonance structure weights) compatible with the experimental diffraction measurements, without the need to introduce information a priori or perform analyses a posteriori. This paper provides a detailed theoretical derivation of the fundamental equations at the basis of the XCSC method and also introduces a further advancement of its original version, mainly consisting in the use of molecular orbitals resulting from XCW calculations at the Hartree–Fock level to describe the inactive electrons in the XCSC computations. Furthermore, extensive test calculations, which have been performed by exploiting high-resolution X-ray diffraction data for salicylic acid and by adopting different basis sets, are presented and discussed. The computational tests have shown that the new technique does not suffer from particular convergence problems. Moreover, all the XCSC calculations provided resonance structure weights, spin-coupled orbitals and global electron densities slightly different from those resulting from the corresponding unconstrained computations. These discrepancies can be ascribed to the capability of the novel strategy to capture the information intrinsically contained in the experimental data used as external constraints.
Collapse
|
16
|
Meyer B, Genoni A. Libraries of Extremely Localized Molecular Orbitals. 3. Construction and Preliminary Assessment of the New Databanks. J Phys Chem A 2018; 122:8965-8981. [PMID: 30339393 DOI: 10.1021/acs.jpca.8b09056] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The fast and reliable determination of wave functions and electron densities of macromolecules has been one of the goals of theoretical chemistry for a long time, and in this context, several linear scaling techniques have been successfully devised over the years. Different approaches have been adopted to tackle this problem, and one of them exploits the fact that, according to the traditional chemical perception, molecules can be seen as constituted of recurring units (e.g., functional groups) with well-defined chemical features. This has led to the development of methods in which the global wave functions or electron densities of macromolecules are obtained by simply transferring density matrices or fuzzy electron densities associated with molecular fragments. In this context, we propose an alternative strategy that aims at quickly reconstructing wave functions and electron densities of proteins through the transfer of extremely localized molecular orbitals (ELMOs), which are orbitals strictly localized on small molecular units and, for this reason, easily transferable from molecule to molecule. To accomplish this task we have constructed original libraries of ELMOs that cover all the possible elementary fragments of the 20 natural amino acids in all their possible protonation states and forms. Our preliminary test calculations have shown that, compared to more traditional methods of quantum chemistry, the transfers from the novel ELMO databanks allow to obtain wave function and electron densities of large polypeptides and proteins at a significantly reduced computational cost. Furthermore, notwithstanding expected discrepancies, the obtained electron distributions and electrostatic potentials are in very good agreement with those obtained at Hartree-Fock and density functional theory (DFT) levels. Therefore, the results encourage to use the new libraries as alternatives to the popular pseudoatom-databases of crystallography in the refinement of crystallographic structures of macromolecules. In particular, in this context, we have already envisaged the coupling of the ELMO databanks with the promising Hirshfeld atom refinement technique to extend the applicability of the latter to very large systems.
Collapse
Affiliation(s)
- Benjamin Meyer
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019 , 1 Boulevard Arago , F-57078 Metz , France
| | - Alessandro Genoni
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019 , 1 Boulevard Arago , F-57078 Metz , France
| |
Collapse
|