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

Matthew Dobson

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

ResearcherUnclaimed source record
math.NAmath-phmath.MPastro-ph.HEastro-ph.SRmath.APmath.PRmath.SPphysics.comp-ph

Catalog footprint

What is connected

15works
9topics
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

15 published item(s)

preprint2020arXiv

SN 2019ehk: A Double-Peaked Ca-rich Transient with Luminous X-ray Emission and Shock-Ionized Spectral Features

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$α$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On assessing the accuracy of defect free energy computations

We develop a rigorous error analysis for coarse-graining of defect-formation free energy. For a one-dimensional constrained atomistic system, we establish the thermodynamic limit of the defect-formation free energy and obtain explicitly the rate of convergence. We then construct a sequence of coarse-grained energies with the same rate but significantly reduced computational cost. We illustrate our analytical results through explicit computations for the case of harmonic potentials and through numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cell List Algorithms for Nonequilibrium Molecular Dynamics

We present two modifications of the standard cell list algorithm for nonequilibrium molecular dynamics simulations of homogeneous, linear flows. When such a flow is modeled with periodic boundary conditions, the simulation box deforms with the flow, and recent progress has been made developing boundary conditions suitable for general 3D flows of this type. For the typical case of short-ranged, pairwise interactions, the cell list algorithm reduces computational complexity of the force computation from O($N^2$) to O($N$), where $N$ is the total number of particles in the simulation box. The new versions of the cell list algorithm handle the dynamic, deforming simulation geometry. We include a comparison of the complexity and efficiency of the two proposed modifications of the standard algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Periodic Boundary Conditions for Long-time Nonequilibrium Molecular Dynamics Simulations of Incompressible Flows

This work presents a generalization of the Kraynik-Reinelt (KR) boundary conditions for nonequilibrium molecular dynamics simulations. In the simulation of steady, homogeneous flows with periodic boundary conditions, the simulation box moves with the flow, and it is possible for particle replicas to become arbitrarily close, causing a breakdown in the simulation. The KR boundary conditions avoid this problem for planar elongational flow and general planar mixed flow [J. Chem. Phys 133, 14116 (2010)] through careful choice of the initial simulation box and by periodically remapping the simulation box in a way that conserves replica locations. In this work, the ideas are extended to a large class of three dimensional flows by using multiple remappings for the simulation box. The simulation box geometry is no longer time-periodic (which was shown to be impossible for uniaxial and biaxial stretching flows in the original work by Kraynik and Reinelt [Int. J. Multiphase Flow 18, 1045 (1992)]). The presented algorithm applies to all flows with nondefective flow matrices, and in particular, to uniaxial and biaxial flows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Derivation of Langevin Dynamics in a Nonzero Background Flow Field

We propose a derivation of a nonequilibrium Langevin dynamics for a large particle immersed in a background flow field. A single large particle is placed in an ideal gas heat bath composed of point particles that are distributed consistently with the background flow field and that interact with the large particle through elastic collisions. In the limit of small bath atom mass, the large particle dynamics converges in law to a stochastic dynamics. This derivation follows the ideas of [D. Dürr, S. Goldstein, and J. L. Lebowitz, 1981 and 1983; P. Calderoni, D. Dürr, and S. Kusuoka, 1989] and provides extensions to handle the nonzero background flow. The derived nonequilibrium Langevin dynamics is similar to the dynamics in [M. McPhie, et al., 2001]. Some numerical experiments illustrate the use of the obtained dynamic to simulate homogeneous liquid materials under flow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

There is no pointwise consistent quasicontinuum energy

Much work has gone into the construction of quasicontinuum energies that reduce the coupling error along the interface between atomistic and continuum regions. The largest consistency errors are typically pointwise $O(\frac{1}{\eps})$ errors, and in some cases this has been reduced to pointwise O(1) errors. In this paper we show that one cannot create a coupling method using a finite-range coupling interface that has o(1)-consistency in the interface, and we use this to give an upper bound on the order of convergence in discrete $w^{1,p}$-norms in 1D.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Accuracy of Quasicontinuum Approximations Near Instabilities

