Would You Choose to be Happy? Tradeoffs Between Happiness and the Other Dimensions of Life in a Large Population Survey

A large literature documents the correlates and causes of subjective well-being, or happiness. But few studies have investigated whether people choose happiness. Is happiness all that people want from life, or are they willing to sacrifice it for other attributes, such as income and health? Tackling this question has largely been the preserve of philosophers. In this article, we find out just how much happiness matters to ordinary citizens. Our sample consists of nearly 13,000 members of the UK and US general populations. We ask them to choose between, and make judgments over, lives that are high (or low) in different types of happiness and low (or high) in income, physical health, family, career success, or education. We find that people by and large choose the life that is highest in happiness but health is by far the most important other concern, with considerable numbers of people choosing to be healthy rather than happy. We discuss some possible reasons for this preference

Jacobi-Predictor-Corrector Approach for the Fractional Ordinary Differential Equations

We present a novel numerical method, called {\tt Jacobi-predictor-corrector approach}, for the numerical solution of fractional ordinary differential equations based on the polynomial interpolation and the Gauss-Lobatto quadrature w.r.t. the Jacobi-weight function $\omega(s)=(1-s)^{\alpha-1}(1+s)^0$. This method has the computational cost O(N) and the convergent order $IN$, where $N$ and $IN$ are, respectively, the total computational steps and the number of used interpolating points. The detailed error analysis is performed, and the extensive numerical experiments confirm the theoretical results and show the robustness of this method.Comment: 24 pages, 5 figure

Spectral Semi-discretisations of Weakly Non-linear Wave Equations over Long Times

The long-time behaviour of spectral semi-discretisations of weakly non-linear wave equations is analysed. It is shown that the harmonic actions are approximately conserved for the semi-discretised system as well. This permits to prove that the energy of the wave equation along the interpolated semi-discrete solution remains well conserved over long times and close to the Hamiltonian of the semi-discrete equation. Although the momentum is no longer an exact invariant of the semi-discretisation, it is shown to be approximately conserved. All these results are obtained with the technique of modulated Fourier expansion

Fourier methods for the perturbed harmonic oscillator in linear and nonlinear Schr\"odinger equations

We consider the numerical integration of the Gross-Pitaevskii equation with a potential trap given by a time-dependent harmonic potential or a small perturbation thereof. Splitting methods are frequently used with Fourier techniques since the system can be split into the kinetic and remaining part, and each part can be solved efficiently using Fast Fourier Transforms. To split the system into the quantum harmonic oscillator problem and the remaining part allows to get higher accuracies in many cases, but it requires to change between Hermite basis functions and the coordinate space, and this is not efficient for time-dependent frequencies or strong nonlinearities. We show how to build new methods which combine the advantages of using Fourier methods while solving the timedependent harmonic oscillator exactly (or with a high accuracy by using a Magnus integrator and an appropriate decomposition).Comment: 12 pages of RevTex4-1, 8 figures; substantially revised and extended versio

Double Time Window Targeting Technique: Real time DMRG dynamics in the PPP model

We present a generalized adaptive time-dependent density matrix renormalization group (DMRG) scheme, called the {\it double time window targeting} (DTWT) technique, which gives accurate results with nominal computational resources, within reasonable computational time. This procedure originates from the amalgamation of the features of pace keeping DMRG algorithm, first proposed by Luo {\it et. al}, [Phys.Rev. Lett. {\bf 91}, 049701 (2003)], and the time-step targeting (TST) algorithm by Feiguin and White [Phys. Rev. B {\bf 72}, 020404 (2005)]. Using the DTWT technique, we study the phenomena of spin-charge separation in conjugated polymers (materials for molecular electronics and spintronics), which have long-range electron-electron interactions and belong to the class of strongly correlated low-dimensional many-body systems. The issue of real time dynamics within the Pariser-Parr-Pople (PPP) model which includes long-range electron correlations has not been addressed in the literature so far. The present study on PPP chains has revealed that, (i) long-range electron correlations enable both the charge and spin degree of freedom of the electron, to propagate faster in the PPP model compared to Hubbard model, (ii) for standard parameters of the PPP model as applied to conjugated polymers, the charge velocity is almost twice that of the spin velocity and, (iii) the simplistic interpretation of long-range correlations by merely renormalizing the {\it U} value of the Hubbard model fails to explain the dynamics of doped holes/electrons in the PPP model.Comment: Final (published) version; 39 pages, 13 figures, 1 table; 2 new references adde

On the applicability of constrained symplectic integrators in general relativity

The purpose of this note is to point out that a naive application of symplectic integration schemes for Hamiltonian systems with constraints such as SHAKE or RATTLE which preserve holonomic constraints encounters difficulties when applied to the numerical treatment of the equations of general relativity.Comment: 13 pages, change the title to be more descriptive, typos corrected, added referenc

