Anharmonic Systems Research Articles

Efficient and Adaptive Methods for Computing Accurate Potential Surfaces for Quantum Nuclear Effects: Applications to Hydrogen-Transfer Reactions.

We present two sampling measures to gauge critical regions of potential energy surfaces. These sampling measures employ (a) the instantaneous quantum wavepacket density, an approximation to the (b) potential surface, its (c) gradients, and (d) a Shannon information theory based expression that estimates the local entropy associated with the quantum wavepacket. These four criteria together enable a directed sampling of potential surfaces that appears to correctly describe the local oscillation frequencies, or the local Nyquist frequency, of a potential surface. The sampling functions are then utilized to derive a tessellation scheme that discretizes the multidimensional space to enable efficient sampling of potential surfaces. The sampled potential surface is then combined with four different interpolation procedures, namely, (a) local Hermite curve interpolation, (b) low-pass filtered Lagrange interpolation, (c) the monomial symmetrization approximation (MSA) developed by Bowman and co-workers, and (d) a modified Shepard algorithm. The sampling procedure and the fitting schemes are used to compute (a) potential surfaces in highly anharmonic hydrogen-bonded systems and (b) study hydrogen-transfer reactions in biogenic volatile organic compounds (isoprene) where the transferring hydrogen atom is found to demonstrate critical quantum nuclear effects. In the case of isoprene, the algorithm discussed here is used to derive multidimensional potential surfaces along a hydrogen-transfer reaction path to gauge the effect of quantum-nuclear degrees of freedom on the hydrogen-transfer process. Based on the decreased computational effort, facilitated by the optimal sampling of the potential surfaces through the use of sampling functions discussed here, and the accuracy of the associated potential surfaces, we believe the method will find great utility in the study of quantum nuclear dynamics problems, of which application to hydrogen-transfer reactions and hydrogen-bonded systems is demonstrated here.

Anomalous interfacial temperature profile induced by phonon localization.

Through the integration of the power spectral density, we obtain temperature profiles of both multisegment harmonic and anharmonic systems, showing the presence of an anomalous negative temperature gradient inside the interfacial segment. Via investigating patterns of the power spectral density, we found that the counterintuitive phenomenon comes from the presence of interfacial localized phonon modes. Two out-band localized modes of the harmonic model, which make no contributions to local temperature due to the absence of phonon interactions, result in the concave temperature profile and overcooling effect. For the anharmonic model, thanks to the phonon-phonon interactions, the localized modes are excited and make considerable contributions to interfacial temperature, which is clearly shown by examining the temperature accumulation function. When anharmonicity is considerably large, the negative temperature gradient is absent since the localized phonon modes are fully mixed. The presence of localized modes are evidently demonstrated by the inverse participation ratio and normal mode analysis for the isolated harmonic model. The localized modes make contribution to interfacial temperature profiles of the harmonic system when they are excited in initial conditions of simulations.

Harmonic-phase path-integral approximation of thermal quantum correlation functions.

We present an approximation to the thermal symmetric form of the quantum time-correlation function in the standard position path-integral representation. By transforming to a sum-and-difference position representation and then Taylor-expanding the potential energy surface of the system to second order, the resulting expression provides a harmonic weighting function that approximately recovers the contribution of the phase to the time-correlation function. This method is readily implemented in a Monte Carlo sampling scheme and provides exact results for harmonic potentials (for both linear and non-linear operators) and near-quantitative results for anharmonic systems for low temperatures and times that are likely to be relevant to condensed phase experiments. This article focuses on one-dimensional examples to provide insights into convergence and sampling properties, and we also discuss how this approximation method may be extended to many-dimensional systems.

Anharmonic quantum mechanical systems do not feature phase space trajectories

Phase space dynamics in classical mechanics is described by transport along trajectories. Anharmonic quantum mechanical systems do not allow for a trajectory-based description of their phase space dynamics. This invalidates some approaches to quantum phase space studies. We first demonstrate the absence of trajectories in general terms. We then give an explicit proof for all quantum phase space distributions with negative values: we show that the generation of coherences in anharmonic quantum mechanical systems is responsible for the occurrence of singularities in their phase space velocity fields, and vice versa. This explains numerical problems repeatedly reported in the literature, and provides deeper insight into the nature of quantum phase space dynamics.

