Next Article in Journal
Exploring Bioactivities and Peptide Content of Body Mucus from the Lusitanian Toadfish Halobatrachus didactylus
Next Article in Special Issue
Thermal Conductivity Calculation in Organic Liquids: Application to Poly-α-Olefin
Previous Article in Journal
The Preparation of a Challenging Superconductor Nb3Al by Exploiting Nano Effect
Previous Article in Special Issue
Exploring Routes to Enhance the Calculation of Free Energy Differences via Non-Equilibrium Work SQM/MM Switching Simulations Using Hybrid Charge Intermediates between MM and SQM Levels of Theory or Non-Linear Switching Schemes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nuclear Motion Is Classical: Spectrum of a Magic Protonated Water Cluster

Theoretical Chemistry, Leibniz University Hannover, Callinstr. 3A, 30167 Hannover, Germany
Molecules 2023, 28(18), 6454; https://doi.org/10.3390/molecules28186454
Submission received: 14 August 2023 / Revised: 1 September 2023 / Accepted: 4 September 2023 / Published: 6 September 2023
(This article belongs to the Special Issue Advances in Computational and Theoretical Chemistry)

Abstract

:
The assumption that nuclear motion is classical explains many phenomena. The problems of Schrödinger’s cat and the EPR paradoxon do not exist in a perfectly deterministic theory. All it needs is to describe nuclear motion classically right from the beginning. To establish this simple idea, it must be tested for as many examples as possible. In the present paper, we use ab initio molecular dynamics to investigate the infrared spectrum of a ‘magic’ protonated water cluster H 3 O + (H 2 O) 20 which exhibits some features that were believed to afford a quantum treatment of nuclear motion. The role of the temperature in contrast to a quantum mechanical description is discussed.

Graphical Abstract

1. Introduction

