Catalog footprint

What is connected

21works
21topics
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

21 published item(s)

preprint2022arXiv

A Scalable Approach to Quantum Simulation via Projection-based Embedding

Owing to the computational complexity of electronic structure algorithms running on classical digital computers, the range of molecular systems amenable to simulation remains tightly circumscribed even after many decades of work. Quantum computers hold the promise of transcending such limitations although in the current era the size and noise of these devices militates against significant progress. Here we describe a new and chemically intuitive approach that permits a subdomain of the electronic structure of a molecule to be calculated accurately on a quantum device, while the rest of the molecule is described at a lower level of accuracy using density functional theory running on a classical computer. We demonstrate that our method produces improved results for molecules that cannot be simulated fully on quantum computers but which can be resolved classically at a lower level of approximation. Our algorithm is tunable, so that the size of the quantum simulation can be adjusted to run on available quantum resources. Therefore, as quantum devices become larger, our method will enable increasingly large subdomains to be studied accurately.

preprint2019arXiv

Ordering of Trotterization: Impact on Errors in Quantum Simulation of Electronic Structure

Trotter-Suzuki decompositions are frequently used in the quantum simulation of quantum chemistry. They transform the evolution operator into a form implementable on a quantum device, while incurring an error---the Trotter error. The Trotter error can be made arbitrarily small by increasing the Trotter number. However, this increases the length of the quantum circuits required, which may be impractical. It is therefore desirable to find methods of reducing the Trotter error through alternate means. The Trotter error is dependent on the order in which individual term unitaries are applied. Due to the factorial growth in the number of possible orderings with respect to the number of terms, finding an optimal strategy for ordering Trotter sequences is difficult. In this paper, we propose three ordering strategies, and assess their impact on the Trotter error incurred. Initially, we exhaustively examine the possible orderings for molecular hydrogen in a STO-3G basis. We demonstrate how the optimal ordering scheme depends on the compatibility graph of the Hamiltonian, and show how it varies with increasing bond length. We then use 44 molecular Hamiltonians to evaluate two strategies based on coloring their incompatibility graphs, while considering the properties of the obtained colorings. We find that the Trotter error for most systems involving heavy atoms, using a reference magnitude ordering, is less than 1 kcal/mol. Relative to this, the difference between ordering schemes can be substantial, being approximately on the order of millihartrees. The coloring-based ordering schemes are reasonably promising, however further work is required. Finally, we consider ordering strategies based on the norm of the Trotter error operator, including an iterative method for generating the new error operator terms added upon insertion of a term into an ordered Hamiltonian.

preprint2015arXiv

A Distributed Multi-agent Market Place for HPC Compute Cycle Resource Trading

Computer simulation is finding a role in an increasing number of scientific disciplines, concomitant with the rise in available computing power. Realizing this inevitably requires access to computational power beyond the desktop, making use of clusters, supercomputers, data repositories, networks and distributed aggregations of these resources. Accessing one such resource entails a number of usability and security problems; when multiple geographically distributed resources are involved, the difficulty is compounded. This presents the user with the problem of how to gain access to suitable resources to run their workloads as they need them. In this paper we present our solutions to this problem, a resource trading platform that allows users to purchase access to resources within a distributed e-infrastructure. We present the implementation of this Resource Allocation Market Place as a distributed multi-agent system, and show how it provides a highly flexible, efficient tool to schedule workflows across high performance computing resources.

preprint2015arXiv

An automated multiscale ensemble simulation approach for vascular blood flow

Cerebrovascular diseases such as brain aneurysms are a primary cause of adult disability. The flow dynamics in brain arteries, both during periods of rest and increased activity, are known to be a major factor in the risk of aneurysm formation and rupture. The precise relation is however still an open field of investigation. We present an automated ensemble simulation method for modelling cerebrovascular blood flow under a range of flow regimes. By automatically constructing and performing an ensemble of multiscale simulations, where we unidirectionally couple a 1D solver with a 3D lattice-Boltzmann code, we are able to model the blood flow in a patient artery over a range of flow regimes. We apply the method to a model of a middle cerebral artery, and find that this approach helps us to fine-tune our modelling techniques, and opens up new ways to investigate cerebrovascular flow properties.

preprint2014arXiv

Assembling Ellipsoidal Particles at Fluid Interfaces using Switchable Dipolar Capillary Interactions

The fabrication of novel soft materials is an important scientific and technological challenge. We investigate the response of magnetic ellipsoidal particles adsorbed at fluid-fluid interfaces to external magnetic fields. By exploiting previously discovered first-order orientation phase transitions we show how to switch on and off dipolar capillary interactions between particles, leading to the formation of distinctive self-assembled structures and allowing dynamic control of the bottom-up fabrication of reconfigurable novel-structured materials

