BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Balsa Terzic

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

ResearcherUnclaimed source record
astro-phDistributed, Parallel, and Cluster ComputingPerformancephysics.acc-ph

Catalog footprint

What is connected

6works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

6 published item(s)

preprint2022arXiv

m-CUBES An efficient and portable implementation of multi-dimensional integration for gpus

The task of multi-dimensional numerical integration is frequently encountered in physics and other scientific fields, e.g., in modeling the effects of systematic uncertainties in physical systems and in Bayesian parameter estimation. Multi-dimensional integration is often time-prohibitive on CPUs. Efficient implementation on many-core architectures is challenging as the workload across the integration space cannot be predicted a priori. We propose m-Cubes, a novel implementation of the well-known Vegas algorithm for execution on GPUs. Vegas transforms integration variables followed by calculation of a Monte Carlo integral estimate using adaptive partitioning of the resulting space. m-Cubes improves performance on GPUs by maintaining relatively uniform workload across the processors. As a result, our optimized Cuda implementation for Nvidia GPUs outperforms parallelization approaches proposed in past literature. We further demonstrate the efficiency of m-Cubes by evaluating a six-dimensional integral from a cosmology application, achieving significant speedup and greater precision than the CUBA library's CPU implementation of VEGAS. We also evaluate m-Cubes on a standard integrand test suite. m-Cubes outperforms the serial implementations of the Cuba and GSL libraries by orders of magnitude speedup while maintaining comparable accuracy. Our approach yields a speedup of at least 10 when compared against publicly available Monte Carlo based GPU implementations. In summary, m-Cubes can solve integrals that are prohibitively expensive using standard libraries and custom implementations. A modern C++ interface header-only implementation makes m-Cubes portable, allowing its utilization in complicated pipelines with easy to define stateful integrals. Compatibility with non-Nvidia GPUs is achieved with our initial implementation of m-Cubes using the Kokkos framework.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Laser Pulsing in Linear Compton Scattering

Previous work on calculating energy spectra from Compton scattering events has either neglected considering the pulsed structure of the incident laser beam, or has calculated these effects in an approximate way subject to criticism. In this paper, this problem has been reconsidered within a linear plane wave model for the incident laser beam. By performing the proper Lorentz transformation of the Klein-Nishina scattering cross section, a spectrum calculation can be created which allows the electron beam energy spread and emittance effects on the spectrum to be accurately calculated, essentially by summing over the emission of each individual electron. Such an approach has the obvious advantage that it is easily integrated with a particle distribution generated by particle tracking, allowing precise calculations of spectra for realistic particle distributions in collision. The method is used to predict the energy spectrum of radiation passing through an aperture for the proposed Old Dominion University inverse Compton source. Many of the results allow easy scaling estimates to be made of the expected spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Empirical models for Dark Matter Halos. I. Nonparametric Construction of Density Profiles and Comparison with Parametric Models

We use techniques from nonparametric function estimation theory to extract the density profiles, and their derivatives, from a set of N-body dark matter halos. We consider halos generated from LCDM simulations of gravitational clustering, as well as isolated, spherical collapses. The logarithmic density slopes gamma = d(log rho)/d(log r) of the LCDM halos are found to vary as power-laws in radius, reaching values of gamma ~ -1 at the innermost resolved radii (~0.01 r_virial). This behavior is significantly different from that of broken power-law models like the NFW profile, but similar to that of models like de Vaucouleurs'. Accordingly, we compare the N-body density profiles with various parametric models to find which provide the best fit. We consider an NFW-like model with arbitrary inner slope; Dehnen & McLaughlin's anisotropic model; Einasto's model (identical in functional form to Sersic's model but fit to the space density); and the density model of Prugniel & Simien that was designed to match the deprojected form of Sersic's R^{1/n} law. Overall, the best-fitting model to the LCDM halos is Einasto's, although the Prugniel-Simien and Dehnen-McLaughlin models also perform well. With regard to the spherical collapse halos, both the Prugniel-Simien and Einasto models describe the density profiles well, with an rms scatter some four times smaller than that obtained with either the NFW-like model or the 3-parameter Dehnen-McLaughlin model. Finally, we confirm recent claims of a systematic variation in profile shape with halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Empirical Models for Dark Matter Halos. II. Inner profile slopes, dynamical profiles, and rho/sigma^3

