Four-body Terms Research Articles

Microsolvation of a Proton by Ar Atoms: Structures and Energetics of ArnH+ Clusters.

We present a computational investigation on the structural arrangements and energetic stabilities of small-size protonated argon clusters, Ar nH +. Using high-level ab initio electronic structure computations, we determined that the linear symmetric triatomic ArH +Ar ion serves as the molecular core for all larger clusters studied. Through harmonic normal-mode analysis for clusters containing up to seven argon atoms, we observed that the proton-shared vibration shifts to lower frequencies, consistent with measurements in gas-phase IRPD and solid Ar-matrix isolation experiments. We explored the sum-of-potentials approach by employing kernel-based machine-learning potential models trained on CCSD(T)-F12 data. These models included expansions of up to two-body, three-body, and four-body terms to represent the underlying interactions as the number of Ar atoms increases. Our results indicate that the four-body contributions are crucial for accurately describing the potential surfaces in clusters with n> 3. Using these potential models and an evolutionary programming method, we analyzed the structural stability of clusters with up to 24 Ar atoms. The most energetically favored Ar nH + structures were identified for magic size clusters at n = 7, 13, and 19, corresponding to the formation of Ar-pentagon rings perpendicular to the ArH +Ar core ion axis. The sequential formation of such regular shell structures is compared to ion yield data from high-resolution mass spectrometry measurements. Our results demonstrate the effectiveness of the developed sum-of-potentials model in describing trends in the nature of bonding during the single proton microsolvation by Ar atoms, encouraging further quantum nuclear studies.

An Atomistic-Based Nonlinear Plate Theory for Hexagonal Boron Nitride.

Through the continuity of the DREIDING force field, we propose, for the first time, the finite-deformation plate theory for the single-layer hexagonal boron nitride (h-BN) to clarify the atomic source of the structure against deformations. Divergent from the classical Föppl-von Karman plate theory, our new theory shows that h-BN’s two in-plane mechanical parameters are independent of two out-of-plane mechanical parameters. The new theory reveals the relationships between the h-BN’s elastic rigidities and the atomic force field: (1) two in-plane elastic rigidities come from the bond stretching and the bond angle bending; (2) the bending rigidity comes from the inversion angle and the dihedral angle torsion; (3) the Gaussian rigidity only comes from the dihedral angle torsion. Mechanical parameters obtained by our theory align with atomic calculations. The new theory proves that two four-body terms in the DREIDING force field are necessary to model the h-BN’s mechanical properties. Overall, our theory establishes a foundation to apply the classical plate theory on the h-BN, and the approach in this paper is heuristic in modelling the mechanical properties of the other two-dimensional nanostructures.

Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD).

We present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-molecular dynamics (MD) classical and nuclear quantum dynamical simulations, in which the energies and forces of the full system are approximated by the two-, three-, and four-body terms of the MBE, and compare the average potential and the vibrational density of states with the full simulation, i.e., the one for which no MBE is used. Our results indicate that the thermally averaged potential energy from the MBE up to the four-body term converges with near-identical behavior to the one from the full simulation. The three-body makes a substantial contribution (∼20%) to the energy, whereas the four-body is necessary for obtaining quantitatively accurate energetics and forces, albeit making a small contribution to each (∼2%). We further show that the harmonic frequencies are reproduced to within a few wavenumbers (cm-1) at the four-body level and that the slowest modes to converge with the MBE rank are those involving the strongest hydrogen bonds. Anharmonicity exacerbates this effect, so that a four-body description of the energies and forces is needed to achieve accurate anharmonic vibrational frequencies in the hydrogen-bonded OH-stretching region. We also discuss the asymptotic scaling of the MBE-MD protocol with respect to the cost of the underlying potential energy evaluation, suggesting that electronic structure methods that scale at least as N4, N being the size of the system, are needed to result in savings over the traditional full MD simulation. We anticipate that the MBE-MD protocol can evolve into a powerful and practical method, which will allow for highly accurate ab initio MD simulations on a much broader range of molecular systems than can be currently handled.

Accessing the applicability of the MBE approach for constructing potential energy surfaces of nitrogen clusters