preprint2014arXiv

Computer simulations reveal complex distribution of haemodynamic forces in a mouse retina model of angiogenesis

There is currently limited understanding of the role played by haemodynamic forces on the processes governing vascular development. One of many obstacles to be overcome is being able to measure those forces, at the required resolution level, on vessels only a few micrometres thick. In the current paper, we present an in silico method for the computation of the haemodynamic forces experienced by murine retinal vasculature (a widely used vascular development animal model) beyond what is measurable experimentally. Our results show that it is possible to reconstruct high-resolution three-dimensional geometrical models directly from samples of retinal vasculature and that the lattice-Boltzmann algorithm can be used to obtain accurate estimates of the haemodynamics in these domains. We generate flow models from samples obtained at postnatal days (P) 5 and 6. Our simulations show important differences between the flow patterns recovered in both cases, including observations of regression occurring in areas where wall shear stress gradients exist. We propose two possible mechanisms to account for the observed increase in velocity and wall shear stress between P5 and P6: i) the measured reduction in typical vessel diameter between both time points, ii) the reduction in network density triggered by the pruning process. The methodology developed herein is applicable to other biomedical domains where microvasculature can be imaged but experimental flow measurements are unavailable or difficult to obtain.

preprint2014arXiv

Deformability-based red blood cell separation in deterministic lateral displacement devices - a simulation study

We show, via three-dimensional immersed-boundary-finite-element-lattice-Boltzmann simulations, that deformability-based red blood cell (RBC) separation in deterministic lateral displacement (DLD) devices is possible. This is due to the deformability-dependent lateral extension of RBCs and enables us to predict a priori which RBCs will be displaced in a given DLD geometry. Several diseases affect the deformability of human cells. Malaria-infected RBCs or sickle cells, for example, tend to become stiffer than their healthy counterparts. It is therefore desirable to design microfluidic devices which can detect those diseases based on the cells' deformability fingerprint, rather than preparing samples using expensive and time-consuming biochemical preparation steps. Our findings should be helpful in the development of new methods for sorting cells and particles by deformability.

preprint2014arXiv

Detachment Energies of Spheroidal Particles from Fluid-Fluid Interfaces

The energy required to detach a single particle from a fluid-fluid interface is an important parameter for designing certain soft materials, for example, emulsions stabilised by colloidal particles, colloidosomes designed for targeted drug delivery, and bio-sensors composed of magnetic particles adsorbed at interfaces. For a fixed particle volume, prolate and oblate spheroids attach more strongly to interfaces because they have larger particle-interface areas. Calculating the detachment energy of spheroids necessitates the difficult measurement of particle-liquid surface tensions, in contrast with spheres, where the contact angle suffices. We develop a simplified detachment energy model for spheroids which depends only on the particle aspect ratio and the height of the particle centre of mass above the fluid-fluid interface. We use lattice Boltzmann simulations to validate the model and provide quantitative evidence that the approach can be applied to simulate particle-stabilized emulsions, and highlight the experimental implications of this validation.

preprint2014arXiv

Impact of blood rheology on wall shear stress in a model of the middle cerebral artery

Perturbations to the homeostatic distribution of mechanical forces exerted by blood on the endothelial layer have been correlated with vascular pathologies including intracranial aneurysms and atherosclerosis. Recent computational work suggests that in order to correctly characterise such forces, the shear-thinning properties of blood must be taken into account. To the best of our knowledge, these findings have never been compared against experimentally observed pathological thresholds. In the current work, we apply the three-band diagram (TBD) analysis due to Gizzi et al. to assess the impact of the choice of blood rheology model on a computational model of the right middle cerebral artery. Our results show that, in the model under study, the differences between the wall shear stress predicted by a Newtonian model and the well known Carreau-Yasuda generalized Newtonian model are only significant if the vascular pathology under study is associated with a pathological threshold in the range 0.94 Pa to 1.56 Pa, where the results of the TBD analysis of the rheology models considered differs. Otherwise, we observe no significant differences.

preprint2014arXiv

Interface deformations affect the orientation transition of magnetic ellipsoidal particles adsorbed at fluid-fluid interfaces

Manufacturing new soft materials with specific optical, mechanical and magnetic properties is a significant challenge. Assembling and manipulating colloidal particles at fluid interfaces is a promising way to make such materials. We use lattice-Boltzmann simulations to investigate the response of magnetic ellipsoidal particles adsorbed at liquid-liquid interfaces to external magnetic fields. We provide further evidence for the first-order orientation phase transition predicted by Bresme and Faraudo [Journal of Physics: Condensed Matter 19 (2007), 375110]. We show that capillary interface deformations around the ellipsoidal particle significantly affect the tilt-angle of the particle for a given dipole-field strength, altering the properties of the orientation transition. We propose scaling laws governing this transition, and suggest how to use these deformations to facilitate particle assembly at fluid-fluid interfaces.

