New path integral simulation algorithms and their application to creep in the quantum sine-Gordon chain

A path integral simulation algorithm which includes a higher-order Trotter approximation (HOA)is analyzed and compared to an approach which includes the correct quantum mechanical pair interaction (effective Propagator (EPr)). It is found that the HOA algorithmconverges to the quantum limit with increasing Trotter number P as P^{-4}, while the EPr algorithm converges as P^{-2}.The convergence rate of the HOA algorithm is analyzed for various physical systemssuch as a harmonic chain,a particle in a double-well potential, gaseous argon, gaseous helium and crystalline argon. A new expression for the estimator for the pair correlation function in the HOA algorithm is derived. A new path integral algorithm, the hybrid algorithm, is developed.It combines an exact treatment of the quadratic part of the Hamiltonian and thehigher-order Trotter expansion techniques.For the discrete quantum sine-Gordon chain (DQSGC), it is shown that this algorithm works more efficiently than all other improved path integral algorithms discussed in this work. The new simulation techniques developed in this work allow the analysis of theDQSGC and disordered model systems in the highly quantum mechanical regime using path integral molecular dynamics (PIMD)and adiabatic centroid path integral molecular dynamics (ACPIMD).The ground state phonon dispersion relation is calculated for the DQSGC by the ACPIMD method.It is found that the excitation gap at zero wave vector is reduced by quantum fluctuations. Two different phases exist: One phase with a finite excitation gap at zero wave vector, and a gapless phase where the excitation gap vanishes.The reaction of the DQSGC to an external driving force is analyzed at T=0.In the gapless phase the system creeps if a small force is applied, and in the phase with a gap the system is pinned. At a critical force, the systems undergo a depinning transition in both phases and flow is induced. The analysis of the DQSGC is extended to models with disordered substrate potentials. Three different cases are analyzed: Disordered substrate potentials with roughness exponent H=0, H=1/2,and a model with disordered bond length. For all models, the ground state phonon dispersion relation is calculated.

Variational Markov Chain Monte Carlo for Bayesian smoothing of non-linear niffusions

In this paper we develop set of novel Markov chain Monte Carlo algorithms for Bayesian smoothing of partially observed non-linear diffusion processes. The sampling algorithms developed herein use a deterministic approximation to the posterior distribution over paths as the proposal distribution for a mixture of an independence and a random walk sampler. The approximating distribution is sampled by simulating an optimized time-dependent linear diffusion process derived from the recently developed variational Gaussian process approximation method. Flexible blocking strategies are introduced to further improve mixing, and thus the efficiency, of the sampling algorithms. The algorithms are tested on two diffusion processes: one with double-well potential drift and another with SINE drift. The new algorithm's accuracy and efficiency is compared with state-of-the-art hybrid Monte Carlo based path sampling. It is shown that in practical, finite sample, applications the algorithm is accurate except in the presence of large observation errors and low observation densities, which lead to a multi-modal structure in the posterior distribution over paths. More importantly, the variational approximation assisted sampling algorithm outperforms hybrid Monte Carlo in terms of computational efficiency, except when the diffusion process is densely observed with small errors in which case both algorithms are equally efficient.

Variational Markov chain Monte Carlo for Bayesian smoothing of non-linear diffusions

In this paper we develop set of novel Markov Chain Monte Carlo algorithms for Bayesian smoothing of partially observed non-linear diffusion processes. The sampling algorithms developed herein use a deterministic approximation to the posterior distribution over paths as the proposal distribution for a mixture of an independence and a random walk sampler. The approximating distribution is sampled by simulating an optimized time-dependent linear diffusion process derived from the recently developed variational Gaussian process approximation method. The novel diffusion bridge proposal derived from the variational approximation allows the use of a flexible blocking strategy that further improves mixing, and thus the efficiency, of the sampling algorithms. The algorithms are tested on two diffusion processes: one with double-well potential drift and another with SINE drift. The new algorithm's accuracy and efficiency is compared with state-of-the-art hybrid Monte Carlo based path sampling. It is shown that in practical, finite sample applications the algorithm is accurate except in the presence of large observation errors and low to a multi-modal structure in the posterior distribution over paths. More importantly, the variational approximation assisted sampling algorithm outperforms hybrid Monte Carlo in terms of computational efficiency, except when the diffusion process is densely observed with small errors in which case both algorithms are equally efficient. © 2011 Springer-Verlag.

Stable algorithm for event detection in event-driven particle dynamics : logical states

Acknowledgments The authors gratefully acknowledge the support of the German Research Foundation (DFG) through the Cluster of Excellence ‘Engineering of Advanced Materials’ at the University of Erlangen-Nuremberg and through Grant Po 472/25.

Influence of Band Non-Parabolicity on Few Ballistic Properties of III-V Quantum Wire Field Effect Transistors Under Strong Inversion