Quantum work fluctuations in connection with the Jarzynski equality

A result of great theoretical and experimental interest, the Jarzynski equality predicts a free energy change ΔF of a system at inverse temperature β from an ensemble average of nonequilibrium exponential work, i.e., 〈e^{-βW}〉=e^{-βΔF}. The number of experimental work values needed to reach a given accuracy of ΔF is determined by the variance of e^{-βW}, denoted var(e^{-βW}). We discover in this work that var(e^{-βW}) in both harmonic and anharmonic Hamiltonian systems can systematically diverge in nonadiabatic work protocols, even when the adiabatic protocols do not suffer from such divergence. This divergence may be regarded as a type of dynamically induced phase transition in work fluctuations. For a quantum harmonic oscillator with time-dependent trapping frequency as a working example, any nonadiabatic work protocol is found to yield a diverging var(e^{-βW}) at sufficiently low temperatures, markedly different from the classical behavior. The divergence of var(e^{-βW}) indicates the too-far-from-equilibrium nature of a nonadiabatic work protocol and makes it compulsory to apply designed control fields to suppress the quantum work fluctuations in order to test the Jarzynski equality.

Non-classicality from the phase-space flow analysis of the Weyl-Wigner quantum mechanics

A fluid analog of the information flux in the phase-space associated to purity and von Neumann entropy is identified in the Weyl-Wigner formalism of quantum mechanics. Once constrained by symmetry and positiveness, the encountered continuity equations provide novel quantifiers for non-classicality (non-Liouvillian fluidity) given in terms of quantum decoherence, purity and von Neumann entropy fluxes. Through definitions in the Weyl-Wigner formalism, one can identify the quantum fluctuations that distort the classical-quantum coincidence regime, and the corresponding quantum information profile, whenever some bounded x − p volume of the phase-space is specified. The dynamics of anharmonic systems is investigated in order to illustrate such a novel paradigm for describing quantumness and classicality through the flux of quantum information in the phase-space.

Thermal transport through a spin-phonon interacting junction: A nonequilibrium Green's function method study

Using the nonequilibrium Green's function method, we consider heat transport in an insulating ferromagnetic spin chain model with spin-phonon interaction under an external magnetic field. Employing the Holstein-Primakoff transformation to the spin system, we treat the resulted magnon-phonon interaction within the self-consistent Born approximation. We find the magnon-phonon coupling can change qualitatively the magnon thermal conductance in the high-temperature regime. At a spectral mismatched ferromagnetic-normal insulator interface, we also find thermal rectification and negative differential thermal conductance due to the magnon-phonon interaction. We show that these effects can be effectively tuned by the external applied magnetic field, a convenient advantage absent in anharmonic phonon and electron-phonon systems studied before.

Wigner\u2019s quantum phase-space current in weakly-anharmonic weakly-excited two-state systems

.There are no phase-space trajectories for anharmonic quantum systems, but Wigner’s phase-space representation of quantum mechanics features Wigner current J . This current reveals fine details of quantum dynamics —finer than is ordinarily thought accessible according to quantum folklore invoking Heisenberg’s uncertainty principle. Here, we focus on the simplest, most intuitive, and analytically accessible aspects of J. We investigate features of J for bound states of time-reversible, weakly-anharmonic one-dimensional quantum-mechanical systems which are weakly-excited. We establish that weakly-anharmonic potentials can be grouped into three distinct classes: hard, soft, and odd potentials. We stress connections between each other and the harmonic case. We show that their Wigner current fieldline patterns can be characterised by J’s discrete stagnation points, how these arise and how a quantum system’s dynamics is constrained by the stagnation points’ topological charge conservation. We additionally show that quantum dynamics in phase space, in the case of vanishing Planck constant ℏ or vanishing anharmonicity, does not pointwise converge to classical dynamics.