We have recently shown that both the Prugniel-Simien model and Sersic's function (hereafter referred to as the Einasto model when applied to internal density profiles) describe simulated dark matter halos better than an NFW-like model with equal number of parameters. Here we provide analytical expressions for the logarithmic slopes of these models, and compare them with data from real galaxies. Depending on the Einasto parameters of the dark matter halo, one can expect an extrapolated, inner (0.01--1 kpc), logarithmic profile slope ranging from -0.2 to -1.5, with a typical value at 0.1 kpc around -0.7. Application of this (better fitting) model therefore alleviates some of the past disagreement with observations on this issue. We additionally provide useful expressions for the concentration and assorted scale radii: r_s, r_{-2}, r_e, R_e, r_virial, and r_max -- the radius where the circular velocity profile has its maximum value. We also present the circular velocity profiles and the radial behavior of rho(r)/sigma(r)^3 for both the Einasto and Prugniel-Simien models. We find this representation of the phase-space density profile to be well approximated by a power-law with slope slightly shallower than -2 near r=r_{-2}.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Empirical Models for Dark Matter Halos. III. The Kormendy relation and the log(rho_e)-log(R_e) relation

We have recently shown that the 3-parameter density-profile model from Prugniel & Simien provides a better fit to simulated, galaxy- and cluster-sized, dark matter halos than an NFW-like model with arbitrary inner profile slope gamma (Paper I). By construction, the parameters of the Prugniel-Simien model equate to those of the Sersic R^{1/n} function fitted to the projected distribution. Using the Prugniel-Simien model, we are therefore able to show that the location of simulated (10^{12} M_sun) galaxy-sized dark matter halos in the <mu>_e-log(R_e) diagram coincides with that of brightest cluster galaxies, i.e., the dark matter halos appear consistent with the Kormendy relation defined by luminous elliptical galaxies. These objects are also seen to define the new, and equally strong, relation log(rho_e) = 0.5 - 2.5log(R_e), in which rho_e is the internal density at r=R_e. Simulated (10^{14.5} M_sun) cluster-sized dark matter halos and the gas component of real galaxy clusters follow the relation log(rho_e) = 2.5[1 - log(R_e)]. Given the shapes of the various density profiles, we are able to conclude that while dwarf elliptical galaxies and galaxy clusters can have dark matter halos with effective radii of comparable size to the effective radii of their baryonic component, luminous elliptical galaxies can not. For increasingly large elliptical galaxies, with increasingly large profile shapes `n', to be dark matter dominated at large radii requires dark matter halos with increasingly large effective radii compared to the effective radii of their stellar component.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Density-potential pairs for spherical stellar systems with Sersic light-profiles and (optional) power-law cores

Popular models for describing the luminosity-density profiles of dynamically hot stellar systems (e.g., Jaffe, Hernquist, Dehnen) were constructed with the desire to match the deprojected form of an R^{1/4} light-profile. Real galaxies, however, are now known to have a range of different light-profile shapes that scale with mass. Consequently, although highly useful, the above models have implicit limitations, and this is illustrated here through their application to a number of real galaxy density profiles. On the other hand, the analytical density profile given by Prugniel & Simien (1997) closely matches the deprojected form of Sersic R^{1/n} light-profiles - including deprojected exponential light-profiles. It is thus applicable for describing bulges in spiral galaxies, dwarf elliptical galaxies, and both ordinary and giant elliptical galaxies. Here we provide simple equations, in terms of elementary and special functions, for the gravitational potential and force associated with this density profile. Furthermore, to match galaxies with partially depleted cores, and better explore the supermassive black hole / galaxy connection, we have added a power-law core to this density profile and derived similar expressions for the potential and force of this hybrid profile. Expressions for the mass and velocity dispersion, assuming isotropy, are also given. These models may also prove appropriate for describing the dark matter distribution in halos formed from LCDM cosmological simulations.

0 citations0 reviewsNo code linkedOpen access