Researcher profile

Matthew Dobson

Matthew Dobson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 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.

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.

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.

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.

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.

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.

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.

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 &#34;function spaces&#34; 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.

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.

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&#39;&#39; 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.