Transition state theory for activated systems with driven anharmonic barriers.

Classical transition state theory has been extended to address chemical reactions across barriers that are driven and anharmonic. This resolves a challenge to the naive theory that necessarily leads to recrossings and approximate rates because it relies on a fixed dividing surface. We develop both perturbative and numerical methods for the computation of a time-dependent recrossing-free dividing surface for a model anharmonic system in a solvated environment that interacts strongly with an oscillatory external field. We extend our previous work, which relied either on a harmonic approximation or on periodic force driving. We demonstrate that the reaction rate, expressed as the long-time flux of reactive trajectories, can be extracted directly from the stability exponents, namely, Lyapunov exponents, of the moving dividing surface. Comparison to numerical results demonstrates the accuracy and robustness of this approach for the computation of optimal (recrossing-free) dividing surfaces and reaction rates in systems with Markovian solvation forces. The resulting reaction rates are in strong agreement with those determined from the long-time flux of reactive trajectories.

Fine tuning classical and quantum molecular dynamics using a generalized Langevin equation.

Generalized Langevin Equation (GLE) thermostats have been used very effectively as a tool to manipulate and optimize the sampling of thermodynamic ensembles and the associated static properties. Here we show that a similar, exquisite level of control can be achieved for the dynamical properties computed from thermostatted trajectories. We develop quantitative measures of the disturbance induced by the GLE to the Hamiltonian dynamics of a harmonic oscillator, and show that these analytical results accurately predict the behavior of strongly anharmonic systems. We also show that it is possible to correct, to a significant extent, the effects of the GLE term onto the corresponding microcanonical dynamics, which puts on more solid grounds the use of non-equilibrium Langevin dynamics to approximate quantum nuclear effects and could help improve the prediction of dynamical quantities from techniques that use a Langevin term to stabilize dynamics. Finally we address the use of thermostats in the context of approximate path-integral-based models of quantum nuclear dynamics. We demonstrate that a custom-tailored GLE can alleviate some of the artifacts associated with these techniques, improving the quality of results for the modeling of vibrational dynamics of molecules, liquids, and solids.

Wavepacket revivals via complex trajectory propagation

Complex-valued semiclassical methods hold out the promise of treating classically allowed and classically forbidden processes on the same footing. In addition, they provide a natural way to describe optical excitation with complex fields within the trajectory framework. Despite their promise, these methods have until now been limited to short time propagation, due to the numerical difficulties introduced by the complexification. Using a new Final Value Representation of the Coherent State Propagator (FINCO), combined with an analysis of the complex classical phase space, we achieve accurate wavepacket propagation all the way to the revival time of a strongly anharmonic system.

Manipulating Fock states of a harmonic oscillator while preserving its linearity

We present a new scheme for controlling the quantum state of a harmonic oscillator by coupling it to an anharmonic multilevel system (MLS) with first to second excited state transition frequency on-resonance with the oscillator. In this scheme that we call "ef-resonant", the spurious oscillator Kerr non-linearity inherited from the MLS is very small, while its Fock states can still be selectively addressed via an MLS transition at a frequency that depends on the number of photons. We implement this concept in a circuit-QED setup with a microwave 3D cavity (the oscillator, with frequency 6.4 GHz and quality factor QO=2E-6) embedding a frequency tunable transmon qubit (the MLS). We characterize the system spectroscopically and demonstrate selective addressing of Fock states and a Kerr non-linearity below 350 Hz. At times much longer than the transmon coherence times, a non-linear cavity response with driving power is also observed and explained.

ChemInform Abstract: Lead‐Free Tin Chalcogenide Thermoelectric Materials

AbstractReview: 150 refs.

Approaching the strongly anharmonic limit with ab initio calculations of materials’ vibrational properties – a colloquium*

