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.