Source author record

Benjamin D. Goddard

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

ResearcherUnclaimed source record

cond-mat.stat-mech math-ph math.MP cond-mat.soft math.NA cond-mat.mes-hall cs.CY math.AP Numerical Analysis physics.chem-ph physics.flu-dyn physics.soc-ph quant-ph Social and Information Networks

Catalog footprint

What is connected

9works

14topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Improving tobacco social contagion models using agent-based simulations on networks

Over the years, population-level tobacco control policies have considerably reduced smoking prevalence worldwide. However, the rate of decline of smoking prevalence is slowing down. Therefore, there is a need for models that capture the full complexity of the smoking epidemic. These models can then be used as test-beds to develop new policies to limit the spread of smoking. Current models of smoking dynamics mainly use ordinary differential equation (ODE) models, where studying the effect of an individual's contact network is challenging. They also do not consider all the interactions between individuals that can lead to changes in smoking behaviour, implying that they do not consider valuable information on the spread of smoking behaviour. In this context, we develop an agent-based model (ABM), calibrate and then validate it on historical trends observed in the US and UK. Our ABM considers spontaneous terms, interactions between agents, and the agent's contact network. To explore the effect of the underlying network on smoking dynamics, we test the ABM on six different networks, both synthetic and real-world. In addition, we also compare the ABM with an ODE model. Our results suggest that the dynamics from the ODE model are similar to the ABM only when the network structure is fully connected (FC). The FC network performs poorly in replicating the empirical trends in the data, while the real-world network best replicates it amongst the six networks. Further, when information on the real-world network is unavailable, our ABM on Lancichinetti-Fortunato-Radicchi benchmark networks (or networks with a similar average degree as the real-world network) can be used to model smoking behaviour. These results suggest that networks are essential for modelling smoking behaviour and that our ABM can be used to develop network-based intervention strategies and policies for tobacco control.

preprint2022arXiv

MultiShape: A Spectral Element Method, with Applications to Dynamic Density Functional Theory and PDE-Constrained Optimization

A numerical framework is developed to solve various types of PDEs on complicated domains, including steady and time-dependent, non-linear and non-local PDEs, with different boundary conditions that can also include non-linear and non-local terms. This numerical framework, called MultiShape, is a class in Matlab, and the software is open source. We demonstrate that MultiShape is compatible with other numerical methods, such as differential--algebraic equation solvers and optimization algorithms. The numerical implementation is designed to be user-friendly, with most of the set-up and computations done automatically by MultiShape and with intuitive operator definition, notation, and user-interface. Validation tests are presented, before we introduce three examples motivated by applications in Dynamic Density Functional Theory and PDE-constrained optimization, illustrating the versatility of the method.

preprint2016arXiv

Dynamical density functional theory for orientable colloids including inertia and hydrodynamic interactions

Over the last few decades, classical density-functional theory (DFT) and its dynamic extensions (DDFTs) have become powerful tools in the study of colloidal fluids. Recently, previous DDFTs for spherically-symmetric particles have been generalised to take into account both inertia and hydrodynamic interactions, two effects which strongly influence non-equilibrium properties. The present work further generalises this framework to systems of anisotropic particles. Starting from the Liouville equation and utilising Zwanzig's projection-operator techniques, we derive the kinetic equation for the Brownian particle distribution function, and by averaging over all but one particle, a DDFT equation is obtained. Whilst this equation has some similarities with DDFTs for spherically-symmetric colloids, it involves a translational-rotational coupling which affects the diffusivity of the (asymmetric) particles. We further show that, in the overdamped (high friction) limit, the DDFT is considerably simplified and is in agreement with a previous DDFT for colloids with arbitrary shape particles.

preprint2016arXiv

Wave packet dynamics in the optimal superadiabatic approximation

We explain the concept of superadiabatic approximations and show how in the context of the Born- Oppenheimer approximation they lead to an explicit formula that can be used to predict transitions at avoided crossings. Based on this formula, we present a simple method for computing wave packet dynamics across avoided crossings. Only knowledge of the adiabatic electronic energy levels near the avoided crossing is required for the computation. In particular, this means that no diabatization procedure is necessary, the adiabatic energy levels can be computed on the fly, and they only need to be computed to higher accuracy when an avoided crossing is detected. We test the quality of our method on the paradigmatic example of photo-dissociation of NaI, finding very good agreement with results of exact wave packet calculations.

preprint2014arXiv

Fluid structure in the immediate vicinity of an equilibrium three-phase contact line and assessment of disjoining pressure models using density functional theory

We examine the nanoscale behavior of an equilibrium three-phase contact line in the presence of long-ranged intermolecular forces by employing a statistical mechanics of fluids approach, namely density functional theory (DFT) together with fundamental measure theory (FMT). This enables us to evaluate the predictive quality of effective Hamiltonian models in the vicinity of the contact line. In particular, we compare the results for mean field effective Hamiltonians with disjoining pressures defined through (I) the adsorption isotherm for a planar liquid film, and (II) the normal force balance at the contact line. We find that the height profile obtained using (I) shows good agreement with the adsorption film thickness of the DFT-FMT equilibrium density profile in terms of maximal curvature and the behavior at large film heights. In contrast, we observe that while the height profile obtained by using (II) satisfies basic sum rules, it shows little agreement with the adsorption film thickness of the DFT results. The results are verified for contact angles of 20, 40 and 60 degrees.

preprint2012arXiv

Generalized dynamical density functional theory for classical fluids and the significance of inertia and hydrodynamic interactions

We study the dynamics of a colloidal fluid including inertia and hydrodynamic interactions, two effects which strongly influence the non-equilibrium properties of the system. We derive a general dynamical density functional theory (DDFT) which shows very good agreement with full Langevin dynamics. In suitable limits, we recover existing DDFTs and a Navier-Stokes-like equation with additional non-local terms.

preprint2012arXiv

The overdamped limit of dynamic density functional theory: Rigorous results

Consider the overdamped limit for a system of interacting particles in the presence of hydrodynamic interactions. For two-body hydrodynamic interactions and one- and two-body potentials, a Smoluchowski-type evolution equation is rigorously derived for the one-particle distribution function. This new equation includes a novel definition of the diffusion tensor. A comparison with existing formulations of dynamic density functional theory is also made.

preprint2010arXiv

Non-adiabatic transitions through tilted avoided crossings

We investigate the transition of a quantum wave-packet through a one-dimensional avoided crossing of molecular energy levels when the energy levels at the crossing point are tilted. Using superadiabatic representations, and an approximation of the dynamics near the crossing region, we obtain an explicit formula for the transition wavefunction. Our results agree extremely well with high precision ab-initio calculations.

preprint2009arXiv

Superadiabatic transition histories in quantum molecular dynamics

We study the dynamics of a molecule's nuclear wave-function near an avoided crossing of two electronic energy levels, for one nuclear degree of freedom. We derive the general form of the Schroedinger equation in the n-th superadiabatic representation for all n, and give some partial results about the asymptotics for large n. Using these results, we obtain closed formulas for the time development of the component of the wave function in an initially unoccupied energy subspace, when a wave packet crosses the transition region. In the optimal superadiabatic representation, which we define, this component builds up monontonically. Finally, we give an explicit formula for the transition wave function away from the crossing, which is in excellent agreement with high precision numerical calculations.