The many-body expansion approach (MBE) is employed to build potential energy surfaces (PESs) for nitrogen clusters, which serve as a possible candidate of high energy density materials. The ground-state global PES of N 4 developed by Truhlar’s group is utilized to construct the two-body, three-body, and four-body interaction potentials. The energy of nitrogen clusters is thus obtained by summing the many-body interaction energies up to four-body term. This work found that such a many-body expansion could give reasonable energy predictions for loose chain nitrogen molecules while produce bad predictions for compact structures. The failure is caused by the neglect of high-order ( n ≥ 5) n -body interactions that are significant in the bonding area. The difficulties in applying the MBE method to construct the potential energy surface of nitrogen clusters are analyzed.

Cross-Investigation on Copper Nitroprusside: Combining XRD and XAS for In-Depth Structural Insights

The emerging energy demand and need to develop sustainable energy storage systems have drawn extensive attention to fundamental and applied research. Anion redox processes were proposed in cathodic materials in addition to traditional transition metal redox to boost the specific capacity and the electrochemical performance. Alternatively, copper nitroprusside (CuNP) features an electroactive nitrosyl ligand alongside the two structural metals (Fe, Cu), representing an alternative to anion redox in layered oxides. Here, a deep structural investigation is carried out on CuNP by complementing the long-range order sensitivity of X-ray diffraction (XRD) and the local atomic probe of X-ray absorption (XAS). Two different CuNP materials are studied, the hydrated and dehydrated forms. A new phase for hydrated CuNP not reported in the literature is solved, and Rietveld refined. The XAS spectra of the two materials at the Cu and Fe K-edges show a similar yet different atomic environment. The extended XAS spectra (EXAFS) analysis is accomplished by considering three- and four-body terms due to the high collinearity of the atomic chains and gives accurate insight into the first-, second-, and third-shell interatomic distances. Both materials are mounted in Li-ion and Na-ion cells to explore the link between structure and electrochemical performance. As revealed by the charge/discharge cycles, the cyclability in Na-ion cells is negatively affected by interstitial water. The similarity in the local environment and the electrochemical differences suggest a long-range structural dependence on the electrochemical performance.

Path-integral calculation of the fourth virial coefficient of helium isotopes.

We use the path-integral Monte Carlo (PIMC) method and state-of-the-art two-body and three-body potentials to calculate the fourth virial coefficients D(T) of 4He and 3He as functions of temperature from 2.6 K to 2000 K. We derive expressions for the contributions of exchange effects due to the bosonic or fermionic nature of the helium isotope; these effects have been omitted from previous calculations. The exchange effects are relatively insignificant for 4He at the temperatures considered, but for 3He, they are necessary for quantitative accuracy below about 4 K. Our results are consistent with previous theoretical work (also with some of the limited and scattered experimental data) for 4He; for 3He, there are no experimental values, and this work provides the first values of D(T) calculated at this level. The uncertainty of the results depends on the statistical uncertainty of the PIMC calculation, the estimated effect of omitting four-body terms in the potential energy, and the uncertainty contribution propagated from the uncertainty of the potentials. At low temperatures, the uncertainty is dominated by the statistical uncertainty of the PIMC calculations, while at high temperatures, the uncertainties related to the three-body potential and omitted higher-order contributions become dominant.

Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage.

The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH4···(H2O)20, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.

A new many-body expansion scheme for atomic clusters: Application to nitrogen clusters

Although the many-body expansion (MBE) approach is widely applied to estimate the energy of large systems containing weak interactions, it is inapplicable to calculating the energies of covalent or metal clusters. In this work, we propose an interaction many-body expansion (IMBE) to calculate the energy of atomic clusters containing covalent bonds. In this approach, the energy of a system is expressed as the sum of the energy of atoms and the interaction energy between the atom and its surrounding atoms. The IMBE method is first applied to calculate the energies of nitrogen clusters, in which the interatomic interactions are truncated to four-body terms. The results show that the IMBE approach could significantly reduce the energy error for nitrogen clusters compared with the traditional MBE method. The weak size and structure dependence of the IMBE error with respect to DFT calculations indicates the IMBE method has good potential application in estimating energy of large covalent systems.

Calculation of the detonation state of HN3 with quantum accuracy.