The formation and motion of lattice defects such as cracks, dislocations, or grain boundaries, occurs when the lattice configuration loses stability, that is, when an eigenvalue of the Hessian of the lattice energy functional becomes negative. When the atomistic energy is approximated by a hybrid energy that couples atomistic and continuum models, the accuracy of the approximation can only be guaranteed near deformations where both the atomistic energy as well as the hybrid energy are stable. We propose, therefore, that it is essential for the evaluation of the predictive capability of atomistic-to-continuum coupling methods near instabilities that a theoretical analysis be performed, at least for some representative model problems, that determines whether the hybrid energies remain stable {\em up to the onset of instability of the atomistic energy}. We formulate a one-dimensional model problem with nearest and next-nearest neighbor interactions and use rigorous analysis, asymptotic methods, and numerical experiments to obtain such sharp stability estimates for the basic conservative quasicontinuum (QC) approximations. Our results show that the consistent quasi-nonlocal QC approximation correctly reproduces the stability of the atomistic system, whereas the inconsistent energy-based QC approximation incorrectly predicts instability at a significantly reduced applied load that we describe by an analytic criterion in terms of the derivatives of the atomistic potential.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Iterative Methods for the Force-based Quasicontinuum Approximation

Force-based atomistic-continuum hybrid methods are the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. For this reason, and due to their algorithmic simplicity, force-based coupling methods have become a popular class of atomistic-continuum hybrid models as well as other types of multiphysics models. However, the recently discovered unusual stability properties of the linearized force-based quasicontinuum (QCF) approximation, especially its indefiniteness, present a challenge to the development of efficient and reliable iterative methods. We present analytic and computational results for the generalized minimal residual (GMRES) solution of the linearized QCF equilibrium equations. We show that the GMRES method accurately reproduces the stability of the force-based approximation and conclude that an appropriately preconditioned GMRES method results in a reliable and efficient solution method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Symplectic schemes for highly oscillatory Hamiltonian systems: the homogenization approach beyond the constant frequency case

We follow up on our previous works which presented a possible approach for deriving symplectic schemes for a certain class of highly oscillatory Hamiltonian systems. The approach considers the Hamilton-Jacobi form of the equations of motion, formally homogenizes it and infers an appropriate symplectic integrator for the original system. In our previous work, the case of a system exhibiting a single constant fast frequency was considered. The present work successfully extends the approach to systems that have either one varying fast frequency or several constant frequencies. Some related issues are also examined.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Spectrum of the Force-Based Quasicontinuum Operator for a Homogeneous Periodic Chain

We show under general conditions that the linearized force-based quasicontinuum (QCF) operator has a positive spectrum, which is identical to the spectrum of the quasinonlocal quasicontinuum (QNL) operator in the case of second-neighbour interactions. Moreover, we establish a bound on the condition number of a matrix of eigenvectors that is uniform in the number of atoms and the size of the atomistic region. These results establish the validity of and improve upon recent conjectures ([arXiv:0907.3861, Conjecture 2] and [arXiv:0910.2013, Conjecture 8]) which were based on numerical experiments. As immediate consequences of our results we obtain rigorous estimates for convergence rates of (preconditioned) GMRES algorithms, as well as a new stability estimate for the QCF method.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An Optimal Order Error Analysis of the One-Dimensional Quasicontinuum Approximation

We derive a model problem for quasicontinuum approximations that allows a simple, yet insightful, analysis of the optimal-order convergence rate in the continuum limit for both the energy-based quasicontinuum approximation and the quasi-nonlocal quasicontinuum approximation. The optimal-order error estimates for the quasi-nonlocal quasicontinuum approximation are given for all strains up to the continuum limit strain for fracture. The analysis is based on an explicit treatment of the coupling error at the atomistic to continuum interface, combined with an analysis of the error due to atomistic and continuum schemes using the stability of the quasicontinuum approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Sharp Stability Estimates for the Force-based Quasicontinuum Method

