BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ander Murua

Ander Murua contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAphysics.comp-phastro-ph.EPmath-phmath.DSmath.MPmath.RAquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2016arXiv

Vibrational resonance: a study with high-order word-series averaging

We study a model problem describing vibrational resonance by means of a high-order averaging technique based on so-called word series. With the tech- nique applied here, the tasks of constructing the averaged system and the associ- ated change of variables are divided into two parts. It is first necessary to build recursively a set of so-called word basis functions and, after that, all the required manipulations involve only scalar coefficients that are computed by means of sim- ple recursions. As distinct from the situation with other approaches, with word- series, high-order averaged systems may be derived without having to compute the associated change of variables. In the system considered here, the construction of high-order averaged systems makes it possible to obtain very precise approxima- tions to the true dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formal series and numerical integrators: some history and some new techniques

This paper provides a brief history of B-series and the associated Butcher group and presents the new theory of word series and extended word series. B-series (Hairer and Wanner 1976) are formal series of functions parameterized by rooted trees. They greatly simplify the study of Runge-Kutta schemes and other numerical integrators. We examine the problems that led to the introduction of B-series and survey a number of more recent developments, including applications outside numerical mathematics. Word series (series of functions parameterized by words from an alphabet) provide in some cases a very convenient alternative to B-series. Associated with word series is a group G of coefficients with a composition rule simpler than the corresponding rule in the Butcher group. From a more mathematical point of view, integrators, like Runge-Kutta schemes, that are affine equivariant are represented by elements of the Butcher group, integrators that are equivariant with respect to arbitrary changes of variables are represented by elements of the word group G.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New families of symplectic splitting methods for numerical integration in dynamical astronomy

We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We derive in a systematic way an independent set of necessary and sufficient conditions to be satisfied by the coefficients of splitting methods to achieve a prescribed order of accuracy. Splitting methods satisfying such (generalized) order conditions are appropriate in particular for the numerical simulation of the Solar System described in Jacobi coordinates. We show that, when using Poincaré Heliocentric coordinates, the same order of accuracy may be obtained by imposing an additional polynomial equation on the coefficients of the splitting method. We construct several splitting methods appropriate for each of the two sets of coordinates by solving the corresponding systems of polynomial equations and finding the optimal solutions. The experiments reported here indicate that the efficiency of our new schemes is clearly superior to previous integrators when high accuracy is required.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Word series for dynamical systems and their numerical integrators

We study word series and extended word series, classes of formal series for the analysis of some dynamical systems and their discretizations. These series are similar to but more compact than B-series. They may be composed among themselves by means of a simple rule. While word series have appeared before in the literature, extended word series are introduced in this paper. We exemplify the use of extended word series by studying the reduction to normal form and averaging of some perturbed integrable problems. We also provide a detailed analysis of the behaviour of splitting numerical methods for those problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient computation of the Zassenhaus formula

A new recursive procedure to compute the Zassenhaus formula up to high order is presented, providing each exponent in the factorization directly as a linear combination of independent commutators and thus containing the minimum number of terms. The recursion can be easily implemented in a symbolic algebra package and requires much less computational effort, both in time and memory resources, than previous algorithms. In addition, by bounding appropriately each term in the recursion, it is possible to get a larger convergence domain of the Zassenhaus formula when it is formulated in a Banach algebra.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High precision Symplectic Integrators for the Solar System

Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These comparisons are made in Jacobi and Heliocentric coordinates and the implementation of the algorithms is fully detailed for practical use. We conclude that high order integrators should be privileged, with a preference for the new $(10,6,4)$ method of (Blanes et al., 2012)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Error analysis of splitting methods for the time dependent Schrodinger equation

A typical procedure to integrate numerically the time dependent Schrö\-din\-ger equation involves two stages. In the first one carries out a space discretization of the continuous problem. This results in the linear system of differential equations $i du/dt = H u$, where $H$ is a real symmetric matrix, whose solution with initial value $u(0) = u_0 \in \mathbb{C}^N$ is given by $u(t) = \e^{-i t H} u_0$. Usually, this exponential matrix is expensive to evaluate, so that time stepping methods to construct approximations to $u$ from time $t_n$ to $t_{n+1}$ are considered in the second phase of the procedure. Among them, schemes involving multiplications of the matrix $H$ with vectors, such as Lanczos and Chebyshev methods, are particularly efficient. In this work we consider a particular class of splitting methods which also involves only products $Hu$. We carry out an error analysis of these integrators and propose a strategy which allows us to construct different splitting symplectic methods of different order (even of order zero) possessing a large stability interval that can be adapted to different space regularity conditions and different accuracy ranges of the spatial discretization. The validity of the procedure and the performance of the resulting schemes are illustrated on several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimized high-order splitting methods for some classes of parabolic equations

We are concerned with the numerical solution obtained by splitting methods of certain parabolic partial differential equations. Splitting schemes of order higher than two with real coefficients necessarily involve negative coefficients. It has been demonstrated that this second-order barrier can be overcome by using splitting methods with complex-valued coefficients (with positive real parts). In this way, methods of orders 3 to 14 by using the Suzuki--Yoshida triple (and quadruple) jump composition procedure have been explicitly built. Here we reconsider this technique and show that it is inherently bounded to order 14 and clearly sub-optimal with respect to error constants. As an alternative, we solve directly the algebraic equations arising from the order conditions and construct methods of orders 6 and 8 that are the most accurate ones available at present time, even when low accuracies are desired. We also show that, in the general case, 14 is not an order barrier for splitting methods with complex coefficients with positive real part by building explicitly a method of order 16 as a composition of methods of order 8.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Splitting methods with complex coefficients

Splitting methods for the numerical integration of differential equations of order greater than two involve necessarily negative coefficients. This order barrier can be overcome by considering complex coefficients with positive real part. In this work we review the composition technique used to construct methods of this class, propose new sixth-order integrators and analyze their main features on a pair of numerical examples, in particular how the errors are propagated along the evolution.

0 citations0 reviewsNo code linkedOpen access