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

Nicolas Bock

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

ResearcherUnclaimed source record
Numerical Analysiscond-mat.mtrl-sciData Structures and Algorithmscond-mat.stat-mechhep-exMathematical Softwarenucl-exphysics.chem-phphysics.comp-phcond-mat.othermath.NAPerformance

Catalog footprint

What is connected

14works
12topics
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

14 published item(s)

preprint2021arXiv

Performance Optimizations of Recursive Electronic Structure Solvers targeting Multi-Core Architectures (LA-UR-20-26665)

As we rapidly approach the frontiers of ultra large computing resources, software optimization is becoming of paramount interest to scientific application developers interested in efficiently leveraging all available on-Node computing capabilities and thereby improving a requisite science per watt metric. The scientific application of interest here is the Basic Math Library (BML) that provides a singular interface for linear algebra operation frequently used in the Quantum Molecular Dynamics (QMD) community. The provisioning of a singular interface indicates the presence of an abstraction layer which in-turn suggests commonalities in the code-base and therefore any optimization or tuning introduced in the core of code-base has the ability to positively affect the performance of the aforementioned library as a whole. With that in mind, we proceed with this investigation by performing a survey of the entirety of the BML code-base, and extract, in form of micro-kernels, common snippets of code. We introduce several optimization strategies into these micro-kernels including 1.) Strength Reduction 2.) Memory Alignment for large arrays 3.) Non Uniform Memory Access (NUMA) aware allocations to enforce data locality and 4.) appropriate thread affinity and bindings to enhance the overall multi-threaded performance. After introducing these optimizations, we benchmark the micro-kernels and compare the run-time before and after optimization for several target architectures. Finally we use the results as a guide to propagating the optimization strategies into the BML code-base. As a demonstration, herein, we test the efficacy of these optimization strategies by comparing the benchmark and optimized versions of the code.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Graph-based linear scaling electronic structure theory

We show how graph theory can be combined with quantum theory to calculate the electronic structure of large complex systems. The graph formalism is general and applicable to a broad range of electronic structure methods and materials, including challenging systems such as biomolecules. The methodology combines well-controlled accuracy, low computational cost, and natural low-communication parallelism. This combination addresses substantial shortcomings of linear scaling electronic structure theory, in particular with respect to quantum-based molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A $N$-Body Solver for Square Root Iteration

We develop the Sparse Approximate Matrix Multiply ($\tt SpAMM$) $n$-body solver for first order Newton Schulz iteration of the matrix square root and inverse square root. The solver performs recursive two-sided metric queries on a modified Cauchy-Schwarz criterion, culling negligible sub-volumes of the product-tensor for problems with structured decay in the sub-space metric. These sub-structures are shown to bound the relative error in the matrix-matrix product, and in favorable cases, to enjoy a reduced computational complexity governed by dimensionality reduction of the product volume. A main contribution is demonstration of a new, algebraic locality that develops under contractive identity iteration, with collapse of the metric-subspace onto the identity's plane diagonal, resulting in a stronger $\tt SpAMM$ bound. Also, we carry out a first order {Fréchet} analyses for single and dual channel instances of the square root iteration, and look at bifurcations due to ill-conditioning and a too aggressive $\tt SpAMM$ approximation. Then, we show that extreme $\tt SpAMM$ approximation and contractive identity iteration can be achieved for ill-conditioned systems through regularization, and we demonstrate the potential for acceleration with a scoping, product representation of the inverse factor.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solvers for $\mathcal{O} (N)$ Electronic Structure in the Strong Scaling Limit

We present a hybrid OpenMP/Charm++ framework for solving the $\mathcal{O} (N)$ Self-Consistent-Field eigenvalue problem with parallelism in the strong scaling regime, $P\gg{N}$, where $P$ is the number of cores, and $N$ a measure of system size, i.e. the number of matrix rows/columns, basis functions, atoms, molecules, etc. This result is achieved with a nested approach to Spectral Projection and the Sparse Approximate Matrix Multiply [Bock and Challacombe, SIAM J.~Sci.~Comput. 35 C72, 2013], and involves a recursive, task-parallel algorithm, often employed by generalized $N$-Body solvers, to occlusion and culling of negligible products in the case of matrices with decay. Employing classic technologies associated with generalized $N$-Body solvers, including over-decomposition, recursive task parallelism, orderings that preserve locality, and persistence-based load balancing, we obtain scaling beyond hundreds of cores per molecule for small water clusters ([H${}_2$O]${}_N$, $N \in \{ 30, 90, 150 \}$, $P/N \approx \{ 819, 273, 164 \}$) and find support for an increasingly strong scalability with increasing system size $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An N-Body Solution to the Problem of Fock Exchange