Despite ab initio computational techniques have opened new possibilities to interpret experimental results and predict the properties of new materials, their applications are limited by the adopted approximative schemes. Consequently, the first-principles calculation of many physical properties and phenomena is hindered and ab initio methods need to be further developed to overcome such limits. For example, the standard harmonic approximation used to assess the vibrational properties of materials often completely breaks down, so that the vibrational properties need to be calculated including strong anharmonic effects. The harmonic approximation has also intrinsic failures as it cannot estimate the lattice thermal conductivity of materials nor the temperature dependence of the phonon frequencies, crucial to account for temperature driven second-order phase transitions. Several methods developed in the last years to account for anharmonicity in the non-perturbative regime and overcome such difficulties are briefly reviewed in this colloquium paper. In particular, the stochastic self-consistent harmonic approximation, a variational method that allows calculating vibrational properties in strongly anharmonic systems, is described in further detail. Applications of the latter method to superconducting palladium, platinum, and sulfur hydrides are discussed, where anharmonicity has a huge impact on their vibrational and superconducting properties.

Quantum thermal transport through anharmonic systems: A self-consistent approach

We propose a feasible and effective approach to study quantum thermal transport through anharmonic systems. The main idea is to obtain an {\it effective} harmonic Hamiltonian for the anharmonic system by applying the self-consistent phonon theory. Using the effective harmonic Hamiltonian we study thermal transport within the framework of nonequilibrium Green's function method using the celebrated Caroli formula. We corroborate our quantum self-consistent approach using the quantum master equation that can deal with anharmonicity exactly, but is limited to the weak system-bath coupling regime. Finally, in order demonstrate its strength we apply the quantum self-consistent approach to study thermal rectification in a weakly coupled two segment anharmonic system.

Semiclassical treatment of quantum propagation with nonlinear classical dynamics: A third-order thawed Gaussian approximation.

The time-dependent WKB approximation for a coherent state is expanded to third order around a guiding real trajectory, allowing for the novel treatment of nonlinearity in its semiclassical dynamics, which is generally present in all physical systems far from the classical regime. The result is a closed-form solution consisting of a linear combination of Airy functions and their derivatives multiplied by an exponential. The expression's ability to capture nonlinearity is demonstrated by examining the autocorrelation of initial coherent states in anharmonic systems with few to many contributing periodic orbits. Its accuracy is compared to the quadratic expansion and found to be superior in regimes of ℏ where the curvature begins to be significant, as expected. Moreover, the expression is shown to be a real-trajectory uniformization over two coalescing saddle points that are emblematic of significant curvature. This extends real-trajectory time-dependent wave-packet semiclassical methods to highly anharmonic systems for the first time and establishes their regime of validity.

Thermodynamics of Anharmonic Systems: Uncoupled Mode Approximations for Molecules.

The partition functions, heat capacities, entropies, and enthalpies of selected molecules were calculated using uncoupled mode (UM) approximations, where the full-dimensional potential energy surface for internal motions was modeled as a sum of independent one-dimensional potentials for each mode. The computational cost of such approaches scales the same with molecular size as standard harmonic oscillator vibrational analysis using harmonic frequencies (HO(hf)). To compute thermodynamic properties, a computational protocol for obtaining the energy levels of each mode was established. The accuracy of the UM approximation depends strongly on how the one-dimensional potentials of each modes are defined. If the potentials are determined by the energy as a function of displacement along each normal mode (UM-N), the accuracies of the calculated thermodynamic properties are not significantly improved versus the HO(hf) model. Significant improvements can be achieved by constructing potentials for internal rotations and vibrations using the energy surfaces along the torsional coordinates and the remaining vibrational normal modes, respectively (UM-VT). For hydrogen peroxide and its isotopologs at 300 K, UM-VT captures more than 70% of the partition functions on average. By contrast, the HO(hf) model and UM-N can capture no more than 50%. For a selected test set of C2 to C8 linear and branched alkanes and species with different moieties, the enthalpies calculated using the HO(hf) model, UM-N, and UM-VT are all quite accurate comparing with reference values though the RMS errors of the HO model and UM-N are slightly higher than UM-VT. However, the accuracies in entropy calculations differ significantly between these three models. For the same test set, the RMS error of the standard entropies calculated by UM-VT is 2.18 cal mol(-1) K(-1) at 1000 K. By contrast, the RMS error obtained using the HO model and UM-N are 6.42 and 5.73 cal mol(-1) K(-1), respectively. For a test set composed of nine alkanes ranging from C5 to C8, the heat capacities calculated with the UM-VT model agree with the experimental values to within a RMS error of 0.78 cal mol(-1) K(-1), which is less than one-third of the RMS error of the HO(hf) (2.69 cal mol(-1) K(-1)) and UM-N (2.41 cal mol(-1) K(-1)) models.