It is perfectly clear that all kind of matter is quantized in the sense that at normal energies everything consists of atoms or, more precisely, of tiny nuclei and an extended electronic wave function. We got used to this picture not only by theoretical considerations, but, thanks to modern microscopy and spectroscopy, also by direct experimental observations of single molecules. Nevertheless, it is not so clear if quantization means that absolutely everything should be described with the Schrödinger equation as proposed by the Copenhagen interpretation. In fact, such a comprehensive use of the Schrödinger equation leads to the well-known paradoxa. These paradoxa are at odds with observation, which is most obvious if one has a look at chemical reactions. We are able to simulate such reactions and to observe how every atom gets from the educt state to the product state without needing any measurement process. The approach we employ is ab initio molecular dynamics (AIMD) using the Car–Parrinello scheme (CPMD) [1,2]. Molecular dynamics (MD) means that we are using Newton dynamics for the nuclei. Ab initio molecular dynamics means that for the description of the electronic cloud the Schrödinger equation is used in a molecular dynamics run. This approach is local and deterministic. Simulating the same reaction two times leads to exactly the same result. Sligthly changing the initial conditions, however, may lead to a completely different product state. We have electronic tunneling and classical chaos for nuclear motion.
The derivation of CPMD can be put in a way as to make no use of the Born–Oppenheimer approximation in order to separate nuclear and electronic motion. The Born–Oppenheimer approximation attributes a wave function to the nuclei and a wave function to the electrons and describes the total wave function as the product of these two wave functions. This would make sense if the problem was separable into an electronic and a nuclear problem. It turns out, it is not. The electronic–nuclear interaction is not negligible and the nuclear–nuclear interaction cannot be omitted if one wants to obtain potential energy surfaces. See, for example, the code in [3].
A more direct derivation starts from assuming the nuclei as classical particles. Note that this does not mean that the nuclei have no structure and no properties at all. An inner structure of the nuclei just does not matter for investigating nuclear motion like it appears in chemical reactions. There is no reason to assume that the nuclear structure obeys the same laws as the electronic structure does. The nuclei are smaller than the lighter electrons, in agreement with the deBroglie relation. However, the heavier the nuclei get in the periodic table, the larger they become while the density stays about the same. A naive application of the deBroglie picture would lead to a clear deviation from the experiment. The nuclei are clearly different from the electrons, also in the sense that they can be discriminated. For the motion in an electrical field generated by the electrons and all other nuclei, we need only their masses and charges.
Matter of any kind consists of atoms, which again consist of electrons and nuclei. The total Hamiltonian for a system with N electrons and M nuclei is:
H ^ t o t a l = T ^ n u c + T ^ e l e c + V ^ n u c n u c + V ^ e l e c n u c + V ^ e l e c e l e c
T ^ n u c is the kinetic energy of the nuclei, T ^ e l e c is the kinetic energy of the electrons, V ^ n u c n u c is the nuclear-nuclear interaction, V ^ e l e c n u c is the electronic–nuclear interaction, and V ^ e l e c e l e c is the electronic–electronic interaction. A normal quantum mechanical calculation uses, for example, the density functional theory (DFT) approximation in the Kohn–Sham (KS) formalism for the description of the electronic structure. In the following, the terms ‘density functional theory’ and ‘Kohn–Sham’ will be used synonymously. In such a DFT calculation, the kinetic energy of the nuclei is put to zero. Since the kinetic energy of the nuclei corresponds to the temperature of the system, this means modelling the situation at zero Kelvin. We obtain H ^ Q C which contains, besides all electronic terms, the nuclear–nuclear interaction V n u c n u c . This term is missing in the traditional derivation according to Born and Oppenheimer as a consequence of separating electronic and nuclear wavefunctions.
H ^ Q C = T ^ e l e c + V ^ n u c n u c + V ^ e l e c n u c + V ^ e l e c e l e c
A separation would lead to the electronic Hamiltonian H ^ e l .
H ^ e l = T ^ e l e c + V ^ e l e c n u c + V ^ e l e c e l e c
H ^ e l , however, does not describe proper Born–Oppenheimer surfaces, in contrast to H ^ Q C . Nevertheless, it is sufficient for deriving equations of motions, because V n u c n u c yields zero when varying the energy expression with respect to the wavefunction.
The derivation of the Car–Parrinello equations starts with an extended Lagrangian using E Q C as obtained from H ^ Q C (capital letters are used for nuclear quantities):
L C P = I 1 2 M I R ˙ ̲ I 2 + μ i ψ ˙ i ψ ˙ i d τ E Q C + i j Λ i j ( ψ i ψ j d τ δ i j )
The four components of the Car–Parrinello Lagrangian are as follows: I 1 2 M I R ˙ ̲ I 2 : kinetic energy of the classical nuclei, μ i ψ ˙ i ψ ˙ i d τ : fictitious kinetic energy of the electrons, E Q C : energy obtained from a quantum chemical calculation of the electronic structure using H ^ Q C , i j Λ i j ( ψ i ψ j d τ δ i j ) : constraints introduced to guarantee orbital orthonormality during the simulation.
Hereby is:
E Q C = E K S + V n u c n u c
E K S is the Kohn–Sham energy computed for the electronic Hamiltonian H ^ e l in density functional approximation, containing the kinetic energy of the electrons, the interaction of the electrons with the nuclei and the electron–electron interaction. E Q C is the energy computed in a normal quantum chemical calculation, here: density functional calculation. E Q C contains all relevant terms except the kinetic energy of the nuclei. It is obtained without making use of a nuclear wave function. The kinetic energy of the nuclei is added in the Car–Parrinello Lagrangian (first term in the Lagrangian) under the assumption of classical particles. Like this, non-zero-temperature situations are modelled. Note that, by adding this classical version of T ^ n u c , all five terms of the energy expression H ^ t o t a l we started from are considered.
The corresponding Euler–Lagrange equations read:
N u c l e i : d d t L R ˙ ̲ I = L R ̲ I
E l e c t r o n s : d d t δ L δ ψ ˙ i = δ L δ ψ i
This results in the following Car–Parrinello equations of motion:
N u c l e i : M I R ¨ ̲ I = R ̲ I E Q C
E l e c t r o n s : μ i ψ ¨ i = [ 1 2 2 + v e f f ( r ̲ ) ] ψ i + j Λ i j ψ j
In deriving the latter equation, the Kohn–Sham approximation was used. It can be read in the following way: The deviation from the Kohn–Sham equations as a consequence of nuclear motion (right side) leads to a force on the orbitals (left side of the equation). Note that the Car–Parrinello equations contain second derivatives both for space and time. While the kinetic energy of the orbitals was introduced ad hoc, it comes as no surprise that something was missing in the time-independent description of electronic structure using the time-independent Schrödinger equation. The well-working Car–Parrinello equations may be similar to a more fundamental equation which is not yet known. It is remarkable that the CPMD equations which contain neither imaginary parts nor single derivatives with respect to time can be solved in a straight-forward way and yield a convincing picture.
The alternative to Car–Parrinello molecular dynamics is so-called Born–Oppenheimer dynamics (BOMD). In this variant of ab initio molecular dynamics, no forces on the orbitals are computed. Instead, the electronic wave function is fully optimized to the ’Born–Oppenheimer’ surface in every time step. It is the more widely used variant, because the implementation is much easier than that of the more efficient Car–Parrinello molecular dynamics. The implementation of CPMD is particularly difficult if Gaussian basis sets are used instead of plane-wave basis sets, due to an additional term taking into account the motion of atom-centered basis functions through space during a dynamics run. For the present consideration, it does not matter if CPMD or BOMD is used, because the treatment of nuclear motion is the same in both approaches.
With the simple idea that nuclear motion is classical, we have arrived at a local and deterministic theory. While it is easy to formulate this new view of quantum chemistry, it is not so easy to show that it is in agreement with experiment. Many phenomena have been attributed to a quantum nuclear motion. This has a long history [4,5,6], but there are also important developments which have been published recently [7,8,9]. It must be shown one by one that such situations can also be explained by nuclear motion at finite temperatures. See our previous publications [10,11,12,13,14,15]. Infrared (IR) spectra are among the most important phenomena that may be attributed to quantum nuclear motion. In the present study we employ CPMD to compute spectra of a water cluster H 3 O + (H 2 O) 20 [16,17] that was found to be particularly stable (Figure 1). Its characteristic IR spectrum contains some features that were attributed to quantum motion of the nuclei. We will show that this assumption is not needed and that temperature is a sufficient concept to explain the observations.