We report an N-Body approach to computing the Fock exchange matrix with and without permutational symmetry. The method achieves an O(N lg N) computational complexity through an embedded metric-query, allowing hierarchical application of direct SCF criteria. The advantages of permutational symmetry are found to be 4-fold for small systems, but decreasing with increasing system size and/or more permissive neglect criteria. This work sets the stage for: (1) the introduction of range queries in multi-level multipole schemes for rank reduction, and (2) recursive task parallelism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal electronic excitations in liquid metals

Thermal electronic excitations in metal crystals are calculated by starting with a reference structure for the nuclei: the crystal structure of the appropriate phase. Here we explain the corresponding theory for metal liquids, starting with an appropriate reference structure for a liquid. We explain the significance of these structures, and we briefly review how to find them and calculate their properties. Then we examine the electronic densities of states for liquid structures of Na, Al, and Cu, comparing them to their crystal forms. Next we explain how to calculate the dominant electronic thermal excitation term, considering issues of accuracy that do not arise in the crystal theory. Finally we briefly discuss the contribution from the interaction between excited electrons and moving nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Optimized Sparse Approximate Matrix Multiply for Matrices with Decay

We present an optimized single-precision implementation of the Sparse Approximate Matrix Multiply (\SpAMM{}) [M. Challacombe and N. Bock, arXiv {\bf 1011.3534} (2010)], a fast algorithm for matrix-matrix multiplication for matrices with decay that achieves an $\mathcal{O} (n \log n)$ computational complexity with respect to matrix dimension $n$. We find that the max norm of the error achieved with a \SpAMM{} tolerance below $2 \times 10^{-8}$ is lower than that of the single-precision {\tt SGEMM} for dense quantum chemical matrices, while outperforming {\tt SGEMM} with a cross-over already for small matrices ($n \sim 1000$). Relative to naive implementations of \SpAMM{} using Intel's Math Kernel Library ({\tt MKL}) or AMD's Core Math Library ({\tt ACML}), our optimized version is found to be significantly faster. Detailed performance comparisons are made for quantum chemical matrices with differently structured sub-blocks. Finally, we discuss the potential of improved hardware prefetch to yield 2--3x speedups.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Energy momentum conservation effects on two-particle correlation functions

Two particle correlations are used to extract information about the characteristic size of the system in proton-proton and heavy ion collisions. The size of the system can be extracted from the Bose-Einstein quantum mechanical effect for identical particles. However there are also long range correlations that shift the baseline of the correlation function from the expected flat behavior. A possible source of these correlations is the conservation of energy and momentum, especially for small systems, where the energy available for particle production is limited. A new technique, first used by the STAR collaboration, of quantifying these long range correlations using energy-momentum conservation considerations is presented in this talk. Using Monte Carlo simulations of proton-proton collisions at 900 GeV, it is shown that the baseline of the two particle correlation function can be described using this technique.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Extended Lagrangian free energy molecular dynamics

Extended free energy Lagrangians are proposed for first principles molecular dynamics simulations at finite electronic temperatures for plane-wave pseudopotential and local orbital density matrix based calculations. Thanks to the extended Lagrangian description the electronic degrees of freedom can be integrated by stable geometric schemes that conserve the free energy. For the local orbital representations both the nuclear and electronic forces have simple and numerically efficient expressions that are well suited for reduced complexity calculations. A rapidly converging recursive Fermi operator expansion method that does not require the calculation of eigenvalues and eigenfunctions for the construction of the fractionally occupied density matrix is discussed. An efficient expression for the Pulay force that is valid also for density matrices with fractional occupation occurring at finite electronic temperatures is also demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast Multiplication of Matrices with Decay