HN3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN3 with the accuracy of Kohn-Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN3. We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn-Sham density-functional theory.

Two-Dimensional Quantum-Link Lattice Quantum Electrodynamics at Finite Density

We present an unconstrained tree tensor network approach to the study of lattice gauge theories in two spatial dimensions showing how to perform numerical simulations of theories in presence of fermionic matter and four-body magnetic terms, at zero and finite density, with periodic and open boundary conditions. We exploit the quantum link representation of the gauge fields and demonstrate that a fermionic rishon representation of the quantum links allows us to efficiently handle the fermionic matter while finite densities are naturally enclosed in the tensor network description. We explicit perform calculations for quantum electrodynamics in the spin-one quantum link representation on lattice sizes of up to 16x16 sites, detecting and characterizing different quantum regimes. In particular, at finite density, we detect signatures of a phase separation as a function of the bare mass values at different filling densities. The presented approach can be extended straightforwardly to three spatial dimensions.

Long-Range Multibody Interactions and Three-Body Antiblockade in a Trapped Rydberg Ion Chain.

Trapped Rydberg ions represent a flexible platform for quantum simulation and information processing that combines a high degree of control over electronic and vibrational degrees of freedom. The possibility to individually excite ions to high-lying Rydberg levels provides a system where strong interactions between pairs of excited ions can be engineered and tuned via external laser fields. We show that the coupling between Rydberg pair interactions and collective motional modes gives rise to effective long-range and multibody interactions consisting of two, three, and four-body terms. Their shape, strength, and range can be controlled via the ion trap parameters and strongly depends on both the equilibrium configuration and vibrational modes of the ion crystal. By focusing on an experimentally feasible quasi one-dimensional setup of ^{88}Sr^{+} Rydberg ions, we demonstrate that multibody interactions are enhanced by the emergence of soft modes associated with, e.g., a structural phase transition. This has a striking impact on many-body electronic states and results-for example-in a three-body antiblockade effect that can be employed as a sensitive probe to detect structural phase transitions in Rydberg ion chains. Our study unveils the possibilities offered by trapped Rydberg ions for studying exotic phases of matter and quantum dynamics driven by enhanced multibody interactions.

Development and initial applications of an e-ReaxFF description of Ag nanoclusters.

Metal nanocrystals are of considerable scientific interest because of their uses in electronics, catalysis, and spectroscopy, but the mechanisms by which nanocrystals nucleate and grow to achieve selective shapes are poorly understood. Ab initio calculations and experiments have consistently shown that the lowest energy isomers for small silver nanoparticles exhibit two-dimensional (2D) configurations and that a transition into three-dimensional (3D) configurations occurs with the addition of only a few atoms. We parameterized an e-ReaxFF potential for Ag nanoclusters (N ≤ 20 atoms) that accurately reproduces the 2D-3D transition observed between the Ag5 and Ag7 clusters. This potential includes a four-body dihedral term that imposes an energetic penalty to 3D structures that is significant for small clusters but is overpowered by the bond energy from out-of-plane Ag-Ag bonds in larger 3D clusters. The potential was fit to data taken from density-functional theory and coupled-cluster calculations and compared to an embedded atom method potential to gauge its quality. We also demonstrate the potential of e-ReaxFF to model redox reactions in silver halides and plasmon motion using molecular dynamics simulations. This is the first case in which e-ReaxFF is used to describe metals. Furthermore, the inclusion of a bond-order dependent dihedral angle in this force field is a unique solution to modeling the 2D-3D transition seen in small metal nanoclusters.

Mediation of resonance energy transfer by two polarisable particles.