preprint2014arXiv

Weighted decomposition in high-performance lattice-Boltzmann simulations: are some lattice sites more equal than others?

Obtaining a good load balance is a significant challenge in scaling up lattice-Boltzmann simulations of realistic sparse problems to the exascale. Here we analyze the effect of weighted decomposition on the performance of the HemeLB lattice-Boltzmann simulation environment, when applied to sparse domains. Prior to domain decomposition, we assign wall and in/outlet sites with increased weights which reflect their increased computational cost. We combine our weighted decomposition with a second optimization, which is to sort the lattice sites according to a space filling curve. We tested these strategies on a sparse bifurcation and very sparse aneurysm geometry, and find that using weights reduces calculation load imbalance by up to 85%, although the overall communication overhead is higher than some of our runs.

preprint2013arXiv

Analyzing and Modeling the Performance of the HemeLB Lattice-Boltzmann Simulation Environment

We investigate the performance of the HemeLB lattice-Boltzmann simulator for cerebrovascular blood flow, aimed at providing timely and clinically relevant assistance to neurosurgeons. HemeLB is optimised for sparse geometries, supports interactive use, and scales well to 32,768 cores for problems with ~81 million lattice sites. We obtain a maximum performance of 29.5 billion site updates per second, with only an 11% slowdown for highly sparse problems (5% fluid fraction). We present steering and visualisation performance measurements and provide a model which allows users to predict the performance, thereby determining how to run simulations with maximum accuracy within time constraints.

preprint2013arXiv

Choice of boundary condition for lattice-Boltzmann simulation of moderate Reynolds number flow in complex domains

Modeling blood flow in larger vessels using lattice-Boltzmann methods comes with a challenging set of constraints: a complex geometry with walls and inlet/outlets at arbitrary orientations with respect to the lattice, intermediate Reynolds number, and unsteady flow. Simple bounce-back is one of the most commonly used, simplest, and most computationally efficient boundary conditions, but many others have been proposed. We implement three other methods applicable to complex geometries (Guo, Zheng and Shi, Phys Fluids (2002); Bouzdi, Firdaouss and Lallemand, Phys. Fluids (2001); Junk and Yang Phys. Rev. E (2005)) in our open-source application \HemeLB{}. We use these to simulate Poiseuille and Womersley flows in a cylindrical pipe with an arbitrary orientation at physiologically relevant Reynolds (1--300) and Womersley (4--12) numbers and steady flow in a curved pipe at relevant Dean number (100--200) and compare the accuracy to analytical solutions. We find that both the Bouzidi-Firdaouss-Lallemand and Guo-Zheng-Shi methods give second-order convergence in space while simple bounce-back degrades to first order. The BFL method appears to perform better than GZS in unsteady flows and is significantly less computationally expensive. The Junk-Yang method shows poor stability at larger Reynolds number and so cannot be recommended here. The choice of collision operator (lattice Bhatnagar-Gross-Krook vs.\ multiple relaxation time) and velocity set (D3Q15 vs.\ D3Q19 vs.\ D3Q27) does not significantly affect the accuracy in the problems studied.

preprint2013arXiv

Distributed Multiscale Computing with MUSCLE 2, the Multiscale Coupling Library and Environment

We present the Multiscale Coupling Library and Environment: MUSCLE 2. This multiscale component-based execution environment has a simple to use Java, C++, C, Python and Fortran API, compatible with MPI, OpenMP and threading codes. We demonstrate its local and distributed computing capabilities and compare its performance to MUSCLE 1, file copy, MPI, MPWide, and GridFTP. The local throughput of MPI is about two times higher, so very tightly coupled code should use MPI as a single submodel of MUSCLE 2; the distributed performance of GridFTP is lower, especially for small messages. We test the performance of a canal system model with MUSCLE 2, where it introduces an overhead as small as 5% compared to MPI.

preprint2013arXiv

Flexible composition and execution of high performance, high fidelity multiscale biomedical simulations

Multiscale simulations are essential in the biomedical domain to accurately model human physiology. We present a modular approach for designing, constructing and executing multiscale simulations on a wide range of resources, from desktops to petascale supercomputers, including combinations of these. Our work features two multiscale applications, in-stent restenosis and cerebrovascular bloodflow, which combine multiple existing single-scale applications to create a multiscale simulation. These applications can be efficiently coupled, deployed and executed on computers up to the largest (peta) scale, incurring a coupling overhead of 1 to 10% of the total execution time.