A sharp stability analysis of atomistic-to-continuum coupling methods is essential for evaluating their capabilities for predicting the formation and motion of lattice defects. We formulate a simple one-dimensional model problem and give a detailed analysis of the stability of the force-based quasicontinuum (QCF) method. The focus of the analysis is the question whether the QCF method is able to predict a critical load at which fracture occurs. Numerical experiments show that the spectrum of a linearized QCF operator is identical to the spectrum of a linearized energy-based quasi-nonlocal quasicontinuum operator (QNL), which we know from our previous analyses to be positive below the critical load. However, the QCF operator is non-normal and it turns out that it is not generally positive definite, even when all of its eigenvalues are positive. Using a combination of rigorous analysis and numerical experiments, we investigate in detail for which choices of "function spaces" the QCF operator is stable, uniformly in the size of the atomistic system. Force-based multi-physics coupling methods are popular techniques to circumvent the difficulties faced in formulating consistent energy-based coupling pproaches. Even though the QCF method is possibly the simplest coupling method of this kind, we anticipate that many of our observations apply more generally.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stability, Instability, and Error of the Force-based Quasicontinuum Approximation

Due to their algorithmic simplicity and high accuracy, force-based model coupling techniques are an exciting development in computational physics. For example, the force-based quasicontinuum approximation is the only known pointwise consistent quasicontinuum approximation for coupling a general atomistic model with a finite element continuum model. In this paper, we analyze the stability of the force-based quasicontinuum approximation. We then use our stability result to obtain an optimal order error analysis of this coupling method that provides theoretical justification for the high accuracy of the force-based quasicontinuum approximation -- the computational efficiency of continuum modeling can be utilized without the loss of significant accuracy if defects are captured in the atomistic region. The main challenge we need to overcome is the fact (which we prove) that the linearized quasicontinuum operator is typically not positive definite. Moreover, we prove that no uniform inf-sup stability condition holds for discrete versions of the $W^{1,p}$-$W^{1,q}$ "duality pairing" with $1/p+1/q=1$, if $1 \leq p < \infty$. We must therefore derive an inf-sup stability condition for a discrete version of the $W^{1,\infty}$-$W^{1,1}$ "duality pairing" which then leads to optimal order error estimates in a discrete $W^{1,\infty}$-norm.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

An Analysis of the Effect of Ghost Force Oscillation on Quasicontinuum Error

The atomistic to continuum interface for quasicontinuum energies exhibits nonzero forces under uniform strain that have been called ghost forces. In this paper, we prove for a linearization of a one-dimensional quasicontinuum energy around a uniform strain that the effect of the ghost forces on the displacement nearly cancels and has a small effect on the error away from the interface. We give optimal order error estimates that show that the quasicontinuum displacement converges to the atomistic displacement at the optimal rate O($h$) in the discrete $\ell^\infty$ norm and O($h^{1/p}$) in the $w^{1,p}$ norm for $1 \leq p < \infty.$ where $h$ is the interatomic spacing. We also give a proof that the error in the displacement gradient decays away from the interface to O($h$) at distance O($h|\log h|$) in the atomistic region and distance O($h$) in the continuum region. E, Ming, and Yang previously gave a counterexample to convergence in the $w^{1,\infty}$ norm for a harmonic interatomic potential. Our work gives an explicit and simplified form for the decay of the effect of the atomistic to continuum coupling error in terms of a general underlying interatomic potential and gives the estimates described above in the discrete $\ell^\infty$ and $w^{1,p}$ norms.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Analysis of a Force-Based Quasicontinuum Approximation

We analyze a force-based quasicontinuum approximation to a one-dimensional system of atoms that interact by a classical atomistic potential. This force-based quasicontinuum approximation is derived as the modification of an energy-based quasicontinuum approximation by the addition of nonconservative forces to correct nonphysical ``ghost'' forces that occur in the atomistic to continuum interface. We prove that the force-based quasicontinuum equations have a unique solution under suitable restrictions on the loads. For Lennard-Jones next-nearest-neighbor interactions, we show that unique solutions exist for loads in a symmetric region extending nearly to the tensile limit. We give an analysis of the convergence of the ghost force iteration method to solve the equilibrium equations for the force-based quasicontinuum approximation. We show that the ghost force iteration is a contraction and give an analysis for its convergence rate.

0 citations0 reviewsNo code linkedOpen access