2. Results and Discussions

We performed a total of 27 CPMD simulation runs (see Table 1 and the Section 4 below).
We find an Eigen structure, which means that we have in essence H 3 O + instead of H 5 O 2 + . This is in agreement with previous calculations [17]. See also [18]. Zundel structures, H 5 O 2 + , are observed as transition states when H 3 O + exchanges a proton with a neighbouring water molecule, which happens several times in the simulation runs.
We started with nine simulations that tested the effect of the super-cell size as well as the addition of a chlorine atom as counter ion. (If there is no counter ion, CPMD uses a background charge to describe the periodically repeated charged system). Each of these scenarios was computed for plane-wave cutoffs of 50 and 60 Rydberg. In part of the simulations, a dispersion correction was added. See Table 1. The result is described in Figure 2. Three phenomena are most prominent: the H 2 O libration below 1000 cm 1 , the water peak at about 1600 cm 1 , and the OH stretch vibration at 3000 cm 1 and above. All three scenarios (plane-wave cutoff of 50 and 60 Rydberg and dispersion correction), do not change the picture strongly. An exception is the OH stretch vibration. Here, increasing the cutoff and introducing a dispersion correction brings the peaks farther away from experiment. This is attributed to shortcomings of the BLYP functional. It should be emphasized, however, that overall BLYP is doing a great job, particularly if one considers the cost/value ratio.
The super cell size and the addition of a counter ion have little effect. In other words, the cell is big enough and the background charge applied to the charged system does not influence the result.
Most important is the comparison of the weak H 3 O + peaks to experiment, because these peaks are believed to be influenced by quantum mechanical nuclear motion [17]. With our classical approach, we observe weak features close to 1200 cm 1 and around 2000 cm 1 . However, this observation is a bit accidental. In some of the simulations, some of these peaks are missing. There is no clear trend, hence we assume that this is due to the accidental observation of a particular vibration on the picosecond timescale. It is simply unlikely, that all relevant vibrations are excited during the short simulation time span in a single cluster.
In a second suite of nine simulations, we tested the effect of temperature (Table 1 and Figure 3).
Additionally, here, most spectra are similar, whereby the spectra computed close to 0 K show the largest deviations. In these simulations, the water peak does not appear at all. In contrast, the weak features close to 1200 cm 1 and 2000 cm 1 are obtained correctly. A further exception represents one of the runs at 400 K (simulation 17): At this high temperature, the water cluster can decompose on a picosecond timescale, resulting in a clearly different spectrum with the O-H stretch vibration becoming most important.
To summarize, the change of all these parameters affected the results relatively little. Sometimes a peak is missing. Obviously not all vibrational excitations can be observed for a single cluster in the gas phase on the picosecond timescale. Consequently, we put all simulations together in a single graph (Figure 4), omitting the outlier at 400 K. A similar procedure was applied before to compute the spectra of HCl and NH 3 [15]. We get a spectrum that compares very well to the experimental spectrum. All relevant features are observed. This procedure has some justification by the fact that we simulate isolated water clusters on the picosecond timescale. Without energy exchange with other clusters, even a way larger amount of CPU time would not lead to complete spectra. We can only do simulations with different initial conditions and simulation parameters.
To check if this procedure is reproducible, a third suite of simulations was run (Table 1). This time only the temperature was varied, all other parameters stayed the same. The result (Figure 5) is in beautiful agreement with experiment. In particular, the weak features around 1200 cm 1 and around 2000 cm 1 are now perfectly reproduced.