Proton dynamics in hydrogen-bonded systems

ABSTRACTThe proton mean kinetic energy Ke(H) in various systems was calculated between 5 and 320 K using a semi-empirical (SE) approach. The SE calculation relies on the harmonic approximation and decoupling between the various modes, where the input data of the internal and external frequencies were taken from inelastic neutron scattering (INS) and IR/Raman experiments. The studied systems included ordinary H2O phases, water of crystallisation in sulphate salts, adsorbed water and water confined in various samples of pore dimensions less than 20 Å. These included some zeolites, periodic mesoporous organosilicas (PMOs), beryl, Bikitaite and single- and double-wall carbon nanotubes. All SE calculated Ke(H) values were close to that of pure ice/liquid water, for which a good agreement was found with the deep inelastic neutron scattering (DINS) measurements. However, for water in Beryl at 5 K and ice in carbon nanotubes, at 170 K, large deviations from DINS results were found. Some insight into this problem may be gained by comparing those deviations with recently studied anharmonic systems involving proton double well potentials: Rb3H(SO4)2 and KH2PO4, where an excellent agreement was obtained between SE calculations and DINS measurements.

Vibrational spectroscopy via the Caldeira-Leggett model with anharmonic system potentials.

The Caldeira-Leggett (CL) model, which describes a system bi-linearly coupled to a harmonic bath, has enjoyed popularity in condensed phase spectroscopy owing to its utmost simplicity. However, the applicability of the model to cases with anharmonic system potentials, as it is required for the description of realistic systems in solution, is questionable due to the presence of the invertibility problem [F. Gottwald et al., J. Phys. Chem. Lett. 6, 2722 (2015)] unless the system itself resembles the CL model form. This might well be the case at surfaces or in the solid regime, which we here confirm for a particular example of an iodine molecule in the atomic argon environment under high pressure. For this purpose we extend the recently proposed Fourier method for parameterizing linear generalized Langevin dynamics [F. Gottwald et al., J. Chem. Phys. 142, 244110 (2015)] to the non-linear case based on the CL model and perform an extensive error analysis. In order to judge on the applicability of this model in advance, we give practical empirical criteria and discuss the effect of the potential renormalization term. The obtained results provide evidence that the CL model can be used for describing a potentially broad class of systems.

Underlying mechanisms for normal heat transport in one-dimensional anharmonic oscillator systems with a double-well interparticle interaction

Previous studies have suggested a crossover from superdiffusive to normal heat transport in one-dimensional (1D) anharmonic oscillator systems with a double-well type interatomic interaction like , when the system temperature is varied. In order to better understand this unusual manner of thermal transport, here we perform a direct dynamics simulation to examine how the spreading processes of the three physical quantities, i.e. the heat, the total energy and the momentum, would depend on temperature. We find three main points that are worth noting. (i) The crossover from superdiffusive to normal heat transport is well verified from a new perspective of heat spread. (ii) The spreading of the total energy is found to be very distinct from heat diffusion, especially under some temperature regimes, energy is strongly localized, while heat can be superdiffusive. So one should take care to derive a general connection between the heat conduction and energy diffusion. (iii) In a narrow range of temperatures, the spreading of momentum implies clear unusual non-ballistic behaviors; however, such unusual transport of momentum cannot be directly related to the normal transport of heat. An analysis of phonon spectra suggests that one should also take the effects of phonon softening into account. All of these results may provide insights into establishing the connection between the macroscopic heat transport and the underlying dynamics in 1D systems.