Open Boundaries for the Nonlinear Schrodinger Equation

We present a new algorithm, the Time Dependent Phase Space Filter (TDPSF) which is used to solve time dependent Nonlinear Schrodinger Equations (NLS). The algorithm consists of solving the NLS on a box with periodic boundary conditions (by any algorithm). Periodically in time we decompose the solution into a family of coherent states. Coherent states which are outgoing are deleted, while those which are not are kept, reducing the problem of reflected (wrapped) waves. Numerical results are given, and rigorous error estimates are described. The TDPSF is compatible with spectral methods for solving the interior problem. The TDPSF also fails gracefully, in the sense that the algorithm notifies the user when the result is incorrect. We are aware of no other method with this capability.Comment: 21 pages, 4 figure

Exponential propagators for the Schrödinger equation with a time-dependent potential

[EN] We consider the numerical integration of the Schrodinger equation with a time-dependent Hamiltonian given as the sum of the kinetic energy and a time-dependent potential. Commutator-free (CF) propagators are exponential propagators that have shown to be highly efficient for general time-dependent Hamiltonians. We propose new CF propagators that are tailored for Hamiltonians of the said structure, showing a considerably improved performance. We obtain new fourth-and sixth-order CF propagators as well as a novel sixth-order propagator that incorporates a double commutator that only depends on coordinates, so this term can be considered as cost-free. The algorithms require the computation of the action of exponentials on a vector similar to the well-known exponential midpoint propagator, and this is carried out using the Lanczos method. We illustrate the performance of the new methods on several numerical examples. Published by AIP Publishing.We wish to acknowledge Fernando Casas for his help in the construction of the methods Upsilon3[6]. The authors acknowledge Ministerio de Economia y Competitividad (Spain) for financial support through Project No. MTM2016-77660-P (AEI/FEDER, UE). Additionally, Kopylov has been partly supported by Grant No. GRISOLIA/2015/A/137 from the Generalitat Valenciana.Bader, P.; Kopylov, N.; Blanes Zamora, S. (2018). Exponential propagators for the Schrödinger equation with a time-dependent potential. The Journal of Chemical Physics. 149(24):1-7. https://doi.org/10.1063/1.5036838S1714924Bader, P., Iserles, A., Kropielnicka, K., & Singh, P. (2014). Effective Approximation for the Semiclassical Schrödinger Equation. Foundations of Computational Mathematics, 14(4), 689-720. doi:10.1007/s10208-013-9182-8Bader, P., Iserles, A., Kropielnicka, K., & Singh, P. (2016). Efficient methods for linear Schrödinger equation in the semiclassical regime with time-dependent potential. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 472(2193), 20150733. doi:10.1098/rspa.2015.0733Lubich, C. (2008). From Quantum to Classical Molecular Dynamics: Reduced Models and Numerical Analysis. doi:10.4171/067Feit, M. ., Fleck, J. ., & Steiger, A. (1982). Solution of the Schrödinger equation by a spectral method. Journal of Computational Physics, 47(3), 412-433. doi:10.1016/0021-9991(82)90091-2Tremblay, J. C., & Carrington, T. (2004). Using preconditioned adaptive step size Runge-Kutta methods for solving the time-dependent Schrödinger equation. The Journal of Chemical Physics, 121(23), 11535-11541. doi:10.1063/1.1814103Sanz‐Serna, J. M., & Portillo, A. (1996). Classical numerical integrators for wave‐packet dynamics. The Journal of Chemical Physics, 104(6), 2349-2355. doi:10.1063/1.470930Kormann, K., Holmgren, S., & Karlsson, H. O. (2008). Accurate time propagation for the Schrödinger equation with an explicitly time-dependent Hamiltonian. The Journal of Chemical Physics, 128(18), 184101. doi:10.1063/1.2916581Peskin, U., Kosloff, R., & Moiseyev, N. (1994). The solution of the time dependent Schrödinger equation by the (t,t’) method: The use of global polynomial propagators for time dependent Hamiltonians. The Journal of Chemical Physics, 100(12), 8849-8855. doi:10.1063/1.466739Castro, A., Marques, M. A. L., & Rubio, A. (2004). Propagators for the time-dependent Kohn–Sham equations. The Journal of Chemical Physics, 121(8), 3425-3433. doi:10.1063/1.1774980Blanes, S., Casas, F., & Murua, A. (2017). Symplectic time-average propagators for the Schrödinger equation with a time-dependent Hamiltonian. The Journal of Chemical Physics, 146(11), 114109. doi:10.1063/1.4978410Blanes, S., Casas, F., & Murua, A. (2015). An efficient algorithm based on splitting for the time integration of the Schrödinger equation. Journal of Computational Physics, 303, 396-412. doi:10.1016/j.jcp.2015.09.047Gray, S. K., & Verosky, J. M. (1994). Classical Hamiltonian structures in wave packet dynamics. The Journal of Chemical Physics, 100(7), 5011-5022. doi:10.1063/1.467219McLachlan, R. I., & Quispel, G. R. W. (2002). Splitting methods. Acta Numerica, 11, 341-434. doi:10.1017/s0962492902000053Neuhauser, C., & Thalhammer, M. (2009). On the convergence of splitting methods for linear evolutionary Schrödinger equations involving an unbounded potential. BIT Numerical Mathematics, 49(1), 199-215. doi:10.1007/s10543-009-0215-2Thalhammer, M. (2008). High-Order Exponential Operator Splitting Methods for Time-Dependent Schrödinger Equations. SIAM Journal on Numerical Analysis, 46(4), 2022-2038. doi:10.1137/060674636Thalhammer, M. (2012). Convergence Analysis of High-Order Time-Splitting Pseudospectral Methods for Nonlinear Schrödinger Equations. SIAM Journal on Numerical Analysis, 50(6), 3231-3258. doi:10.1137/120866373Gray, S. K., & Manolopoulos, D. E. (1996). Symplectic integrators tailored to the time‐dependent Schrödinger equation. The Journal of Chemical Physics, 104(18), 7099-7112. doi:10.1063/1.471428Sanz-Serna, J. M., & Calvo, M. P. (1994). Numerical Hamiltonian Problems. doi:10.1007/978-1-4899-3093-4Saad, Y. (1992). Analysis of Some Krylov Subspace Approximations to the Matrix Exponential Operator. SIAM Journal on Numerical Analysis, 29(1), 209-228. doi:10.1137/0729014Park, T. J., & Light, J. C. (1986). Unitary quantum time evolution by iterative Lanczos reduction. The Journal of Chemical Physics, 85(10), 5870-5876. doi:10.1063/1.451548Blanes, S., & Moan, P. C. (2000). Splitting methods for the time-dependent Schrödinger equation. Physics Letters A, 265(1-2), 35-42. doi:10.1016/s0375-9601(99)00866-xMagnus, W. (1954). On the exponential solution of differential equations for a linear operator. Communications on Pure and Applied Mathematics, 7(4), 649-673. doi:10.1002/cpa.3160070404Blanes, S., Casas, F., Oteo, J. A., & Ros, J. (2009). The Magnus expansion and some of its applications. Physics Reports, 470(5-6), 151-238. doi:10.1016/j.physrep.2008.11.001Blanes, S., & Moan, P. C. (2006). Fourth- and sixth-order commutator-free Magnus integrators for linear and non-linear dynamical systems. Applied Numerical Mathematics, 56(12), 1519-1537. doi:10.1016/j.apnum.2005.11.004Thalhammer, M. (2006). A fourth-order commutator-free exponential integrator for nonautonomous differential equations. SIAM Journal on Numerical Analysis, 44(2), 851-864. doi:10.1137/05063042Alvermann, A., & Fehske, H. (2011). High-order commutator-free exponential time-propagation of driven quantum systems. Journal of Computational Physics, 230(15), 5930-5956. doi:10.1016/j.jcp.2011.04.006Auer, N., Einkemmer, L., Kandolf, P., & Ostermann, A. (2018). Magnus integrators on multicore CPUs and GPUs. Computer Physics Communications, 228, 115-122. doi:10.1016/j.cpc.2018.02.019Munthe–Kaas, H., & Owren, B. (1999). Computations in a free Lie algebra. Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 357(1754), 957-981. doi:10.1098/rsta.1999.0361Bader, P., Blanes, S., Ponsoda, E., & Seydaoğlu, M. (2017). Symplectic integrators for the matrix Hill equation. Journal of Computational and Applied Mathematics, 316, 47-59. doi:10.1016/j.cam.2016.09.041P.V. Koseleff, “Formal calculus for Lie methods in Hamiltonian mechanics,” Ph.D. thesis, Lawrence Berkeley Laboratory, 1994.Chin, S. A. (1997). Symplectic integrators from composite operator factorizations. Physics Letters A, 226(6), 344-348. doi:10.1016/s0375-9601(97)00003-0Omelyan, I. P., Mryglod, I. M., & Folk, R. (2002). Construction of high-order force-gradient algorithms for integration of motion in classical and quantum systems. Physical Review E, 66(2). doi:10.1103/physreve.66.026701Walker, R. B., & Preston, R. K. (1977). Quantum versus classical dynamics in the treatment of multiple photon excitation of the anharmonic oscillator. The Journal of Chemical Physics, 67(5), 2017. doi:10.1063/1.435085Iserles, A., Kropielnicka, K., & Singh, P. (2018). Magnus--Lanczos Methods with Simplified Commutators for the Schrödinger Equation with a Time-Dependent Potential. SIAM Journal on Numerical Analysis, 56(3), 1547-1569. doi:10.1137/17m114983