The molecular quantum electrodynamics theory is employed to calculate the matrix element and Fermi golden rule rate for resonant transfer of electronic excitation energy between a donor and an acceptor in the vicinity of two neutral electric dipole polarizable particles, which play the role of bridging species. The emitter and absorber couple linearly to the electric displacement field via their electric dipole moments, while each mediator interacts quadratically with this field through its dynamic polarizability. This form of interaction Hamiltonian enables fourth-order perturbation theory to be used to compute the probability amplitude together with summation over 24 time-ordered diagrams representing a single virtual photon exchange between each pair of coupled particles. Expressions for the migration rate mediated by two inert molecules are obtained for an arbitrary arrangement of the four species that are in fixed mutual orientation or are freely tumbling. These formulae are valid for all interparticle separation distances outside the orbital overlap region. From the general result, rate equations applicable to an equidistant collinear configuration of the four bodies are evaluated. Near- and far-zone limiting forms of the transfer rate for the relay pathway are also calculated and exhibit inverse sixth and inverse square dependences on relative separation distances between pairs of particles, confirming the short-range (radiationless) and long-range (radiative) energy transfer mechanisms associated with two-body theory. The distance behavior of interference terms between two-, three-, and four-body terms is also examined, and the relative importance of each contribution to the total transfer rate is discussed.

First-Principles Study on Crystal Configuration and Many-Body Cohesive Energy of Solid Argon

Based on the first-principles, using CCSD(T) ab initio calculation method, many-body potential energy of solid argon are accurately calculated with the atomic distance R from 2.0Å to 3.6Å at T=300K, and firstly establish and discuss the face-centered cubic (fcc) atomic crystal configurations of two-, three-, and four-body terms by geometry optimization. The results shows that the total number of (Ar)2 clusters is 903, which belongs to 12 different geometric configurations, the total number of (Ar)3 clusters is 861, which belongs to 25 different geometric configurations, and the total number of (Ar)4 clusters of is 816 which belongs to 27 different geometric configurations. We find that the CCSD(T) with the aug-cc-pVQZ basis set is most accurate and practical by comprehensive consideration. The total potential energy Un reachs saturation at R>2.0Å when the only two-and three-body interaction energy are considered. When R≤2.0Å, the total potential energy Un must consider four-and higher-body interaction energy to achieve saturation. Many-body expansion potential of fcc solid argon is an exchange convergent series.

Difficulties and Virtues in Assessing the Potential Energy Surfaces of Carbon Clusters via DMBE Theory: Stationary Points of Cκ (κ = 2-10) at the Focal Point.

The previously reported potential energy surfaces (PESs) of ground-state singlet C3 and triplet C4 are here utilized as input for the construction of approximate cluster expansions for larger Cκ (κ = 5-10) species. Relying primarily on the double many-body expansion (DMBE) approach, global potentials are obtained by summing the total interaction energies of all atomic subclusters up to four-body terms. The notable capability of the final forms in predicting good estimates of the linear global minima and their thermochemical/structural properties may provide important insights into the structure-determining nature of the (2 + 3 + 4) terms. The main difficulties and virtues in assessing Cκ's via DMBE theory are analyzed and new prospects given for the construction of global reliable PESs for the target species.

Four-body Interaction and Equation of State for Solid Neon from Ab Initio Calculation

Using ab initio Hartree-Fock self-consistent field method combined with many-body expansion method, the investigation is based on the first-principles. We have considered two-, three- and four-body potential energies of face-centered cubic (fcc) solid neon of which the atomic distance ranges from 1.6 to 3.0 . By discussing the truncation and convergence of many-body potential of solid neon, we obtain the cohesive energy, the zero-point vibration energy and equation of state (EOS). The results show that, when the number of neighboring atoms increases, two-body, three-body, and four-body potential energy tend to a saturation value for a certain atomic distance ( ). The even many-body contributions to the cohesive energy, such as two-, four-body terms and so on, are positive, whereas the odd many-body contributions to the cohesive energy, such as three-, five-body terms and so on, are negative. The zero-point vibration energy of solid neon is only 6% of the total atomic interaction energy, but should not be neglected. Compared with the experimental data, the importance of the four-body interactions in compressed solid neon is emphasized. Only taking into account the two-body term, the pressure is overestimated, and our calculated results are in good agreement with the experimental values at the low-pressure regions (<15GPa). Adding three-body term up to 55GPa, considering the four-body term, it has a good consistency at the experimentally studied pressure range of 0~237GPa, and maybe helpful to accurately explain the phenomenon of the experiment above 237GPa when the higher many-body effects are considered.