preprint2013arXiv

Rheology of Cubic Blue Phases

We study the behaviour of cubic blue phases under shear flow via lattice Boltzmann simulations. We focus on the two experimentally observed phases, Blue Phase I (BPI) and Blue Phase II (BPII). The disclination network of Blue Phase II continuously breaks and reforms under steady shear, leading to an oscillatory stress response in time. For larger shear rates, the structure breaks up into a Grandjean texture with a cholesteric helix lying along the flow gradient direction. Blue Phase I leads to a very different response. Here, oscillations are only possible for intermediate shear rates -- very slow flow causes a transition of the initially ordered structure into an amorphous network with an apparent yield stress. Larger shear rates lead to another amorphous state with different structure of the defect network. For even larger flow rates the same break-up into a Grandjean texture as for Blue Phase II is observed. At the highest imposed flow rates both cubic blue phases adopt a flow-aligned nematic state. Our results provide the first theoretical investigation of sheared blue phases in large systems, and are relevant to understanding the bulk rheology of these materials.

preprint2013arXiv

Survey of Multiscale and Multiphysics Applications and Communities

Multiscale and multiphysics applications are now commonplace, and many researchers focus on combining existing models to construct combined multiscale models. Here we present a concise review of multiscale applications and their source communities. We investigate the prevalence of multiscale projects in the EU and the US, review a range of coupling toolkits they use to construct multiscale models and identify areas where collaboration between disciplines could be particularly beneficial. We conclude that multiscale computing has become increasingly popular in recent years, that different communities adopt very different approaches to constructing multiscale simulations, and that simulations on a length scale of a few metres and a time scale of a few hours can be found in many of the multiscale research domains. Communities may receive additional benefit from sharing methods that are geared towards these scales.

preprint2012arXiv

Coalesced communication: a design pattern for complex parallel scientific software

We present a new design pattern for high-performance parallel scientific software, named coalesced communication. This pattern allows for a structured way to improve the communication performance through coalescence of multiple communication needs using two communication management components. We apply the design pattern to several simulations of a lattice-Boltzmann blood flow solver with streaming visualisation which engenders a reduction in the communication overhead of approximately 40%.

preprint2004arXiv

Energy controlled insertion of polar molecules in dense fluids

We present a method to search low energy configurations of polar molecules in the complex potential energy surfaces associated with dense fluids. The search is done in the configurational space of the translational and rotational degrees of freedom of the molecule, combining steepest-descent and Newton-Raphson steps which embed information on the average sizes of the potential energy wells obtained from prior inspection of the liquid structure. We perform a molecular dynamics simulation of a liquid water shell which demonstrates that the method enables fast and energy-controlled water molecule insertion in aqueous environments. The algorithm finds low energy configurations of incoming water molecules around three orders of magnitude faster than direct random insertion. This method is an important step towards dynamic simulations of open systems and it may also prove useful for energy-biased ensemble average calculations of the chemical potential.

preprint1998arXiv

Finite-difference methods for simulation models incorporating non-conservative forces

We discuss algorithms applicable to the numerical solution of second-order ordinary differential equations by finite-differences. We make particular reference to the solution of the dissipative particle dynamics fluid model, and present extensive results comparing one of the algorithms discussed with the standard method of solution. These results show the successful modeling of phase separation and surface tension in a binary immiscible fluid mixture.

preprint1998arXiv

Inverse Chapman-Enskog Derivation of the Thermohydrodynamic Lattice-BGK Model for the Ideal Gas

A thermohydrodynamic lattice-BGK model for the ideal gas was derived by Alexander et al. in 1993, and generalized by McNamara et al. in the same year. In these works, particular forms for the equilibrium distribution function and the transport coefficients were posited and shown to work, thereby establishing the sufficiency of the model. In this paper, we rederive the model from a minimal set of assumptions, and thereby show that the forms assumed for the shear and bulk viscosities are also necessary, but that the form assumed for the thermal conductivity is not. We derive the most general form allowable for the thermal conductivity, and the concomitant generalization of the equilibrium distribution. In this way, we show that it is possible to achieve variable (albeit density-dependent) Prandtl number even within a single-relaxation-time lattice-BGK model. We accomplish this by demanding analyticity of the third moments and traces of the fourth moments of the equilibrium distribution function. The method of derivation demonstrates that certain undesirable features of the model -- such as the unphysical dependence of the viscosity coefficients on temperature -- cannot be corrected within the scope of lattice-BGK models with constant relaxation time.