Source author record

Aleksandar Belic

Aleksandar Belic appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

16 published item(s)

preprint2019arXiv

AEDGE: Atomic Experiment for Dark Matter and Gravity Exploration in Space

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

preprint2016arXiv

Hybrid OpenMP/MPI programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

We present hybrid OpenMP/MPI (Open Multi-Processing/Message Passing Interface) parallelized versions of earlier published C programs (D. Vudragovic et al., Comput. Phys. Commun. 183, 2021 (2012)) for calculating both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation in three spatial dimensions. The GP equation describes the properties of dilute Bose-Einstein condensates at ultra-cold temperatures. Hybrid versions of programs use the same algorithms as the C ones, involving real- and imaginary-time propagation based on a split-step Crank-Nicolson method, but consider only a fully-anisotropic three-dimensional GP equation, where algorithmic complexity for large grid sizes necessitates parallelization in order to reduce execution time and/or memory requirements per node. Since distributed memory approach is required to address the latter, we combine MPI programing paradigm with existing OpenMP codes, thus creating fully flexible parallelism within a combined distributed/shared memory model, suitable for different modern computer architectures. The two presented C/OpenMP/MPI programs for real- and imaginary-time propagation are optimized and accompanied by a customizable makefile. We present typical scalability results for the provided OpenMP/MPI codes and demonstrate almost linear speedup until inter-process communication time starts to dominate over calculation time per iteration. Such a scalability study is necessary for large grid sizes in order to determine optimal number of MPI nodes and OpenMP threads per node.

preprint2013arXiv

Scaling exponents and phase separation in a nonlinear network model inspired by the gravitational accretion

We study dynamics and scaling exponents in a nonlinear network model inspired by the formation of planetary systems. Dynamics of this model leads to phase separation to two types of condensate, light and heavy, distinguished by how they scale with mass. Light condensate distributions obey power laws given in terms of several identified scaling exponents that do not depend on initial conditions. The analyzed properties of heavy condensates have been found to be scale-free. Calculated mass distributions agree well with more complex models, and fit observations of both our own Solar System, and the best observed extra-solar planetary systems.

preprint2011arXiv

Fast Converging Path Integrals for Time-Dependent Potentials I: Recursive Calculation of Short-Time Expansion of the Propagator

In this and subsequent paper arXiv:1011.5185 we develop a recursive approach for calculating the short-time expansion of the propagator for a general quantum system in a time-dependent potential to orders that have not yet been accessible before. To this end the propagator is expressed in terms of a discretized effective potential, for which we derive and analytically solve a set of efficient recursion relations. Such a discretized effective potential can be used to substantially speed up numerical Monte Carlo simulations for path integrals, or to set up various analytic approximation techniques to study properties of quantum systems in time-dependent potentials. The analytically derived results are numerically verified by treating several simple models.

preprint2011arXiv

Fast Converging Path Integrals for Time-Dependent Potentials II: Generalization to Many-body Systems and Real-Time Formalism

Based on a previously developed recursive approach for calculating the short-time expansion of the propagator for systems with time-independent potentials and its time-dependent generalization for simple single-particle systems, in this paper we present a full extension of this formalism to a general quantum system with many degrees of freedom in a time-dependent potential. Furthermore, we also present a recursive approach for the velocity-independent part of the effective potential, which is necessary for calculating diagonal amplitudes and partition functions, as well as an extension from the imaginary-time formalism to the real-time one, which enables to study the dynamical properties of quantum systems. The recursive approach developed here allows an analytic derivation of the short-time expansion to orders that have not been accessible before, using the implemented SPEEDUP symbolic calculation code. The analytically derived results are extensively numerically verified by treating several models in both imaginary and real time.

preprint2011arXiv

SPEEDUP Code for Calculation of Transition Amplitudes via the Effective Action Approach

We present Path Integral Monte Carlo C code for calculation of quantum mechanical transition amplitudes for 1D models. The SPEEDUP C code is based on the use of higher-order short-time effective actions and implemented to the maximal order $p$=18 in the time of propagation (Monte Carlo time step), which substantially improves the convergence of discretized amplitudes to their exact continuum values. Symbolic derivation of higher-order effective actions is implemented in SPEEDUP Mathematica codes, using the recursive Schroedinger equation approach. In addition to the general 1D quantum theory, developed Mathematica codes are capable of calculating effective actions for specific models, for general 2D and 3D potentials, as well as for a general many-body theory in arbitrary number of spatial dimensions.

preprint2009arXiv

Properties of Quantum Systems via Diagonalization of Transition Amplitudes I: Discretization Effects

We analyze the method for calculation of properties of non-relativistic quantum systems based on exact diagonalization of space-discretized short-time evolution operators. In this paper we present a detailed analysis of the errors associated with space discretization. Approaches using direct diagonalization of real-space discretized Hamiltonians lead to polynomial errors in discretization spacing $Δ$. Here we show that the method based on the diagonalization of the short-time evolution operators leads to substantially smaller discretization errors, vanishing exponentially with $1/Δ^2$. As a result, the presented calculation scheme is particularly well suited for numerical studies of few-body quantum systems. The analytically derived discretization errors estimates are numerically shown to hold for several models. In the followup paper [1] we present and analyze substantial improvements that result from the merger of this approach with the recently introduced effective-action scheme for high-precision calculation of short-time propagation.