3. Discussion

AIMD simulations are able to describe IR spectra [19,20,21,22]. There is no need for a quantum mechanical approach. Of course, the electrons must be described using the Schrödinger equation, but there is no reason to do that also for the nuclei. In contrast, applying the Schrödinger equation to the nuclei leads to paradoxa. One of the most obvious examples that there is something wrong with the Copenhagen interpretation is actually the structure of ammonia. The ammonia inversion is described by a double-well potential. If the nuclei had a wavefunction in the sense of Schrödinger, both minima were equally occupied in the ground state, leading to a D 3 h structure, while the experiment clearly proves a C 3 v structure with a permanent dipole moment.
With the CPMD code we have a powerful tool in our hands, which, however, is not perfect. Deviations from the experiment are mainly attributed to shortcomings in the description of the electronic many-bodied problem, not in the classical description of nuclear motion. In the present example, the O-H stretch vibrations are shifted to lower wave numbers. This can be attributed to the use of the BLYP functional. BLYP, as a GGA (generalized gradient approximation) functional, has a very good cost-value ratio, but clearly does not have the accuracy of a hybrid functional. We are limited to the use of GGA functionals in our simulations because we use large plane-wave basis sets which makes the use of exact exchange problematic. Nevertheless, the BLYP error is preferable to the problems one encounters when moving atom-centered basis functions through space during dynamics.
In some graphs, not all features of the IR spectrum appear. We clearly have a problem with ergodicity stemming from limited CPU time. Nevertheless, from summing all simulations, it can be seen that the classical description of nuclear motion is in good agreement with the experimental spectrum. There is no peak that could not be explained with ab initio molecular dynamics. In view of practical calculations, however, it has some justification to use a quantum mechanical description of the vibrational wave function and to obtain a complete spectrum in a single calculation, instead of doing many similar CPMD simulations at non-zero temperatures. Additionally, of course, for small systems it is appealing not to be restricted to the accuracy of GGA functionals.
This observation does not alter the statement that nuclear motion is classical. At normal energies, i.e., as long as relativistic effects are not to be expected, Newton’s equations of motion are sufficient. An extension in the sense of special relativity is straight-forward, but there we get into a regime which is not covered by the daily experience of a chemist. Hence, we ignore relativity for the moment. This concerns of course the motion of the nuclei, not the electronic structure, where relativistic effects may be important.
What does a classical motion of the nuclei mean in practice? Our theory is local and deterministic. On this basis, some quantum mechanical paradoxa can be explained. Consider the thought experiment of separating a diatomic such as molecular hydrogen or hydrogen fluoride by irradiation with light. The products are separated from each other in closed boxes. In our local and deterministic theory, it is perfectly clear, at every point in time, which atom is on which side. This does not depend on the question, if the diatomic consists of like atoms such as in molecular hydrogen or of different atoms such as in hydrogen fluoride. The centers of these atoms move deterministicly, independent of a particular experiment. The electronic wavefunction surrounding these centers forms one big wavefunction with all the electrons of the boxes containing the atoms. It is not possible to define the box as something that would not let information through. It may be possible in a mathematical construct, to imagine that the inner of the box is not part of our universal wavefunction. In any physical experiment, the box and its contents are clearly part of our physical universe at every point in time. Hence, it is always clear, which atom is where.
Likewise, every atom of Schrödinger’s cat moves deterministicly at any point in time. From the heat flux through the box it is always clear, if it is dead or alive.
What about the double slit experiment? Feynman’s thought experiment [23] is particularly difficult to understand even if one is open to understanding the Copenhagen interpretation. Actually, it turns out that what Feynman describes is difficult to realize in an experiment, even if it is meanwhile possible to observe and manipulate single molecules. A pure double slit for atoms or molecules has never been constructed; the deBroglie wavelength is simply too small. Zeilinger’s experiments on fullerenes deal with a 100 nm grating, not with a double slit [24]. Fullerenes have a size in the order of 1 nm and this cannot be changed by any experiment, unless the molecules are destroyed. Fullerenes have certainly quantum properties, in the sense that one cannot cut them into halves without strongly changing their properties. In addition, they certainly have an electronic wavefunction surrounding them. However, if they are shot against a double slit, every atom will pass either the one slit or the other—unless the slits are smaller than 1 nm, then the fullerenes will hardly reach the other side. So, while Zeilinger’s experiment is interesting, it must be clearly stated, that the scenario is not (and is not intended to be) a confirmation of Feynman’s double-slit thought experiment.
What about the double-slit experiment for electrons? There is surprisingly little known in the literature, while lots of double-slit experiments have been performed for light. What would the thought experiment look like in a modern many-particle view? Electrons moving through a double slit become, at least for a short time, part of the electronic wave function of the double slit or rather of the material surrounding the double slit. It makes little sense to ask if they went through the one or the other slit. In a perfectly symmetric scenario, they may move through both slits with equal amounts at the same time. However, that is not the full story: depending on the electron energy, on the double-slit material, and on the double-slit geometry, a larger amount of electron density will pass through the surrounding material. Depending on the conditions, the electron may also be swallowed completely by the sea of electrons forming the double slit, which for simplicity is assumed to be grounded.
To conclude, the notion that nuclear motion is classical while the electronic cloud must be described with the time-independent Schrödinger equation, leads to a simple and compelling model of matter at normal energies. For other situations, the Maxwell equations may represent an important alternative. We still do not have the one single equation that describes everything. In any case, further calculations are necessary to show where the Schrödinger equation must be used for describing nuclear motion and where classical motion at finite temperatures is an alternative concept. Till now, we always found an alternative explanation, see, for example, [13].

