BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Siu A. Chin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This work shows that, in one dimension, due to its topology, a closed-loop product of short-time propagators is always positive, despite the fact that each anti-symmetric free fermion propagator can be of either sign.
This work gives a Lie operator derivation of various Boris solvers via a detailed study of trajectory errors in a constant magnetic field. These errors in the gyrocenter location and the gyroradius are the foundational basis for why Boris solvers existed, independent of any finite-difference schemes. This work shows that there are two distinct ways of eliminating these errors so that the trajectory of a charged particle in a constant magnetic field is exactly on the cyclotron orbit. One way reproduces the known second-order symmetric Boris solver. The other yields a previously unknown, but also on-orbit solver, not derivable from finite-difference schemes. By revisiting some historical calculations, it is found that many publications do not distinguish the poorly behaved leap-frog Boris solver from the symmetric second-order Boris algorithm. This symmetric second-order Boris solver's trajectory is much more accurate and remains close to the exact orbit in a combined $nonuniform$ electric and magnetic field at time-steps greater than the cyclotron period. Finally, this operator formalism showed that Buneman's cycloid fitting scheme is mathematically identical to Boris' on-orbit solver and that Boris' E-B splitting is unnecessary.
We investigate higher-order breathers of the cubic nonlinear Schrödinger equation on an elliptic background. We find that, beyond first-order, any arbitrarily constructed breather is a single-peaked solitary wave on a disordered background. These "quasi-rogue waves" are also common on periodic backgrounds. We assume the higher-order breather is constructed out of constituent first-order breathers with commensurate periods (i.e., higher-order harmonic waves). In that case, one obtains "quasi-periodic" breathers with distorted side-peaks. Fully periodic breathers are obtained when their wavenumbers are harmonic multiples of the background and each other. They are truly rare, requiring finely-tuned parameters. Thus, on a periodic background, we arrive at the paradoxical conclusion that the apparent higher-order rogue waves are rather common, while the truly periodic breathers are exceedingly rare.
It is well known that the use of the primitive second-order propagator in Path Integral Monte Carlo calculations of many-fermion systems leads to the sign problem. In this work, we show that by using the similarity-transformed Fokker-Planck propagator, it is possible to solve for the ground state of a large quantum dot, with up to 100 polarized electrons, without solving the sign problem. These similarity-transformed propagators naturally produce rotational symmetry-breaking ground state wave functions previously used in the study of quantum dots and quantum Hall effects. However, instead of localizing the electrons at positions which {\it minimize} the potential energy, this derivation shows that they should be located at positions which {\it maximize} the bosonic ground state wave function. Further improvements in the energy can be obtained by using these as initial wave functions in a Ground State Path-Integral Monte Carlo calculation with second and fourth-order propagators.
It is well known that Akhmediev breathers of the nonlinear cubic Schrodinger equation can be superposed nonlinearly via the Darboux transformation to yield breathers of higher order. Surprisingly, we find that the peak height of each Akhmediev breather only adds {\it linearly} to form the peak height of the final breather. Using this new peak-height formula, we show that at any given periodicity, there exist a unique high-order breather of maximal intensity. Moreover, these high-order breathers form a continuous hierarchy, growing in intensity with increasing periodicity. For any such higher-order breather, a simple initial wave function can be extracted from the Darboux transformation to dynamically generate that breather from the nonlinear Schrodinger equation.
By invoking Bogoliubov's spectrum, we show that for the nonlinear Schrodinger equation, the modulation instability (MI) of its n = 1 Fourier mode on a finite background automatically triggers a further cascading instability, forcing all the higher modes to grow exponentially in locked-step with the n = 1 mode. This fundamental insight, the enslavement of all higher modes to the n = 1 mode, explains the formation of a triangular-shaped spectrum which generates the Akhmediev breather, predicts its formation time analytically from the initial modulation amplitude, and shows that the Fermi-Pasta-Ulam (FPU) recurrence is just a matter of energy conservation with a period twice the breather's formation time. For higher order MI with more than one initial unstable modes, while most evolutions are expected to be chaotic, we show that it is possible to have isolated cases of "super-recurrence", where the FPU period is much longer than that of a single unstable mode.
The conventional second-order Path Integral Monte Carlo method is plagued with the sign problem in solving many-fermion systems. This is due to the large number of anti-symmetric free fermion propagators that are needed to extract the ground state wave function at large imaginary time. In this work, we show that optimized fourth-order Path Integral Monte Carlo methods, which use no more than 5 free-fermion propagators, can yield accurate quantum dot energies for up to 20 polarized electrons with the use of the Hamiltonian energy estimator.
An elementary proof of Bertrand's theorem is given by examining the radial orbit equation, without needing to solve complicated equations or integrals.
Conventional finite-difference schemes for solving partial differential equations are based on approximating derivatives by finite-differences. In this work, an alternative theory is proposed which view finite-difference schemes as systematic ways of matching up to the operator solution of the partial differential equation. By completely abandon the idea of approximating derivatives directly, the theory provides a unified description of explicit finite-difference schemes for solving a general linear partial differential equation with constant coefficients to any time-marching order. As a result, the stability of the first-order algorithm for an entire class of linear equations can be determined all at once. Because the method is based on solution-matching, it can also be used to derive any order schemes for solving the general nonlinear advection equation.
Prior to the recent development of symplectic integrators, the time-stepping operator $\e^{h(A+B)}$ was routinely decomposed into a sum of products of $\e^{h A}$ and $\e^{hB}$ in the study of hyperbolic partial differential equations. In the context of solving Hamiltonian dynamics, we show that such a decomposition give rises to both {\it even} and {\it odd} order Runge-Kutta and Nyström integrators. By use of Suzuki's forward-time derivative operator to enforce the time-ordered exponential, we show that the same decomposition can be used to solve non-autonomous equations. In particular, odd order algorithms are derived on the basis of a highly non-trivial {\it time-asymmetric} kernel. Such an operator approach provides a general and unified basis for understanding structure non-preserving algorithms and is especially useful in deriving very high-order algorithms via {\it analytical} extrapolations. In this work, algorithms up to the 100th order are tested by integrating the ground state wave function of the hydrogen atom. For such a singular Coulomb problem, the multi-product expansion showed uniform convergence and is free of poles usually associated with structure-preserving methods. Other examples are also discussed.
We present a new class of high-order imaginary time propagators for path-integral Monte Carlo simulations by subtracting lower order propagators. By requiring all terms of the extrapolated propagator be sampled uniformly, the subtraction only affects the potential part of the path integral. The negligible violation of positivity of the resulting path integral at small time steps has no discernable affect on the accuracy of our method. Thus in principle arbitrarily high order algorithms can be devised for path-integral Monte Carlo simulations. We verify this claim is by showing that fourth, sixth, and eighth order convergence can indeed be achieved in solving for the ground state of strongly interacting quantum many-body systems such as bulk liquid $^4$He.
Since the kinetic and the potential energy term of the real time nonlinear Schrödinger equation can each be solved exactly, the entire equation can be solved to any order via splitting algorithms. We verified the fourth-order convergence of some well known algorithms by solving the Gross-Pitaevskii equation numerically. All such splitting algorithms suffer from a latent numerical instability even when the total energy is very well conserved. A detail error analysis reveals that the noise, or elementary excitations of the nonlinear Schrödinger, obeys the Bogoliubov spectrum and the instability is due to the exponential growth of high wave number noises caused by the splitting process. For a continuum wave function, this instability is unavoidable no matter how small the time step. For a discrete wave function, the instability can be avoided only for $\dt k_{max}^2{<\atop\sim}2 π$, where $k_{max}=π/Δx$.
By implementing the exact density matrix for the rotating anisotropic harmonic trap, we derive a class of very fast and accurate fourth order algorithms for evolving the Gross-Pitaevskii equation in imaginary time. Such fourth order algorithms are possible only with the use of {\it forward}, positive time step factorization schemes. These fourth order algorithms converge at time-step sizes an order-of-magnitude larger than conventional second order algorithms. Our use of time-dependent factorization schemes provides a systematic way of devising algorithms for solving this type of nonlinear equations.
We show that the method of splitting the operator ${\rm e}^{ε(T+V)}$ to fourth order with purely positive coefficients produces excellent algorithms for solving the time-dependent Schrödinger equation. These algorithms require knowing the potential and the gradient of the potential. One 4th order algorithm only requires four Fast Fourier Transformations per iteration. In a one dimensional scattering problem, the 4th order error coefficients of these new algorithms are roughly 500 times smaller than fourth order algorithms with negative coefficient, such as those based on the traditional Ruth-Forest symplectic integrator. These algorithms can produce converged results of conventional second or fourth order algorithms using time steps 5 to 10 times as large. Iterating these positive coefficient algorithms to 6th order also produced better converged algorithms than iterating the Ruth-Forest algorithm to 6th order or using Yoshida's 6th order algorithm A directly.