preprint2009arXiv

Properties of Quantum Systems via Diagonalization of Transition Amplitudes II: Systematic Improvements of Short-time Propagation

In this paper, building on a previous analysis [1] of exact diagonalization of the space-discretized evolution operator for the study of properties of non-relativistic quantum systems, we present a substantial improvement to this method. We apply recently introduced effective action approach for obtaining short-time expansion of the propagator up to very high orders to calculate matrix elements of space-discretized evolution operator. This improves by many orders of magnitude previously used approximations for discretized matrix elements and allows us to numerically obtain large numbers of accurate energy eigenvalues and eigenstates using numerical diagonalization. We illustrate this approach on several one and two-dimensional models. The quality of numerically calculated higher order eigenstates is assessed by comparison with semiclassical cumulative density of states.

preprint2008arXiv

Fast Convergence of Path Integrals for Many-body Systems

We generalize a recently developed method for accelerated Monte Carlo calculation of path integrals to the physically relevant case of generic many-body systems. This is done by developing an analytic procedure for constructing a hierarchy of effective actions leading to improvements in convergence of $N$-fold discretized many-body path integral expressions from 1/N to $1/N^p$ for generic $p$. In this paper we present explicit solutions within this hierarchy up to level $p=5$. Using this we calculate the low lying energy levels of a two particle model with quartic interactions for several values of coupling and demonstrate agreement with analytical results governing the increase in efficiency of the new method. The applicability of the developed scheme is further extended to the calculation of energy expectation values through the construction of associated energy estimators exhibiting the same speedup in convergence.

preprint2005arXiv

Asymptotic Properties of Path Integral Ideals

We introduce and analyze a new quantity, the path integral ideal, governing the flow of generic discrete theories to the continuum limit and greatly increasing their convergence. The said flow is classified according to the degree of divergence of the potential at spatial infinity. Studying the asymptotic behavior of path integral ideals we isolate the dominant terms in the effective potential that determine the behavior of a generic theory for large discrete time steps.

preprint2005arXiv

Generalization of Euler's Summation Formula to Path Integrals

A recently developed analytical method for systematic improvement of the convergence of path integrals is used to derive a generalization of Euler's summation formula for path integrals. The first $p$ terms in this formula improve convergence of path integrals to the continuum limit from 1/N to $1/N^p$, where $N$ is the coarseness of the discretization. Monte Carlo simulations performed on several different models show that the analytically derived speedup holds.

preprint2005arXiv

Jaggedness of Path Integral Trajectories

We define and investigate the properties of the jaggedness of path integral trajectories. The new quantity is shown to be scale invariant and to satisfy a self-averaging property. Jaggedness allows for a classification of path integral trajectories according to their relevance. We show that in the continuum limit the only paths that are not of measure zero are those with jaggedness 1/2, i.e. belonging to the same equivalence class as random walks. The set of relevant trajectories is thus narrowed down to a specific subset of non-differentiable paths. For numerical calculations, we show that jaggedness represents an important practical criterion for assessing the quality of trajectory generating algorithms. We illustrate the obtained results with Monte Carlo simulations of several different models.

preprint2005arXiv

Systematic Speedup of Path Integrals of a Generic $N$-fold Discretized Theory

We present and discuss a detailed derivation of a new analytical method that systematically improves the convergence of path integrals of a generic $N$-fold discretized theory. We develop an explicit procedure for calculating a set of effective actions $S^{(p)}$, for $p=1,2,3,...$ which have the property that they lead to the same continuum amplitudes as the starting action, but that converge to that continuum limit ever faster. Discretized amplitudes calculated using the $p$ level effective action differ from the continuum limit by a term of order $1/N^p$. We obtain explicit expressions for the effective actions for levels $p\le 9$. We end by analyzing the speedup of Monte Carlo simulations of two different models: an anharmonic oscillator with quartic coupling and a particle in a modified Pöschl-Teller potential.

preprint2005arXiv

Systematically Accelerated Convergence of Path Integrals

We present a new analytical method that systematically improves the convergence of path integrals of a generic $N$-fold discretized theory. Using it we calculate the effective actions $S^{(p)}$ for $p\le 9$ which lead to the same continuum amplitudes as the starting action, but that converge to that continuum limit as $1/N^p$. We checked this derived speedup in convergence by performing Monte Carlo simulations on several different models.

preprint1998arXiv

An Improved Gaussian Approximation for Quantum Field Theory

We present a new approximation technique for quantum field theory. The standard one-loop result is used as a seed for a recursive formula that gives a sequence of improved Gaussian approximations for the generating functional. In a different setting, the basic idea of this recursive scheme is used in the second part of the paper to substantialy speed up the standard Monte Carlo algorithm.

preprint1998arXiv

Improved Gaussian Approximation

In a recently developed approximation technique for quantum field theory the standard one-loop result is used as a seed for a recursive formula that gives a sequence of improved Gaussian approximations for the generating functional. In this paper we work with the generic $ϕ^3+ϕ^4$ model in $d=0$ dimensions. We compare the first, and simplest, approximation in the above sequence with the one-loop and two-loop approximations, as well as the exact results (calculated numericaly).