4. Materials and Methods

Car–Parrinello molecular dynamics simulations using the CPMD code [1,25,26] have been performed in the NVE ensemble using the Becke–Lee–Yang–Parr (BLYP) functional [27,28]. In part of the simulations, the Grimme dispersion correction was applied [29]. The time step was chosen as 1 a.u. (0.0242 fs) and the fictitious electron mass as 100 a.u. These low values were used in order to guarantee an accurate motion of the light hydrogen atoms. Troullier–Martins pseudopotentials as optimized for the BLYP functional were employed for describing the core electrons [30,31]. The plane-wave cutoff which determines the size of the basis set, was set to 50.0 or 60.0 Rydberg. The simulation cell size was 24 × 24 × 24 a.u. 3 (12.70 × 12.70 × 12.70 Å 3 ) or 30 × 30 × 30 a.u. 3 (15.87 × 15.87 × 15.87 Å 3 ). Additionally, these values are relatively large for the small systems investigated. The aim is not to lose accuracy through limited basis sets and cell sizes. After geometry optimization, simulation runs were performed with initial kinetic energies of 300 K (simulations 1 to 9, see Table 1). Trajectories were determined for 1,000,000 steps per simulation run, that is, for 24 ps. In a second suite of simulations, the temperature was varied between 0 and 400 K, while in part of the calculations, the pseudopotential cutoff was changed and the Grimme correction was applied (simulations 10 to 18). The last suite of simulations (simulations 19 to 27) was performed as a check for temperatures varying between 0 and 400 K, while all other parameters remained unchanged. The TRAVIS program [22,32] was used to compute the spectra from the trajectories. The TRAVIS code computes the spectra from the autocorrelation of the dipole moments [21,33]. Both the contributions of nuclei and orbitals to the dipole moment are taken into account. The orbital contribution is included via localized Wannier orbitals, as obtained from AIMD simulations [34,35,36]. For the computation of the dipole moments, only the centers of these localized orbitals are needed.
Sample input files are provided for CPMD and for TRAVIS in the Supplementary Material. Sample job files for running CPMD in a parallel environment are also added.
Table 1 lists the parameters used in the single simulation runs.

5. Conclusions

If nuclear motion is classical, it must be possible to describe infrared spectra on the basis of ab initio molecular dynamics simulations. The present investigation of H 3 O + (H 2 O) 20 shows once more that ab initio molecular dynamics simulations present a reliable basis for the determination of infrared spectra. This confirms previous results [15]. It is important to state that a single simulation at the pico-second time scale is hardly sufficient to simulate a spectrum. Instead the results of several simulation runs with different starting conditions must be added together. This was found for gas-phase systems with relatively few degrees of freedom. For large systems in the condensed phase, a single simulation may be enough. However, it is also recommended in this case to test the effect of temperature. The accuracy of our calculations is limited by the use of density functionals in the GGA approximation. Nevertheless, for infrared spectra of large and complex systems, ab initio molecular dynamics simulations using GGA functionals can even be the best choice.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28186454/s1, sample CPMD version 4.3 input files in1 (equilibration) and in2 (production), sample job files job1 and job2, sample TRAVIS input file input.txt.

Funding

Part of the calculations were performed on the local cluster of the Leibniz University of Hannover at the LUIS, project number 424969120.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

More input files for CPMD version 4.3 are available on demand from the corresponding author.

Acknowledgments