A theoretical study of energy transfer in Ar(1S) + SO2( ') collisions: Cross sections and rate coefficients for vibrational transitions.

Vibrational transitions, induced by collisions between rare-gas atoms and molecules, play a key role in many problems of interest in physics and chemistry. A theoretical investigation of the translation-to-vibration (T-V) energy transfer process in argon atom and sulfur dioxide molecule collisions is presented here. For such a purpose, the framework of the quasi-classical trajectory (QCT) methodology was followed over the range of translational energies 2 ≤ Etr/kcal mol-1 ≤ 100. A new realistic potential energy surface (PES) for the ArSO2 system was developed using pairwise addition for the four-body energy term within the double many-body expansion. The topological features of the obtained function are compared with a previous one reported by Hippler et al. [J. Phys. Chem. 90, 6158 (1986)]. To test the accuracy of the PES, additional coupled cluster singles and doubles method with a perturbative contribution of connected triples calculations were carried out for the global minimum configuration. From dynamical calculations, the cross sections for the T-V excitation process indicate a barrier-type mechanism due to strong repulsive interactions between SO2 molecules and the Ar atom. Corrections to zero-point energy leakage in QCT were carried out using vibrational energy quantum mechanical threshold of the complex and variations. Rate coefficients and cross sections are calculated for some vibrational transitions using pseudo-quantization approaches of the vibrational energy of products. Main attributes of the title molecular collision are discussed and compared with available information in the literature.

Cn (n=2-4): current status.

The major aspects of the C2, C3 and C4 elemental carbon clusters are surveyed. For C2, a brief analysis of its current status is presented. Regarding C3, the most recent results obtained in our group are reviewed with emphasis on modelling its potential energy surface which is particularly complicated due to the presence of multiple conical intersections. As for C4, the most stable isomeric forms of both triplet and singlet spin states and their possible interconversion pathways are examined afresh by means of accurate ab initio calculations. The main strategies for modelling the ground triplet C4 potential are also discussed. Starting from a truncated cluster expansion and a previously reported DMBE form for C3, an approximate four-body term is calibrated from the ab initio energies. The final six-dimensional global DMBE form so obtained reproduces all known topographical aspects while providing an accurate description of the C4 linear-rhombic isomerization pathway. It is therefore commended for both spectroscopic and reaction dynamics studies.This article is part of the theme issue 'Modern theoretical chemistry'.

Understanding the many-body expansion for large systems. III. Critical role of four-body terms, counterpoise corrections, and cutoffs.

Papers I and II in this series [R. M. Richard et al., J. Chem. Phys. 141, 014108 (2014); K. U. Lao et al., ibid. 144, 164105 (2016)] have attempted to shed light on precision and accuracy issues affecting the many-body expansion (MBE), which only manifest in larger systems and thus have received scant attention in the literature. Many-body counterpoise (CP) corrections are shown to accelerate convergence of the MBE, which otherwise suffers from a mismatch between how basis-set superposition error affects subsystem versus supersystem calculations. In water clusters ranging in size up to (H2O)37, four-body terms prove necessary to achieve accurate results for both total interaction energies and relative isomer energies, but the sheer number of tetramers makes the use of cutoff schemes essential. To predict relative energies of (H2O)20 isomers, two approximations based on a lower level of theory are introduced and an ONIOM-type procedure is found to be very well converged with respect to the appropriate MBE benchmark, namely, a CP-corrected supersystem calculation at the same level of theory. Results using an energy-based cutoff scheme suggest that if reasonable approximations to the subsystem energies are available (based on classical multipoles, say), then the number of requisite subsystem calculations can be reduced even more dramatically than when distance-based thresholds are employed. The end result is several accurate four-body methods that do not require charge embedding, and which are stable in large basis sets such as aug-cc-pVTZ that have sometimes proven problematic for fragment-based quantum chemistry methods. Even with aggressive thresholding, however, the four-body approach at the self-consistent field level still requires roughly ten times more processors to outmatch the performance of the corresponding supersystem calculation, in test cases involving 1500-1800 basis functions.

A transmon quantum annealer: decomposing many-body Ising constraints into pair interactions

Adiabatic quantum computing is an analogue quantum computing scheme with various applications in solving optimisation problems. In the parity picture of quantum optimization, the problem is encoded in local fields that act on qubits that are connected via local four-body terms We present an implementation of a parity annealer with Transmon qubits with a specifically tailored Ising interaction from Josephson ring modulators.