A fast algorithm for the approximate multiplication of matrices with decay is introduced; the Sparse Approximate Matrix Multiply (SpAMM) reduces complexity in the product space, a different approach from current methods that economize within the matrix space through truncation or rank reduction. Matrix truncation (element dropping) is compared to SpAMM for quantum chemical matrices with approximate exponential and algebraic decay. For matched errors in the electronic total energy, SpAMM is found to require fewer to far fewer floating point operations relative to dropping. The challenges and opportunities afforded by this new approach are discussed, including the potential for high performance implementations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Femtoscopy and energy-momentum conservation effects in proton-proton collisions at 900 GeV in ALICE

Two particle correlations are used to extract information about the characteristic size of the system for proton-proton collisions at 900 GeV measured by the ALICE (A Large Ion Collider experiment) detector at CERN. The correlation functions obtained show the expected Bose-Einstein effect for identical particles, but there are also long range correlations present that shift the baseline from the expected flat behavior. A possible source of these correlations is the conservation of energy and momentum, especially for small systems, where the energy available for particle production is limited. A new technique, first introduced by the STAR collaboration, of quantifying these long range correlations using energy-momentum conservation considerations is presented here. It is shown that the baseline of the two particle correlation function can be described using this technique.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Liquid state properties from first principles DFT calculations: Static properties

In order to test the Vibration-Transit (V-T) theory of liquid dynamics, ab initio density functional theory (DFT) calculations of thermodynamic properties of Na and Cu are performed and compared with experimental data. The calculations are done for the crystal at T = 0 and T_m, and for the liquid at T_m. The key theoretical quantities for crystal and liquid are the structural potential and the dynamical matrix, both as function of volume. The theoretical equations are presented, as well as details of the DFT computations. The properties compared with experiment are the equilibrium volume, the isothermal bulk modulus, the internal energy and the entropy. The agreement of theory with experiment is uniformly good. Our primary conclusion is that the application of DFT to V-T theory is feasible, and the resulting liquid calculations achieve the same level of accuracy as does ab initio lattice dynamics for crystals. Moreover, given the well established reliability of DFT, the present results provide a significant confirmation of V-T theory itself.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An improved model for the transit entropy of monatomic liquids

In the original formulation of vibration-transit (V-T) theory for monatomic liquid dynamics, the transit contribution to entropy was taken to be a universal constant, calibrated to the constant-volume entropy of melting. This model suffers two deficiencies: (a) it does not account for experimental entropy differences of 2% among elemental liquids, and (b) it implies a value of zero for the transit contribution to internal energy. The purpose of this paper is to correct these deficiencies. To this end, the V-T equation for entropy is fitted to an overall accuracy of 0.1% to the available experimental high temperature entropy data for elemental liquids. The theory contains two nuclear motion contributions: (a) the dominant vibrational contribution S_{vib}(T/θ_0), where T is temperature and θ_0 is the vibrational characteristic temperature, and (b) the transit contribution S_{tr}(T/θ_{tr}), where θ_{tr} is a scaling temperature for each liquid. The appearance of a common functional form of S_{tr} for all the liquids studied is a property of the experimental data, when analyzed via the V-T formula. The resulting S_{tr} implies the correct transit contribution to internal energy. The theoretical entropy of melting is derived, in a single formula applying to normal and anomalous melting alike. An ab initio calculation of θ_0, based on density functional theory, is reported for liquid Na and Cu. Comparison of these calculations with the above analysis of experimental entropy data provides verification of V-T theory. In view of the present results, techniques currently being applied in ab initio simulations of liquid properties can be employed to advantage in the further testing and development of V-T theory.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Density Analysis of Network Community Divisions

We present a compact matrix formulation of the modularity, a commonly used quality measure for the community division in a network. Using this formulation we calculate the density of modularities, a statistical measure of the probability of finding a particular modularity for a random but valid community division into $C$ communities. We present our results for some well--known and some artificial networks, and we conclude that the general features of the modularity density are quite similar for the different networks. From a simple model of the modularity we conclude that all nnected networks must show similar shapes of their modularity densities. The general features of this density may give valuable information in the search for good optimization schemes of the modularity.

0 citations0 reviewsNo code linkedOpen access