The idea of the current project stems from a fruitful discussion with Kiyoshi Yagi which is gratefully acknowledged.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Sample Availability

Not applicable.

Abbreviations

The following abbreviations are used in this manuscript:
AIMDAb Initio Molecular Dynamics
BLYPBecke–Lee–Yang–Parr density functional
BOMDBorn–Oppenheimer Molecular Dynamics
CPMDCar–Parrinello Molecular Dynamics
DFTDensity Functional Theory
GGAGeneralized Gradient Approximation
KSKohn–Sham

References

  1. Car, R.; Parrinello, M. Unified Approach for Molecular Dynamics and Density-Functional Theory. Phys. Rev. Lett. 1985, 55, 2471–2474. [Google Scholar] [CrossRef]
  2. Marx, D.; Hutter, J. Ab Initio Molecular Dynamics: Theory and Implementation. In Modern Methods and Algorithms of Quantum Chemistry; Grotendorst, J., Ed.; NIC Series; Forschungszentrum Jülich: Jülich, Germany, 2000; Volume 1. [Google Scholar]
  3. Szabo, A.; Ostlund, N.S. Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory; Dover: Mineola, NY, USA, 1996. [Google Scholar]
  4. Marx, D.; Parrinello, M. Ab initio path-integral molecular dynamics. Z. Phys. B 1994, 95, 143. [Google Scholar] [CrossRef]
  5. Marx, D.; Parrinello, M. Structural quantum effects and three-centre two-electron bonding in CH 5 + . Nature 1995, 375, 216. [Google Scholar] [CrossRef]
  6. Marx, D.; Parrinello, M. The effect of quantum and thermal fluctuations on the structure of the floppy molecule C2 H 3 + . Science 1996, 271, 179. [Google Scholar] [CrossRef]
  7. Liu, Y.; Long, R.; Fang, W.H.; Prezhdo, O.V. Nuclear Quantum Effects Prolong Charge Carrier Lifetimes in Hybrid Organic-Inorganic Perovskites. J. Am. Chem. Soc. 2023, 145, 14112. [Google Scholar] [CrossRef]
  8. Tikhonov, D.S.; Vishnevskiy, J.V. Describing nuclear quantum effects in vibrational properties using molecular dynamics with Wigner sampling. Phys. Chem. Chem. Phys. 2023, 25, 18406. [Google Scholar] [CrossRef]
  9. Ple, T.; Mauger, N.; Adjoua, O.; Inizan, T.J.; Lagardere, L.; Huppert, S.; Piquemal, J.P. Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: From Force Fields to Machine Learning Potentials. J. Chem. Theory Comput. 2023, 19, 1432. [Google Scholar] [CrossRef]
  10. Frank, I.; Genuit, S.; Matz, F.; Oschinski, H. Ammonia, water, and hydrogen: Can nuclear motion be described classically? Int. J. Quantum Chem. 2020, 120, e26142. [Google Scholar] [CrossRef]
  11. Frank, I. Classical motion of the nuclei in a molecule: A concept without alternatives. Chem. Sel. 2020, 5, 1872. [Google Scholar] [CrossRef]
  12. Büchel, R.C.; Rudolph, D.A.; Frank, I. Deterministic quantum mechanics: The role of the Maxwell-Boltzmann distribution. Int. J. Quantum Chem. 2021, 121, e26555. [Google Scholar] [CrossRef]
  13. Rohloff, E.; Rudolph, D.A.; Strolka, O. Classical nuclear motion: Does it fail to explain reactions and spectra in certain cases? Int. J. Quantum Chem. 2022, 122, e26902. [Google Scholar] [CrossRef]
  14. Frank, I. Classical nuclear motion: Comparison to approaches with quantum mechanical nuclear motion. Hydrogen 2023, 4, 11–21. [Google Scholar] [CrossRef]
  15. Frank, I. Nuclear motion is classical: Spectra of hydrogen chloride and ammonia. Hydrogen 2023, 4, 287–294. [Google Scholar] [CrossRef]
  16. Yu, Q.; Pavosevic, F.; Hammes-Schiffer, S. Development of nuclear basis sets for multicomponent quantum chemistry methods. J. Chem. Phys. 2020, 152, 244123. [Google Scholar] [CrossRef]
  17. Liu, J.; Yang, J.; Zeng, X.C.; Xantheas, S.S.; Yagi, K.; He, X. Towards complete assignment of the infrared spectrum of the protonated water cluster H+(H2O)21. Nat. Commun. 2021, 12, 6141. [Google Scholar] [CrossRef]
  18. Yagi, K.; Thomson, B. Infrared spectra of protonated water clusters, H+(H2O)4, in Eigen and Zundel forms studied by vibrational quasi-degenerated perturbation theory. J. Phys. Chem. A 2017, 121, 2386. [Google Scholar] [CrossRef] [PubMed]
  19. Silvestrelli, P.L.; Bernasconi, M.; Parrinello, M. Ab initio infrared spectrum of liquid water. Chem. Phys. Lett. 1997, 277, 478–482. [Google Scholar] [CrossRef]
  20. Schmitz, M.; Tavan, P. Vibrational spectra from atomic fluctuations in dynamics simulations. I. Theory, limitations, and a sample application. J. Chem. Phys. 2004, 121, 12233. [Google Scholar] [CrossRef] [PubMed]
  21. Thomas, M.; Brehm, M.; Fligg, R.; Vöhringer, P.; Kirchner, B. Computing vibrational spectra from ab initio molecular dynamics. Phys. Chem. Chem. Phys. 2013, 15, 6608. [Google Scholar] [CrossRef] [PubMed]
  22. Brehm, M.; Thomas, M.; Gehrke, S.; Kirchner, B. TRAVIS—A free analyzer for trajectories from molecular simulation. J. Chem. Phys. 2020, 152, 164105. [Google Scholar] [CrossRef] [PubMed]
  23. Feynman, R.P.; Leighton, R.B.; Sands, M. The Feynman Lectures on Physics, 2nd ed.; Pearson Education: London, UK, 2006. [Google Scholar]
  24. Arndt, M.; Nairz, O.; Vos-Andreae, J.; Keller, C.; van der Zouw, G.; Zeilinger, A. Wave-particle duality of C60 molecules. Nature 1999, 401, 680. [Google Scholar] [CrossRef] [PubMed]
  25. Marx, D.; Hutter, J. Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  26. Hutter, J. CPMD Version 4.3, Copyright IBM Corp 1990–2015, Copyright MPI für Festkörperforschung Stuttgart 1997–2001. Available online: https://github.com/CPMD-code/CPMD/releases/tag/4.3 (accessed on 3 September 2023).
  27. Becke, A. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098. [Google Scholar] [CrossRef] [PubMed]
  28. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle–Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [PubMed]
  29. Grimme, S. Accurate Description of van der Waals Complexes by Density Functional Theory Including Empirical Corrections. J. Comput. Chem. 2004, 25, 1463. [Google Scholar] [CrossRef]
  30. Troullier, N.; Martins, J.L. Efficient Pseudopotentials for Plane-Wave Calculations. Phys. Rev. B 1991, 43, 1993. [Google Scholar] [CrossRef] [PubMed]
  31. Boero, M.; Parrinello, M.; Terakura, K.; Weiss, H. Car–Parrinello study of Ziegler-Natta heterogeneous catalysis: Stability and destabilization problems of the active site models. Mol. Phys. 2002, 100, 2935–2940. [Google Scholar] [CrossRef]
  32. Brehm, M.; Kirchner, B. TRAVIS—A free analyzer and visualizer for monte carlo and molecular dynamics trajectories. J. Chem. Inf. Model. 2011, 51, 2007. [Google Scholar] [CrossRef]
  33. Thomas, M.; Brehm, M.; Kirchner, B. Voronoi dipole moments for the simulation of bulk phase vibrational spectra. Phys. Chem. Chem. Phys. 2015, 17, 3207. [Google Scholar] [CrossRef]
  34. Marzari, N.; Vanderbilt, D. Maximally localized generalized Wannier functions for composite energy bands. Phys. Rev. B 1997, 56, 12847. [Google Scholar] [CrossRef]
  35. Ikeda, T.; Boero, M. Role of van der Waals corrections in first principles simulations of alkali metal ions in aqueous solutions. J. Chem. Phys. 2015, 143, 194510. [Google Scholar] [CrossRef]
  36. Silvestrelli, P.L.; Ambrosetti, A. Van der Waals interactions in DFT using Wannier functions without empirical parameters. J. Chem. Phys. 2019, 150, 164109. [Google Scholar] [CrossRef]
