# Computational Physics (physics.comp-ph)

• The finite-difference time-domain (FDTD) method is a well established method for solving the time evolution of Maxwell's equations. Unfortunately the scheme introduces numerical dispersion and therefore phase and group velocities which deviate from the correct values. The solution to Maxwell's equations in more than one dimension results in non-physical predictions such as numerical dispersion or numerical Cherenkov radiation emitted by a relativistic electron beam propagating in vacuum. Improved solvers, which keep the staggered Yee-type grid for electric and magnetic fields, generally modify the spatial derivative operator in the Maxwell-Faraday equation by increasing the computational stencil. These modified solvers can be characterized by different sets of coefficients, leading to different dispersion properties. In this work we introduce a norm function to rewrite the choice of coefficients into a minimization problem. We solve this problem numerically and show that the minimization procedure leads to phase and group velocities that are considerably closer to $c$ as compared to schemes with manually set coefficients available in the literature. Depending on a specific problem at hand (e.g. electron beam propagation in plasma, high-order harmonic generation from plasma surfaces, etc), the norm function can be chosen accordingly, for example, to minimize the numerical dispersion in a certain given propagation direction. Particle-in-cell simulations of an electron beam propagating in vacuum using our solver are provided.
• Peridynamics (PD) represents a new approach for modelling fracture mechanics, where a continuum domain is modelled through particles connected via physical bonds. This formulation allows us to model crack initiation, propagation, branching and coalescence without special assumptions. Up to date, anisotropic materials were modelled in the PD framework as different isotropic materials (for instance, fibre and matrix of a composite laminate), where the stiffness of the bond depends on its orientation. A non-ordinary state-based formulation will enable the modelling of generally anisotropic materials, where the material properties are directly embedded in the formulation. Other material models include rocks, concrete and biomaterials such as bones. In this paper, we implemented this model and validated it for anisotropic composite materials. A composite damage criterion has been employed to model the crack propagation behaviour. Several numerical examples have been used to validate the approach, and compared to other benchmark solution from the finite element method (FEM) and experimental results when available.
• Atomic properties such as field shift constants, magnetic dipole and electric quadrupole hyperfine structure constants, Landé $g_J$ factors, and electric quadrupole moments that are described by electronic operators with different ranks and radial behaviors are studied and the role of electron correlation effects in their determination are investigated. We have adopted the Dirac-Hartree-Fock method, the second- and third-order relativistic many-body perturbation theories, and an all-order relativistic many-body method in the coupled-cluster theory framework considering only the linearized terms and also all the non-linearized terms in the singles and doubles with partial triples excitations approximation to carry out these analyses. Variations in the propagation of electron correlation effects with operators having same angular factors but different radial behaviors and with different ranks are highlighted. Corrections from the higher-order relativistic corrections due to the Breit and quantum electrodynamics interactions to all these properties are also estimated. Understanding of trends of electron correlation effects in these properties can be useful to establish accuracies in the theoretical results of different atomic properties and to substantiate validity of an approximated many-body method.

### Recent comments

Chris Ferrie May 01 2017 12:52 UTC

I've posted a public referee report on this paper here: https://nbviewer.jupyter.org/github/csferrie/openreviews/blob/master/arxiv.1703.10743/arxiv.1703.10743.ipynb

Salvatore Mandrà Nov 21 2014 13:50 UTC

The manuscript has been widely revised to focus the reader's attention on the proposed method and its application in presence of local disorder.

Best,

Salvatore, Gian Giacomo and Alán

Salvatore Mandrà Aug 01 2014 19:11 UTC

Thanks Dr. Hastings for your comment. It is true that the transverse field Ising model does not satisfy the requirements to apply our method with an exponential reduction. Indeed, the opposite would be quite impressive since the random Ising model is a NP-Hard problem and we ourselves would be prett

...(continued)
Matt Hastings Aug 01 2014 16:43 UTC

The "quite general" conditions seem not to include the transverse field Ising model, the subject of most of the intensive numerical work previously. Incidentally, the terms "Lanczos" and "Krylov subspace" might be helpful.

Steve Flammia Jan 10 2014 08:01 UTC

This paper is an important first step towards making numerical and computational studies more easily checkable and reproducible. Very interesting stuff.

Noon van der Silk Jan 10 2014 16:26 UTC

Interestingly, another thematically-similar paper popped up today - https://scirate3.herokuapp.com/1401.2134 (posting the link here as I assume not many people are subbed to cs.DL)