A simplified yet analytical approach on few ballistic properties of III-V quantum wire transistor has been presented by considering the band non-parabolicity of the electrons in accordance with Kane's energy band model using the Bohr-Sommerfeld's technique. The confinement of the electrons in the vertical and lateral directions are modeled by an infinite triangular and square well potentials respectively, giving rise to a two dimensional electron confinement. It has been shown that the quantum gate capacitance, the drain currents and the channel conductance in such systems are oscillatory functions of the applied gate and drain voltages at the strong inversion regime. The formation of subbands due to the electrical and structural quantization leads to the discreetness in the characteristics of such 1D ballistic transistors. A comparison has also been sought out between the self-consistent solution of the Poisson's-Schrodinger's equations using numerical techniques and analytical results using Bohr-Sommerfeld's method. The results as derived in this paper for all the energy band models gets simplified to the well known results under certain limiting conditions which forms the mathematical compatibility of our generalized theoretical formalism.

Barrier crossing in one and three dimensions by a long chain

We consider the Kramers problem for a long chain polymer trapped in a biased double-well potential. Initially the polymer is in the less stable well and it can escape from this well to the other well by the motion of its N beads across the barrier to attain the configuration having lower free energy. In one dimension we simulate the crossing and show that the results are in agreement with the kink mechanism suggested earlier. In three dimensions, it has not been possible to get an analytical `kink solution' for an arbitrary potential; however, one can assume the form of the solution of the nonlinear equation as a kink solution and then find a double-well potential in three dimensions. To verify the kink mechanism, simulations of the dynamics of a discrete Rouse polymer model in a double well in three dimensions are carried out. We find that the time of crossing is proportional to the chain length, which is in agreement with the results for the kink mechanism. The shape of the kink solution is also in agreement with the analytical solution in both one and three dimensions.

Aerodynamic pressure variation over SMA wire integrated morphing aerofoil

An analysis of the pressure variation over an aerofoil with integrated Shape Memory Alloy (SMA) wire is reported. A computational model based on finite elements and potential flow computation is proposed to obtain the deflections of the upper and the lower skins of the aerofoil subjected to aerodynamic pressure and hysteretic deformation of the SMA wire. The computational model couples a one-dimensional phenomenological constitutive model of SMA wire with the laminar incompressible aerodynamic pressure induced deformation of the aerofoil skins. The SMA wires are actuated by thermoelectric control system with auxiliary compensator feeding the piezoelectric stack actuators to adjust the hysteretic dynamics of the SMA wire. At each step of this coupled deformation process, the deflected/morphed shape of the aerofoil is d while recalculating to get the pressure distribution. Panel method based on incompressible and inviscid flow is employed for this purpose. The aerodynamic lift is then obtained from the pressure distributions. Numerical results on the variation of coefficient of pressure are reported.

Structural-modulation-driven low-temperature glassy behavior in SrRuO3

We present a critical investigation on the structural, magnetic, and magnetotransport properties of two sets of polycrystalline SrRuO3 samples with uniquely defined ferromagnetic transition temperatures. The ac magnetic susceptibility study exhibits the remarkable memory effect, a distinct characteristic of glassy behavior, at low temperatures. The transport study suggests a crossover from Fermi-liquid to non-Fermi-liquid behavior. Most strikingly, the temperature-dependent magnetoresistance exhibits two distinct dips (one around ferromagnetic ordering temperature and the other around 50 K), resembling a double-well potential in appearance. In addition, the temperature-dependent coercive field shows a plateau around 50 K. An attempt has been made to employ neutron diffraction to understand the genesis of such unusual low-temperature magnetic features. From the neutron-diffraction study, we find the evidence for changes in the unit-cell lattice parameters around 60 K and, thus, believe that the low-temperature anomalous magnetic response is closely intertwined to lattice-parameter change.

A pseudo-spin model for the wall of cytoskeletal microtubule

A pseudo-spin model is intended to describe the physical dynamics of unbound electrons in the wall of cytoskeletal microtubule (MT). Due to the inherent symmetry of the structure and the electric properties in the MT, one may treat it as a one-dimensional ferroelectric system, and describe the nonlinear dynamics of dimer electric dipoles in one protofilament of the MT by virtue of the double-well potential. Consequently, the physical problem has been mapped onto the pseudo-spin system, and the mean-field approximation has been taken to get some physical results.

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.

Qubit Residence Time Measurements with a Bose-Einstein Condensate

We show that an electrostatic qubit located near a Bose-Einstein condensate trapped in a symmetric double-well potential can be used to measure the duration the qubit has spent in one of its quantum states. The strong, medium, and weak measurement regimes are analyzed. The analogy between the residence and the traversal (tunnelling) times is highlighted.

Unusual decoherence in qubit measurements with a Bose-Einstein condensate

We consider an electrostatic qubit located near a Bose-Einstein condensate (BEC) of noninteracting bosons in a double-well potential, which is used for qubit measurements. Tracing out the BEC variables we obtain a simple analytical expression for the qubit's density matrix. The qubit's evolution exhibits a slow (proportional to 1/root t) damping of the qubit's coherence term, which however turns to be a Gaussian one in the case of static qubit. This is in contrast to the exponential damping produced by most classical detectors. The decoherence is, in general, incomplete and strongly depends on the initial state of the qubit.

Activity Enhancement of Tetrahexahedral Pd Nanocrystals by Bi Decoration towards Ethanol Electrooxidation in Alkaline Media

Tetrahexahedral Pd nanocrystals (THH Pd NCs) were prepared on a glassy carbon electrode using a programmed square-wave potential electrodeposition method, and modified by Bi adatoms with a range of coverages via the cyclic voltammetry method. The reactivity of the catalysts prepared towards ethanol electrooxidation reaction (EOR) was studied in alkaline medium at various temperatures and under other conditions that practical fuel cells operate. Significant activity enhancements were observed for the Bi-modified THH Pd NCs with an optimum Bi coverage (θBi) of around 0.68 being obtained. Furthermore, it was found that increasing temperature from 25 ºC to 60 ºC enhances the reactivity significantly. The general kinetics data of EOR on Bi-decorated and bare THH Pd NCs have also been obtained, from the activation energy calculated based on Arrhenius plots, and compared. At the optimum Bi coverage, an enhancement in the activity of almost 3 times was achieved, and the corresponding activation energy was found to be reduced significantly.

A statistical mechanical approach for the computation of the climatic response to general forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. In this paper we show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied with great success to analyze the climatic response to general forcings. The crucial value of the Ruelle theory lies in the fact that it allows to compute the response of the system in terms of expectation values of explicit and computable functions of the phase space averaged over the invariant measure of the unperturbed state. We choose as test bed a classical version of the Lorenz 96 model, which, in spite of its simplicity, has a well-recognized prototypical value as it is a spatially extended one-dimensional model and presents the basic ingredients, such as dissipation, advection and the presence of an external forcing, of the actual atmosphere. We recapitulate the main aspects of the general response theory and propose some new general results. We then analyze the frequency dependence of the response of both local and global observables to perturbations having localized as well as global spatial patterns. We derive analytically several properties of the corresponding susceptibilities, such as asymptotic behavior, validity of Kramers-Kronig relations, and sum rules, whose main ingredient is the causality principle. We show that all the coefficients of the leading asymptotic expansions as well as the integral constraints can be written as linear function of parameters that describe the unperturbed properties of the system, such as its average energy. Some newly obtained empirical closure equations for such parameters allow to define such properties as an explicit function of the unperturbed forcing parameter alone for a general class of chaotic Lorenz 96 models. We then verify the theoretical predictions from the outputs of the simulations up to a high degree of precision. The theory is used to explain differences in the response of local and global observables, to define the intensive properties of the system, which do not depend on the spatial resolution of the Lorenz 96 model, and to generalize the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable for finite as well as infinite time. Finally, we propose a simple yet general methodology to study general Climate Change problems on virtually any time scale by resorting to only well selected simulations, and by taking full advantage of ensemble methods. The specific case of globally averaged surface temperature response to a general pattern of change of the CO2 concentration is discussed. We believe that the proposed approach may constitute a mathematically rigorous and practically very effective way to approach the problem of climate sensitivity, climate prediction, and climate change from a radically new perspective.

Gas-Phase Nucleophilic and Elimination Reactions in Simple Alkyl Nitrates

There has been increasing interest in the gas-phase reactivity of alkyl nitrates because of their well-known applications as explosives and because of then role in atmospheric and in marine processes This manuscript describes an experimental study by FT-ICR techniques of the gas-phase reactions of OH(-) and F(-) with methyl and ethyl Innate For methyl nitrate, the main reaction channel is found to be an elimination process promoted by abstraction of an a proton from the methyl group. Nucleophilic displacement of nitrate anion through an S(N)2 process at the carbon center Is also found to he an important reaction channel with methyl nitrate In ethyl nitrate, Ruination of NO(3)(-) is greatly enhanced and this is attributed to the ease of an E2-type elimination process promoted by proton abstraction at the beta position of the ethyl group. Theoretical calculations at the MP2/6-311+G(3df,2p)//MP2/6-31+G(d) level of theory ale consistent with the relative importance of the reaction channels and suggest that these reactions proceed through a double well potential The calculations also predict that nucleophilic attack by OH(-) at the nitrogen center (Sn2@N) is energetically the rueful ad pathway but experiments with (18)OH(-) showed no evidence for this channel. Single-point calculations reveal a strong preference for approach to the emboli center and may explain the lack of reactivity at the nitrogen center. Calculations were also carried out or NH(2)(-) and SH(-) to establish the reactivity pattern to provide a better understanding of environmentally relevant nitrate esters.