Figure 1. Structure of the water cluster H 3 O + (H 2 O) 20 . The H 3 O + is at the top of the cluster. Color code: red: oxygen, grey: hydrogen. The smaller one of the two unit cells is depicted.
Figure 1. Structure of the water cluster H 3 O + (H 2 O) 20 . The H 3 O + is at the top of the cluster. Color code: red: oxygen, grey: hydrogen. The smaller one of the two unit cells is depicted.
Molecules 28 06454 g001
Figure 2. First panel: Experimental spectrum. CPMD simulations: BLYP, temperature ions 300 K. Color code: black: plane-wave cutoff 50 Rydberg, red: plane-wave cutoff 60 Rydberg, green: plane-wave cutoff 60 Rydberg, dispersion correction added. Second panel: cell size 12.70 × 12.70 × 12.70 Å 3 ). Third panel: cell size 15.87 × 15.87 × 15.87 Å 3 . Forth panel: cell size 15.87 × 15.87 × 15.87 Å 3 , counterion added.
Figure 2. First panel: Experimental spectrum. CPMD simulations: BLYP, temperature ions 300 K. Color code: black: plane-wave cutoff 50 Rydberg, red: plane-wave cutoff 60 Rydberg, green: plane-wave cutoff 60 Rydberg, dispersion correction added. Second panel: cell size 12.70 × 12.70 × 12.70 Å 3 ). Third panel: cell size 15.87 × 15.87 × 15.87 Å 3 . Forth panel: cell size 15.87 × 15.87 × 15.87 Å 3 , counterion added.
Molecules 28 06454 g002
Figure 3. First panel: Experimental spectrum. CPMD simulations: BLYP, cell size 15.87 × 15.87 × 15.87 Å 3 . Color code: black: plane-wave cutoff 50 Rydberg, red: Plane-wave cutoff 60 Rydberg. Green: plane-wave cutoff 60 Rydberg, dispersion correction added. Second panel: start temperature of the ions is 0 K. Third panel: start temperature of the ions is 200 K. Forth panel: start temperature of the ions is 400 K. The outlier (red curve in the forth panel) is due to a decomposition of the cluster at high temperatures.
Figure 3. First panel: Experimental spectrum. CPMD simulations: BLYP, cell size 15.87 × 15.87 × 15.87 Å 3 . Color code: black: plane-wave cutoff 50 Rydberg, red: Plane-wave cutoff 60 Rydberg. Green: plane-wave cutoff 60 Rydberg, dispersion correction added. Second panel: start temperature of the ions is 0 K. Third panel: start temperature of the ions is 200 K. Forth panel: start temperature of the ions is 400 K. The outlier (red curve in the forth panel) is due to a decomposition of the cluster at high temperatures.
Molecules 28 06454 g003
Figure 4. Colored lines: IR spectra of simulation runs 1 to 18 The outlier 17, which was observed at high temperatures, was omitted. Black line: sum over all these spectra.
Figure 4. Colored lines: IR spectra of simulation runs 1 to 18 The outlier 17, which was observed at high temperatures, was omitted. Black line: sum over all these spectra.
Molecules 28 06454 g004
Figure 5. Colored lines: IR spectra of all simulation 19 to 27. Black line: sum over all these spectra.
Figure 5. Colored lines: IR spectra of all simulation 19 to 27. Black line: sum over all these spectra.
Molecules 28 06454 g005
Table 1. Parameters in the single simulation runs.
Table 1. Parameters in the single simulation runs.
Number of Simulation RunCell SizeTemperaturePlane Wave CutoffDispersion CorrectionCounter Ion
[Ångstrom][Kelvin][Rydberg]
112.7030050nono
212.7030060nono
312.7030060yesno
415.8730050nono
515.8730060nono
615.8730060yesno
715.8730050noyes
815.8730060noyes
915.8730060yesyes
1015.87050nono
1115.87060nono
1215.87060yesno
1315.8720050nono
1415.8720060nono
1515.8720060yesno
1615.8740050nono
1715.8740060nono
1815.8740060yesno
1912.70050nono
2012.705050nono
2112.7010050nono
2212.7015050nono
2312.7020050nono
2412.7025050nono
2512.7030050nono
2612.7035050nono
2712.7040050nono
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Frank, I. Nuclear Motion Is Classical: Spectrum of a Magic Protonated Water Cluster. Molecules 2023, 28, 6454. https://doi.org/10.3390/molecules28186454

AMA Style

Frank I. Nuclear Motion Is Classical: Spectrum of a Magic Protonated Water Cluster. Molecules. 2023; 28(18):6454. https://doi.org/10.3390/molecules28186454

Chicago/Turabian Style

Frank, Irmgard. 2023. "Nuclear Motion Is Classical: Spectrum of a Magic Protonated Water Cluster" Molecules 28, no. 18: 6454. https://doi.org/10.3390/molecules28186454

APA Style

Frank, I. (2023). Nuclear Motion Is Classical: Spectrum of a Magic Protonated Water Cluster. Molecules, 28(18), 6454. https://doi.org/10.3390/molecules28186454

Article Metrics

Back to TopTop